实验四MATLAB程序设计.docVIP

实验四 MATLAB程序设计 1.实验目的 （1）熟练掌握MATLAB的程序流程控制结构掌握M文件的结构熟练掌握函数调用和参数传递2.实验仪器 （1）Matlab6.5应用软件安装版 一套 （3）PC机 一台 3. 实验原理 依据MATLAB的编程语言的原理，进行数字、符号运算，完成曲线绘制。 4. 实验步骤 （1）绘制在不同给定参数时，系统的响应图。 （2）使用plot3、mesh、surf指令绘图。 （3）利用指令plot3绘制绿色实线图。 5. 实验报告内容 （1）分别用for和while循环语句计算的程序，再写出一种避免循环的计算程序。（提示：可考虑利用MATLAB的sum（X,n）函数，实现沿数组X的第n维求和。）保存为M文件. sum(x,1) 列求和；sum(x,2)行求和 （2）将课本例4-8子函数编程及调用演示，三个字图上的圆和多表性绘制在同一坐标系中。并保存为M文件。保存输出结果图。 （3）利用置换指令subs(X,new),例4-18中的脉冲响应在t=[0,18]的曲线。并保存为M文件。保存输出结果图。 （4）试利用feval（）指令计算F(x)+F(x)，其中F可取‘sin’、‘cos’。（提示：先编写以个M函数function y=trif(F,x)实现F(x)+F(x)的计算，在编写调用函数完成F为‘sin’、‘cos’的计算），并保存为M文件。 s=0 for i=0:63 b=pow2(i) s=s+b end i=0 sum=0 while i=63 a=pow2(i) sum=sum+a i=i+1 end a=0:63; a=pow2(a); sum(a,2) function y1=mainfun(a,s) t=(0:a)/a*2*pi y1=subfun(4,s) function y2=subfun(a,s) t=(0:a)/a*2*pi ss=a*exp(i*t) switch s case {base,caller} y2=evalin(s,ss) case self y2=eval(ss) end a=30 t=(0:a)/a*2*pi sss={base,caller,self} for k=1:3 y0=mainfun(8,sss{k}) plot(real(y0),imag(y0),r,LineWidth,4) hold on axis square end