MATLAB系统仿真实验报告范本全套源码及运行结果

2022-2023学年《MATLAB系统仿真》实验报告（范本）院系：班级：学号：姓名：联系电话：实验19.1.2（2）题目代码和答案①（365-522-70）3>>(365-52*2-70)/3ans=63.6667②area=pi*2.5^2>>area=pi*2.5^2area=19.6350③z1=2*sin(85*pi/180)/(1+exp(2))z1=2*sin(85*pi/180)/(1+exp(2))z1=0.2375④x=[21+2i;-0.455]z2=>>x=[21+2i;-0.455];>>z2=1/2*log(x+(1+x^2)^0.5)z2=0.7182-0.0088i0.7461+0.5434i-0.1473-1.5679i1.1555-0.0016i⑤Knownx=3,y=4,calculatez=>>x=3;>>y=4;>>x^2*y^3/(x-y)^2ans=576⑥Assignthefollowingmatrixtom1m1=[162313;511108;97612;414151];Performthefollowingcommands:>>m1(2,3)>>m1(11)>>m1(:,3)>>m1(2:3,1:3)>>m1(1,4)+m1(2,3)+m1(3,2)+m1(4,1)>>m1=[162313;511108;97612;414151];>>m1(2,3)ans=10>>m1(11)ans=6>>m1(:,3)ans=310615>>m1(2:3,1:3)ans=51110976>>m1(1,4)+m1(2,3)+m1(3,2)+m1(4,1)ans=34⑦Performthefollowingcommands:>>x=0:0.1:6*pi;>>y=5*sin(x);>>plot(x,y)>>x=0:0.1:6*pi;>>y=5*sin(x);>>plot(x,y)（3）Matrixoperation①题目代码和答案a)result1=a'>>a=[123;456];b=[24-1;135];c=[10-2];d=[147;852;360];>>result1='a'result1=ab)result2=a*b不合法illegal>>result2=a*b错误使用*内部矩阵维度必须一致。c)result3=a+b>>result3=a+bresult3=3625811d)result4=b*d>>result4=b*dresult4=312222404913e)result5=[b;c']*d不合法illegal>>result5=[b;c']*d错误使用vertcat串联的矩阵的维度不一致。f)result6=a.*bg)result7=a./b>>result6=a.*bresult6=28-341530>>result7=a./bresult7=0.50000.5000-3.00004.00001.66671.2000h)result8=a.*c不合法illegal>>result8=a.*c错误使用.*矩阵维度必须一致。i)result9=a.\b>>result9=a.\bresult9=2.00002.0000-0.33330.25000.60000.8333j)result10=a.^2>>result10=a.^2result10=149162536k)result11=a^2不合法illegal>>result11=a^2错误使用^输入必须为标量和方阵。要按元素进行POWER计算，请改用POWER(.^)。l)result12=2.^a>>result12=2.^aresult12=248163264②SolvethefollowingequationusingMATLAB.输入代码：>>a=[1110;121-1;2-103;335-6];>>b=[1;8;3;5];>>x_y_z_w=a\bx_y_z_w=9.00003.0000-11.0000-4.0000代码：A=[721-2;9153-2;-2-2115;13213];>>rank(A)ans=4>>det(A)ans=1.2568e+04>>inv(A)ans=0.1744-0.0303-0.01250.0270-0.10500.0789-0.01210.00060.00830.01730.0911-0.03110.0095-0.0185-0.01030.0795>>[v,d]=eig(A)v=0.7629+0.0000i0.0919+0.0640i0.0919-0.0640i-0.0299+0.0000i0.6223+0.0000i0.6087+0.0276i0.6087-0.0276i0.2637+0.0000i0.0807+0.0000i-0.7474+0.0000i-0.7474+0.0000i0.6434+0.0000i-0.1554+0.0000i0.0342-0.2374i0.0342+0.2374i0.7180+0.0000id=4.8554+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i12.6460+1.8333i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i12.6460-1.8333i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i15.8526+0.0000i代码：>>a=20;>>b=-2;>>c=0;>>d=1;>>r1=a>br1=1>>r2=a>b&c>dr2=0>>r3=a==b*(-10)r3=1>>r4=~b|cr4=0实验二1.2.2（1）>>a=[3,4,7,2,9,12];>>r=roots(a)r=-0.8612+1.4377i-0.8612-1.4377i0.6737+1.0159i0.6737-1.0159i-0.9583+0.0000i>>plot(r,'*')代码和答案：>>a=[10000-1];>>r=roots(a)r=-0.8090+0.5878i-0.8090-0.5878i0.3090+0.9511i0.3090-0.9511i1.0000+0.0000i>>plot(r,'*')(3）代码和答案：>>f=@(x)x^3-2*x^2*sin(x)+5*x*cos(x)+1/x;>>x0=[0.5,1.52];>>x=fzero(f,x0)x=1.5117>>f=@(x)x^3-2*x^2*sin(x)+5*x*cos(x)+1/x;>>x1=[1.52,4];>>x=fzero(f,x1)x=2.6095（4）代码和答案：>>p=[0.01.12.12.84.25.06.16.98.19.09.9];>>u=[1011131417182224293439];>>P=polyfit(p,u,3);>>a=P(1);b=P(2);c=P(3);d=P(4);>>t=0:0.1:10;>>u_n=a*t.^3+b*t.^2+c*t+d;>>plot(p,u,'*',t,u_n)(5）代码和答案：>>z=0:0.1:4*pi;>>x=cos(z);>>y=sin(z);>>plot3(x,y,z)代码和答案：[x,y]=meshgrid(-3:0.1:3,-3:0.1:3);

z=-x.^2/10+y.^2/10;

>>

mesh(x,y,z)surf(x,y,z)实验三1.3.2（1）while1n=input('Entern,negativequits:');ifn<=0breakenda=n;whilen>1ifrem(n,2)==0n=n/2;elsen=3*n+1;enda=[a,n];endaend运行结果：>>UntitledEntern,negativequits:12a=1263105168421Entern,negativequits:90a=90451366834175226134020105168421Entern,negativequits:-10显然当输入数字为正数时，该程序会持续运行；当输入数字小于等于0时，该程序终止。（2）①代码：sum=0;result=0;n=input('请输入数字：')fori=1:nsum=sum+1/(i*i);endresult=sqrt(sum*6)结果：>>untitled请输入数字：100n=100result=3.1321>>untitled请输入数字：1000n=1000result=3.1406>>untitled请输入数字：10000n=10000result=3.1415②代码：sum=0;n=1;whilen<=10sum=sum+factorial(n);n=n+1;enddisp(sum)结果：>>untitled4037913(3)代码：score=input('请输入分数：');ifscore>=90disp('excellent');elseifscore>=80&score<90disp('good');elseifscore>=70&score<80disp('mean');elseifscore>=60&score<70disp('passed');elseifscore<60disp('failed');elsedisp('成绩有误')end结果：请输入分数：98excellent>>untitled