实验四 MATLAB编程练习.docVIP

实验五　编程练习一　 一．实验目的和要求 1．熟悉Matlab数据结构、控制流、输入输出实验通过熟悉和掌握M文件式建立与使用方法、以及函数与控制程序流程的语句使用，使学生具备一定的编程和调试能力用循环和非循环方法求1~100000间整数的立方和，并比较两者的时间复杂度。 tic; x=1:100000; y=x.^3; sum(y) toc ans = 2.5001e+019 elapsed_time = 0.0160 循环方式: tic; sum=0; for i=1:100000 sum=sum+i^3; end sum toc sum = 2.5001e+019 elapsed_time = 0.3750 2．利用for循环求1！+2！+3！+ 4！+5！的值sum=0; for i=1:5 x=1; for j=1:i x=x*j; end sum=sum+x; end 3．求[100，200]之间所有能被29整除的整数。for i=100:200 x=mod(i,29); if x==0 sum=i,end end 4．从键盘输入成绩百分制成绩x, 将其转换为等级grade为“满分”、“优秀”、“良好”、“及格”、“不及格”从屏幕输出。转换规则如下：x=100，满分x=90，优秀x=80，良好x=60，及格x60，不及格x=input(x=) if x==100 result=满分, elseif x=90x100 result =优秀, elseif x=80x90 result =良好, elseif x=60x80 result =及格, elseif x60 result =不及格,end 实验六　编程练习二　 一．实验目的和要求 1．掌握用MATLAB解决实际电路和系统问题的基本思路； 2．根据已学的电路知识和信号与系统理论，应用MATLAB求解电路和系统相关未知量。 二．实验器材 计算机、MATLAB软件 三．实验内容和要求 1．在如图1所示的电路中，已知R=3, R=5, R=9, R=2, R=6 R= R=3。（1）已知Us＝48V，求Is，I（2）已知I＝2A，求Us，Is。要求进行理论分析，并将MATLAB计算结果和理论计算结果进行比较。 图1 clear all; r1=3;r2=5;r3=9;r4=2;r5=6;r6=1;r7=1; display(解决问题1); a11=r2+r5;a12=-r2;a13=-r5; a21=-r2;a22=r1+r2+r3+r4;a23=-r4; a31=-r5;a32=-r4;a33=r4+r5+r6+r7; b1=1;b2=0;b3=0; us=input(us=), A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]; B=[b1;0;0]; I=A\B*us; Ia=I(1);Ib=I(2);Ic=I(3); Is=Ia+Ib;Io=Ib-Ic; display(解决问题2); a11=0;a12=1;a13=-1; a21=-r2;a22=r1+r2+r3+r4;a23=-r4; a31=-r5;a32=-r4;a33=r4+r5+r6+r7; IO=input(IO=); A=[a11,a12,a13;a21,a22,a23;a31,a32,a33]; c1=1;c2=0;c3=0; C=[c1;0;0]; I=A\C*IO; Ia=I(1);Ib=I(2);Ic=I(3); Is=Ia+Ib; us=(r2+r5)*Ia-Ib*r2-Ic*r5; (1)、Is=12，Io=-3, Us=-32，Is=-8 2.如图2所示，已知,,,电容两端电压,求,，和，。并画出相量图。要求进行理论分析，并将MATLAB计算结果和理论计算结果进行比较。 图2 z1=3*j;z2=5;z3=-2j;uc=10*exp(30j*pi/180); z23=z2*z3/(z2+z3);z=z1+z23; Ir=uc/z2,Ic=uc/z3,I=Ic+Ir,ul=I*z1,us=I*z disp(uc Ir Ic I ul us) disp(幅值),disp(abs([uc,Ir,Ic,I,ul,us])) disp(相角),disp(angle([uc,Ir,Ic,I,ul,us])*180/pi) ha=compass([uc,Ir,Ic,I,ul,us]); set(ha,linewidth,3) 3．已知某系统冲激响应为。以-0.5s开始，宽度为1.5s，幅度为1的矩形脉冲激励该系统，求输出