Matlab实验报告.docx
- 文档编号:4993285
- 上传时间:2022-12-12
- 格式:DOCX
- 页数:21
- 大小:384.51KB
Matlab实验报告.docx
《Matlab实验报告.docx》由会员分享,可在线阅读,更多相关《Matlab实验报告.docx(21页珍藏版)》请在冰豆网上搜索。
Matlab实验报告
MATLAB
学院:
信息科学与技术学院
专业:
信息工程
姓名:
黄成林
学号:
2010********
试验一MATLAB运算基础
1.1
程序代码:
t=0:
0.5:
2.5
y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*((t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2)&(t<3))
结果:
1.3
程序代码:
ch='ABcdefGHd',k=find(ch>='A'&ch<='Z'),ch(k)=[]
结果:
1.5
程序代码:
a11=input('a11='),a12=input('a12='),
a21=input('a21='),a22=input('a22='),
A=[a11,a12;a21,a22],
DA=det(A),IA=inv(A),EA=eig(A)
结果:
实验二M函数与M文件的编写与应用
2.1
程序代码:
n=100;a=ones(1,n);
fori=3:
n
a(i)=a(i-1)+a(i-2);
ifa(i)>10000
a(i),
break;
end;
end,i
结果:
2.3
程序代码:
functiony=ff(t)
y1=exp(-t/3);
y2=exp(-t/3).*sin(3*t);
y=[y1;y2]
结果:
实验三选择与循环结构程序设计
3.1
程序代码:
forn=1:
100
f(n)=1./(2*n-1)
y=sum(f)
ify>=3
my=y-f(n)
mn=n-1
break
end
end
my
Mn
结果:
3.3
程序代码:
per=input('输入成绩:
')
switchfloor(per/10)
case9
grade='A'
case8
grade='B'
case7
grade='C'
case6
grade='D'
casenum2cell(0:
5)
grade='E'
otherwise
ifper==100
grade='A'
elsegrade='error'
end
end
结果:
实验四绘图函数的应用
4.1
程序代码:
t=0:
0.05:
2*pi;
plot(t,sin(t),'r-.o',t,cos(t),'m-s')
legend('sinx','cosx',4)
axis([02*pi-22])
结果:
4.3
程序代码:
x=0:
0.1:
100;
y=log10(x);
subplot(311),loglog(x,y)
grid
title('logloggraph')
subplot(312),semilogx(x,y)
grid
title('semilogxgraph')
subplot(313),semilogy(x,y)
grid
title('semilogygraph')
结果:
4.5
程序代码:
x=[-2:
0.02:
2];
y=1.1*sign(x).*(abs(x)>1.1)+x.*(abs(x)<=1.1);
plot(x,y)
结果:
实验五数值工具箱与符号工具箱的应用
5.1
程序代码:
h=6:
2:
18;
t1=[18202225302824];
t2=[15192428343230];
h1=6.5:
2:
17.5;
t11=interp1(h,t1,h1,'spline')
t22=interp1(h,t2,h1,'spline')
结果:
5.3
程序代码:
clear
fxy=sym('log(x+exp(y))/sqrt(x^2+y^2)')
result=limit(limit(fxy,'x',1),'y',0)
结果:
5.5
程序代码:
clear
symsaxty1y2
y1=sqrt(a+t)
y2=int(y1,t,x,x^2);
diff(y2,x)
结果:
5.7
程序代码:
clear
symsxyz
f=x+2*y-2*sqrt(x*y*z);
fx=diff(f,x);fy=diff(f,y);fz=diff(f,z);
zx=-fx/fz
zy=-fy/fz
结果:
实验六Simulink仿真应用
例5.6
(1)
(2)
程序代码:
%exm5_6.m
clf;holdon
xx=[-2,1;-1,1;0,1;1,1;1,0;1,-1;1,-2];
nxx=size(xx,1);
fork=1:
nxx
opts=simset('initialstate',[xx(k,1),xx(k,2)]);
[t,x,y]=sim('exm5_5',10,opts);
plot(x(:
1),x(:
2));
end
xlabel('x1');ylabel('x2'),grid,holdoff
(3)
结果:
(4)
程序代码:
function[DX1,DX2,DP]=portraitzzy(x1,x2,h)
%PORTRAITZZY
%x1,x2
%h
%DX1,DX2
%DP
opts=simset('solver','ode5','fixedstep',h);
n=length(x1);
X1=zeros(n,n);X2=X1;
forii=1:
n;
forjj=1:
n;
opts=simset(opts,'initialstate',[x1(ii),x2(jj)]);
[t,x,y]=sim('exm5_5',h,opts);
dx1=x(2,1)-x1(ii);
dx2=x(2,2)-x2(jj);
L=sqrt(dx1^2+dx2^2);
Z(jj,ii)=L
ifL>1.e-10
DX1(jj,ii)=dx1/L;DX2(jj,ii)=dx2/L;
end
end
end
DP=Z/h;
(5)
程序代码:
h=0.01;
x1=-2.5:
0.25:
2.5;x2=x1;
k=0.15;
[X1,X2,Z]=portraitzzy(x1,x2,h);
quiver(x1,x2,k*X1,k*X2,0)
xlabel('x1'),ylabel('x2')
结果:
(6)
程序代码:
surfc(x1,x2,Z),view([18,32]),xlabel('x1'),ylabel('x2')
结果:
(7)
程序代码:
xa=trim('exm5_5',[-1,-2]')
xb=trim('exm5_5',[1,2]')
结果:
(8)
程序代码:
Axa=linmod('exm5_5',xa);eig_Axa=(eig(Axa.a))'
Axb=linmod('exm5_5',xb);eig_Axb=(eig(Axb.a))'
结果:
实验七线性系统分析与设计
7.1
程序代码:
num=[0.00010.02181.04369.3599];
den=[0.00060.02680.63656.2711];
sys=tf(num,den);sys1=ss(sys)
sys2=zpk(sys)
结果:
7.3
程序代码:
clc
clear
closeall
num=[00281282];den=[151010510];
sys=tf(num,den)
bode(sys)
figure
nichols(sys),
figure
nyquist(sys)
结果:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 实验 报告