MATLAB第一次实验报告

1、matlab第一次实验报告 matlab第一次试验报告 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x，y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)2); y(n+1) = b*x(n)+a*(y(n)-x(n)2) 程序：function y1,y2=draw(a,b,c,x1,y1) x=; y=; x(1)=x1; y(1)=y1; n=c; for i=1:n x(i+1)=a*x(i)-b*(y(i)-x(i)2); y(i+1)=b*x(i)+a*(y(i)-x(i)2); end y1=x; y2=y; plot(y1,y2)

2、输入：draw（5,1,1,1,2） 2. 编程实现奥运5环图，允许用户输入环的直径。 程序： clc clear r=input('input a mnmber:') color1='b','k','r' color2='y','g' t=linspace(0,2*pi); a1=r+0.2*r; a2=3*r+0.4*r; a3=5*r+0.6*r; a=a1,a2,a3; b1=2*r+0.3*r; b2=4*r+0.5*r; b=b1,b2; for i=1:3 x=r*cos(t)+a(i

3、); y=r*sin(t); plot(x,y,color1(i),'linewidth',5) axis(0,7*r,-2*r,1.5*r); hold on axis equal end for n=1:2 x=r*cos(t)+b(n); y=r*sin(t)-r; plot(x,y,color2(n),'linewidth',5) hold on end axis off 输入：6（圆的半径） 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort函数。 程序： function y=qipaofa(s) s n=length(s);

4、for i=1:n for j=1:n-i if s(j)s(j+1) t=s(j ); s(j)=s(j+1); s(j+1)=t; end end end s 输入：qipaofa(1 4 2 6 8 3) 输出：s = 1 4 2 6 8 3 s = 1 2 3 4 6 8 总结与改进： 本次试验共完成三个程序，试验中的题都要求输入不同的数据所以要求娴熟的使用funtion函数，是程序可以输入不同的值而得到不同的输出结果。还有plot函数画图的使用以及对图像输出颜色，宽度等的调整。还有递归法，循环的使用。 此次试验的改进主要是对于五环的程序的改进，五环的图形是环套环但第一个程序未得到这个

5、效果所以对这个程序进行了改进，让五环的每个环分段输出具有环套环的效果。 改进程序： clc clear r=input('input r:'); lw=8; alp=linspace(0.75*pi,1.75*pi); x=r*cos(alp)-1.2*r; y=r*sin(alp)+r; plot(x,y,'b','linewidth',lw)%1 hold on axis equal axis(-2.8*r,5.4*r,-1.6*r,2.6*r); alp=linspace(pi/2,2.2*pi); x=r*cos(alp)+0.01*r;

6、 y=r*sin(alp); plot(x,y,'y','linewidth',lw)%1 alp=linspace(-0.3*pi,0.8*pi); x=r*cos(alp)-1.2*r; y=r*sin(alp)+r; plot(x,y,'b','linewidth',lw)%2 alp=linspace(0.75*pi,1.75*pi); x=r*cos(alp)+1.25*r; y=r*sin(alp)+r; plot(x,y,'k','linewidth',lw)%1 alp=linspa

7、ce(0.15*pi,0.5*pi); x=r*cos(alp)+0.01*r; y=r*sin(alp); plot(x,y,'y','linewidth',lw)%2 alp=linspace(0.5*pi,2.2*pi); x=r*cos(alp)+2.49*r; y=r*sin(alp); plot(x,y,'g','linewidth',lw)%1 alp=linspace(-0.3*pi,0.8*pi); x=r*cos(alp)+1.25*r; y=r*sin(alp)+r; plot(x,y,'k','linewidth',lw)%2 alp=linspace(0.75,1.75*pi); x=r*cos(alp)+3.7*r; y=r*sin(alp)+r; plot(x,y,'r','linewidth',lw)%1 alp=linspace(0.15*pi,0.5*pi); x=r*cos(alp)+2.49*r; y=r*sin(alp); plot(x,y,'g','linewidth',lw)%2 alp=linspace(-0.3*pi,0.8*