MATLAB上机实验(答案).docVIP

MATLAB上机实验(答案)

00

00.0914+0.4586i00

000.0914-0.4586i0

0000.2275

（2）给出一系列的a值，采用函数画一组椭圆。

程序：

a=0.5:0.5:4.5;%a的绝对值不能大于5

t=[0:pi/50:2*pi];%用参数t表示椭圆方程

X=cos(t)*a;

Y=sin(t)*sqrt(25-a.^2);

plot(X,Y)

结果：

（3）X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2],

（a）写出计算其负元素个数的程序。

程序：

X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2];

L=X0;

A=sum(L)

结果：

A=

5

（b）写出一段程序，使其能够找出向量x中的最大、最小元素。（不能使用min和max命令）

程序：

X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2];

xmin=999;xmax=-999;

fori=1:length(X)

ifxminX(i)

xmin=X(i);

end

ifxmaxX(i)

xmax=X(i);

end

end

[xmin,xmax]

结果：

ans=

-69

方波函数为f(t)=,利用,用MATLAB编程和绘图说明方波是奇次谐波的叠加。

程序：

k=1000;%k值可以改动

x=-pi:0.0001:pi;

y=sin(x);

forn=1:k

y=y+sin((2*n+1)*x)/(2*n+1);

end

plot(x,y)

结果：

MATLAB工具软件实验(2)

应用simulink实现下列系统的仿真，并试试改变信号源后的结果。

略。

（2）利用simulink解二阶微分方程

其中u(t)为单位阶跃函数

框图：

结果：

（3）求矩阵的行列式、逆和特征值。

程序：

symsa11a12a21a22;

A=[a11,a12;a21,a22]

B=det(A),C=inv(A),D=eig(A)

结果：

A=

[a11,a12]

[a21,a22]

B=

a11*a22-a12*a21

C=

[a22/(a11*a22-a12*a21),-a12/(a11*a22-a12*a21)]

[-a21/(a11*a22-a12*a21),a11/(a11*a22-a12*a21)]

D=

1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)

1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)

（4）简化

程序：

symsx;

f=(x^(-3)+6*x^(-2)+12*x^(-1)+8)^(1/3);

g1=simple(f)

g2=simple(g1)

结果：

g1=

(2*x+1)/x

g2=

2+1/x

（5）计算积分

程序：

symsx;

symsareal;%将a定义为实符号变量

f=2^(-(a*x)^2)

int(f,x,2,inf)

结果：

f=

2^(-a^2*x^2)

ans=

-1/2*pi^(1/2)*(-signum(a)+erf(2*a*log(2)^(1/2)))/a/log(2)^(1/2)

MATLAB工具软件实验(3)

（1）对函数，用梯形法、辛普森法计算积分：，并讨论分割区间数n对误差的影响。其中S的精确值取为4.006994。

程序：

formatlong

n=10:10:100;

forii=1:10

x=0:2/n(ii):2;

f=sqrt(1+exp(x));

s=trapz(x,f)

err(ii)=4.006994-s;

end

figure(1)

plot(n,err)

f1=inline(sqrt(1+exp(x)));

quad(f