Matlab实验及答案 下载本文

2.求下列函数的符号积分 (1) y=1/sqrt(1-x^2);

3.求两个函数之间的操作 求和

(1) sin(x)+cos(x)

乘积

(1) exp(-x)*sin(x)

商 (2)

(1) sin(x)/cos(x);

求复合函数 (1) y=exp(u)

u=sin(x)

(10) 总结who,whos,clc,clear命令

三、设计提示

1.初次接触Matlab应该注意函数表达式的文本式描述。

2.在使用图形函数计算器funtool时,注意观察1号和2号窗口中函数的图形。 3. help帮助命令的使用 3. 向量与矩阵的创建

实验二:数组运算及线型方程组的求解

二、实验内容

1.“:”号的用法。用“:”号生成行向量a=[1 2 3 4 5 6 7 8 9 10]、b=[5 3 1 -1 -3 -5];用线性等分命令linspace重新生成上述的a和b向量。另,在100和10000之间用对数等分命令logspace生成10维的向量c。ak=logspace(2,4,10) linspace(1,10,10) linspace(5,-5,6) 2. 已知多项式a(x)=x2+2x+3,b(x)=4x2+5x+6,求a,b的积并微分。 >> a=[1,2,3];b=[4,5,6];polyder(a,b) ans =

16 39 56 27 >> poly2str(ans,'x') ans =

16 x^3 + 39 x^2 + 56 x + 27

3.生成下列矩阵,取出方框内的数组元素

a(2,2:3) a(2:4,4) a(4:5,1:3)

4. 生成一个9×9维的魔方矩阵,提取其中心的3×3维子矩阵M,利用sum函数检验其各行和各列的和是否相等。并且实现上述中心矩阵左旋90°或右旋90°,左右翻转,上下翻转 a=magic(9) >> b=a(4:6,4:6) b =

20 31 42 30 41 52 40 51 62 >> sum(b,1) ans =

90 123 156 >> sum(b,2) ans = 93 123 153 rot90(b) rot90(b,-1)

fliplr(b) flipud(b)

5.已知a=[1 2 3;4 5 6;7 8 0],求其特征多项式并求其根、特征值和特征多项式 >> a=[1 2 3;4 5 6;7 8 0] a =

1 2 3 4 5 6 7 8 0 >> poly(a) ans =

1.0000 -6.0000 -72.0000 -27.0000 >> poly2sym(ans) ans =

x^3-6*x^2-72*x-7599824371187741/281474976710656 >> [d v]=eig(a) d =

-0.2998 -0.7471 -0.2763 -0.7075 0.6582 -0.3884 -0.6400 -0.0931 0.8791 v =

12.1229 0 0 0 -0.3884 0 0 0 -5.7345

6. 计算二重不定积分 >> syms x y

>> f=int(int(x*exp(-x*y),x),y)