trapz

概括：这道题是王冠吵同学的课后练习题，主要是关于trapz，指导老师为魏老师。核苷酸抗代谢物，是一些碱基、氨基酸或叶酸等的类似物，它们以多种方式干扰或阻断核苷酸的合成代谢，从而进一步阻止核酸及蛋白质的生物合成，这些代谢物具有抗肿瘤作用。

题目：trapz

采用梯形法求积分

下面两种用法等效

h=pi/200;

t=0:h:pi/2;

y=sin(t)+cos(2*t);

z1=trapz(t,y)

z2=trapz(y)*h

z1 =

1.0000

z2 =

1.0000

举一反三

例1： Matlab中getframe函数的用法比如这句F=getframe(2,[85,58,20,20])；这句中的各个数字表示什么意思,忘知道的大侠解答下,感激不尽!

matlab help:

F = getframe(h,rect) specifies

a rectangular area from which to copy the pixmap.rect is

relative to the lower left corner of the figure or axes h,

in pixel units.rect is a four-element vector in the form [left

bottom width height],where width and height define

the dimensions of the rectangle.

本函数是要从图形对象中截取一块出来.其中h是图像句柄,rect是一个四元素向量,代表[a b c d],其中a表示截取区域到图像左端的像素数,b表示截取区域到图像底端的像素数,c表示截取区域水平像素数,d表示截取区域垂直向像素数.

得到的F是一个结构体(struct),包含两个元素：cdata和colormap,其中cdata是截取区域的数据,是一个三维数组,其中第三维的长度是3,即截取到的是区域中RGB的值.colormap是调色板,这个没啥意思.

F = getframe(2,[85,58,20,20]),显然 F.cdata 是一个 20*20*3的数组.

尝试以下代码：

Z = peaks; surf(Z);

F = getframe(gcf,[150 150 100 100]);

figure,imshow(F.cdata);

例2： matlab中unifrnd函数用法

函数功能：

生成(连续)均匀分布的随机数

使用方法：

R = unifrnd(A,B)

生成被A和B指定上下端点[A,B]的连续均匀分布的随机数组R.

如果A和B是数组,R(i,j)是生成的被A和B对应元素指定连续均匀分布的随机数.

如果N或P是标量,则被扩展为和另一个输入有相同维数的数组.

R = unifrnd(A,B,m,n,...) 或 R = unifrnd(A,B,[m,n,...])

返回m*n*...数组.

如果A和B是标量,R中所有元素是相同分布产生的随机数.

如果A或B是数组,则必须是m*n*...数组.

应用举例 ：

a = 0; b = 1:5;

r1 = unifrnd(a,b)

r1 =

0.8147 1.8116 0.3810 3.6535 3.1618

B = repmat(b,5,1);

R = unifrnd(a,B)

R =

0.0975 0.3152 0.4257 2.6230 3.7887

0.2785 1.9412 1.2653 0.1428 3.7157

0.5469 1.9143 2.7472 3.3965 1.9611

0.9575 0.9708 2.3766 3.7360 3.2774

0.9649 1.6006 2.8785 2.7149 0.8559

r2 = unifrnd(a,b(2),1,5)

r2 =

1.4121 0.0637 0.5538 0.0923 0.1943

例3： matlab中simple函数的用法

simple('s'),s是字符串,如果是对某个字符函数如y=(2+x)/x;化简,直接用simple(y)就得到2/x+1;

有时表达式比较复杂,一次化简结果不太理想,就再次化简,可以用simple(simple(y))

例4： 【matlab中strcat函数的用法filename=strcat('D:\thesis\data\gamit\coordinate\',yr,'\',model,yr,'.org')的意思是啥啊,求说明括号里面的意思.】

先明白strcat函数的定义：

定义

strcat 即 Strings Catenate,横向连接字符串.

语法

combinedStr= strcat(s1,s2,...,sN)

描述

将数组 s1,s2,...,sN 水平地连接成单个字符串,并保存于变量combinedStr中.如果任一参数是元胞数组,那么结果 combinedStr 是一个元胞数组,否则,combinedStr是一个字符数组.

其实它的作用很简单就是将多个字符连接成单个字符串,关键在于这个语句中括号里面的内容,

'D:\thesis\data\gamit\coordinate\' 加了单引号的,就说明单引号里面的内容为：字符串 ,下面的'\' 与'.org' 相同都只是字符串而已,原样输出即可.置于接下来的 yr 和 model 肯定是先前人为定义的变量,里面肯定储存了某些字符串信息的.

那么这个语句最终实现的结果是：

将“ D:\thesis\data\gamit\coordinate\yr的值\model的值.org” 存于变量filename 中

看程序部能光看一句,这样这个语句事实上完成了 一个文件的创建,便于接下去程序中实现 调用路径中的文件数据 或 修改.

ps：如果你学过C语言这个对你来说应该不难

例5： 【MATLAB中feval函数的用法给个简单的例子就可以不要matlab中feval的英文,因为已经看过了,看不太明白】

feval就是把已知的数据或符号带入到一个定义好的函数句柄中,你看看下面的例子

syms t

f=@(x,y) x^2+y^2

k1=feval(f,1,t)

k2=f(1,t)

k3=feval(f,1,1)

k4=f(1,1)

k1 =

t^2 + 1

k2 =

t^2 + 1

k3 =

2

k4 =

2

相关思考练习题：

题1：matlab中的cumtrapz是什么意思？

点拨：trapz是利用梯形法求数值积分 要求输入的是一个数列，或者是两个一样长的对应数列 trapz(Y)或者trapz(X,Y); 就是计算由X，Y两组数据决定的曲线下的积分面积 而当只输入一个变量时，函数默认其为Y变量， 而X变量缺省默认为等间距为1的等差数列，...

题2：matlab trapz的使用

点拨：只有一组数据直接用trapz(data)横坐标默认横坐标等间增量为1， 相当于trapz(1:length(data),data) 如果有横坐标x 那么用trapz(x,data) x和data要一样长，一一对应，x不需要等间距

题3：MATLAB TRAPZ() 用法

点拨：你把Y中x换成X，我试了试没错 >> X = 0:1/160:1; >> Y = exp(10*cos(pi*X)); >> Q = trapz(X,Y) Q = 2.8157e+03

题4：matlab数值积分trapz 出现错误，不知道在哪

点拨：你的程序中f=trapz(x,s)，s是一个数，不是矢量，所以，提示错误。

题5：matlab trapz是派什么用处的

点拨：matlab trapz用处 只有一组数据直接用trapz(data)横坐标默认横坐标等间增量为1， 相当于trapz(1:length(data),data) 如果有横坐标x 那么用trapz(x,data) x和data要一样长，一一对应，x不需要等间距