matlab作业

第一题

求解线性方程组AX=B 其中20

900

83

6

215.8733322

3107

4

45.06.307925

.11---=A ，165204

3-=B

程序：a=[1 1.5 2 9 7;0 3.6 0.5 -4 4;7 10 -3 22 33;3 7 8.5 21 6;3 8 0

90 -20];

b=[3 -4 20 5 16]'; A=det(a); for i=1:5 D=a;

D(:,i)=b;

x(i)=det(D)/A; end x

运行结果：

第二题

x (x<1) 计算函数值 2x-1 (1=10) 程序：display('第二题') x=input('输入x=') if x<1 x

elseif x<10 2*x-1 else

3*x-11

End

运行结果：

第三题

求一元二次方程02

=++c bx ax

程序：a=input('a=');

b=input('b='); c=input('c='); d=b*b-4*a*c;

x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];

disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);

x = 3.5653 -0.9255 -0.2695 0.1435 0.0101

第二题 输入x=2 x = 2 ans = 3

运行结果：

第四题

输入三角形的三条边，求面积。 程序：

disp('请输入三角形的三条边') a=input('a='); b=input('b='); c=input('c=');

if a+b>c & b+c>a & c+a>b d=(a+b+c)/2;

s=sqrt(d*(d-a)*(d-b)*(d-c)); disp(['三角形的面积是：',num2str(s)]) else a+b<=c|b+c<=a|c+a<=b

disp('不能构成三角形') end

运行结果：

第五题

输入20个数，求其中最大数和最小数。要求分别用循环结构和调用MATLAB 的max 函数、min

函数来实现。 程序：

a=rand(1,20)

m=a(1);%最大值 n=a(1);%最小值 for i=2:20 if a(i)>m m=a(i); end

if a(i)

n

运行结果：

a=2 b=3 c=4

x1=-0.75+1.199i,x2=-0.75-1.199i

请输入三角形的三条边

a=8 b=7 c=8 三角形的面积是：25.1781 a =

Columns 1 through 19

0.8147 0.9058 0.1270 0.9134 0.6324

0.0975 0.2785 0.5469 0.9575 0.9649

0.1576 0.9706 0.9572 0.4854 0.8003

0.1419 0.4218 0.9157 0.7922

Column 20 0.9595 m =0.9706 n =0.0975

第六题

随机产生一个一行20列的数组，并将数列的元素从大到小排列后输出。 程序：x=rand(1,20) for i=1:20 for j=1:20

if x(i)>x(j) min=x(j); x(j)=x(i); x(i)=min; end end end

disp(x)

运行结果：

第七题

求Fibonacci 数列(1)大于4000的最小项。(2)5000之内的项数。 程序：f(1)=1; f(2)=1; for i=3:30

f(i)=f(i-1)+f(i-2); end f x=f

y=length(x); for i=1:y

if x(i)>4000&x(i)<5000 b=x(i);

disp('min='),disp(b(1)) disp('i='),disp(i) end end

x =

Columns 1 through 13 0.3111 0.9234 0.4302 0.1848 0.9049 0.9797 0.4389 0.1111 0.2581 0.4087 0.5949 0.2622 0.6028 Columns 14 through 20 0.7112 0.2217 0.1174 0.2967 0.3188 0.4242 0.5079 Columns 1 through 13 0.9797 0.9234 0.9049 0.7112 0.6028 0.5949 0.5079 0.4389 0.4302 0.4242 0.4087 0.3188 0.3111 Columns 14 through 20 0.2967 0.2622 0.2581 0.2217 0.1848 0.1174 0.1111

运行结果：

第八题

绘制下列曲线:

（1）

2

1100x +；（2）2

221x e -

π

；（3）???==32

5t

y t x

程序（1）：x=linspace(-10,10,1000); y=100./(1+x.^2); plot(x,y) 运行结果：

程序（2）：x=linspace(-10,10,1000);

y=(1/(2*pi))*(exp(-0.5.*x.^2));

f =

Columns 1 through 11

1 1

2

3 5 8 13 21 3

4 5

5 89

Columns 12 through 22

144 233 377 610 987 1597

2584 4181 6765 10946 17711

Columns 23 through 30

28657 46368 75025 121393

196418 317811

514229 832040

x =

Columns 1 through 11

1 1

2

3 5 8 13

21 34 55 89

Columns 12 through 22

144 233 377 610 987 1597

2584 4181 6765 10946 17711

Columns 23 through 30

28657 46368 75025 121393 196418 317811

514229 832040

min=

4181

i=

19

plot(x,y) 运行结果：

程序（3）：t=linspace(-10,10,1000); x=t.^2; y=5.*t.^3; plot(x,y) 运行结果：

第九题

在同一图形窗口采用子图形式分别绘制正方形、圆、三角形和六边形。 程序： 运行结果：

t1=0:pi/10:2*pi;

t2=0:pi/2:2*pi; t3=0:2*pi/3:2*pi; t4=0:pi/3:2*pi;

subplot(2,2,1),plot(exp(i*t1)) axis equal

subplot(2,2,2),plot(exp(i*t2)) axis equal

subplot(2,2,3),plot(exp(i*t3)) axis equal

subplot(2,2,4),plot(exp(i*t4)) axis equal

绘制下列极坐标图。（1）4cos 5+=θρ；（2）

θ

12

；（3）23

θπ

ρ=

程序：t=0:pi/60:2*pi;

polar(t,5*cos(t)+4) ,text(5,5,'5cos(\theta)+4') hold on

polar(t,12./(sqrt(t)));text(-5,-5,'12/sqrt(\theta)') polar(t,pi/3.*(t.^2));text(-5,3,'\pi\theta^2/3') 运行结果：

第十一题

某工厂2005年度各季度生产值（单位：万元）分别为：450.6、395.9、410.2、450.9，试绘制折线图和饼图，并说明图形的实际意义。 程序：x=[450.6 395.9 410.2 450.9] subplot(1,2,1);plot(x)

xlabel('季度');ylabel('生产值')

subplot(1,2,2);pie([450.6 395.9 410.2 450.9]);legend('第一季度','第二季度','第三季度','第四季度') 运行结果：

根据1252

2

22=-+

a y a x 绘制平面曲线，并分析参数对其形状的影响。 程序：a=[1 2 3 4]; for i=1:4

subplot(2,2,i)

ezplot(['(x.^2)./(',num2str(a(i)),'.^2)+(y.^2)/(25-(',num2str(a(i )),'.^2))=1']); axis equal end 运行结果： 第十三题

绘制下列三维图形。（1）?????===t z t y t x sin cos ；（2）??

?

??=+=+=u z v u y v u x sin sin )cos 1(cos )cos 1(；(3) 半径为10的球面

程序（1）：t=0:pi/30:8*pi;x=cos(t);y=sin(t);z=t;

plot3(x,y,z,'color','k','markersize',20)

运行结果：

程序（2）：u=0:pi/30:2*pi v=0:pi/30:2*pi

x=(1+sin(u)).*cos(v)

y=(1+cos(v)).*sin(u) z=sin(u)

plot3(x,y,z) 运行结果：

程序（3）：[u,v]=meshgrid(0:pi/100:4*pi);运行结果： r=10

x=r.*cos(u).*cos(v); y=r.*cos(u).*sin(v); z=r.*sin(u); surf(x,y,z) axis equal grid on

上机考试程序：

1、绘制李莎如图形：

程序：t=linspace(0,2*pi,100);

e={['=0'],['=\pi/4'],['=\pi/2'],['=3\pi/4'],['=\pi']};%e={'0','\p i/4','3\pi/4','\pi'}; f=['\phi']; for i=1:5 x=cos(t);

y=cos(t+(i-1)*pi/4);

subplot(2,3,i),plot(x,y) xlabel('x'),ylabel('y') ff=e(i); %title(e(i)) h=[f,ff]; title(h)

subplot(2,3,6),plot(x,y) hold on title('总图')

xlabel('x'),ylabel('y') End 运行结果：

2、程序

t=linspace(0,3*pi,100); xx=t-sin(t); yy=1-cos(t); x=sin(t); y=cos(t); for i=1:100

plot(xx(1:i),yy(1:i)) hold on

plot(xx(i),yy(i),'r.','markersize',20)

plot(t(i)+x,y+1,t,0,'b');plot([0 13],[0 0],'r') axis([0,10,-1,1]); axis equal ; hold off

m(i)=getframe; end ;

movie(m); 运行结果：

3、程序：t=0:0.1:3*pi; v=50;g=9.8;

theta=[pi/6,pi/4,pi/3 ]; e={'=30','=45','=60'}; f={'\theta'}; c=['g','r','b']; for i=1:3

x=v*cos(theta(i)).*t;

y=v*sin(theta(i)).*t-g/2.*t.*t; plot(x,y,c(i))

axis([0,300,0,150]) hold on ff=e(i) h=[f,ff]

text(v*cos(theta(i))*(2*max(y)/g)^0.5,max(y),h) end

title('不同抛射角的斜抛运动轨道') xlabel('水平距离X') ylabel('

垂直距离Y') 运行结果：

考试程序：

问题：如下图所示，在一单位长度的线段上，对其三等分，将中间段直线换成一个去掉底边的三角形，再在每条直线上重复以上操作，如此进行下去直到无穷，就得到分形曲线科赫曲线。

思路及说明：在将科赫雪花曲线程序弄懂以后，编写此程序。本程序与科赫雪花曲线有相似之处，所以在方法上是引用科赫雪花曲线程序的方法。

程序：new=[0,i];%建立一个数组

subplot(3,3,1);

plot(new);%画出第一幅图形

old=new;%准备画第二幅图，故将第一幅图的点转为旧的点

diff=(old(2)-old(1))/3;%计算第一幅图中1/3处点的坐标

new(1)=old(1);%画出第二幅图中第一个点

new(2)=old(1)+diff;%画出第二幅图中第一、二个点

new(3)=new(2)+diff*((1+sqrt(3)*i)/2);%画出第二幅图中第一、二、三个点

new(4)=old(1)+2*diff;%画出第二幅图中第一、二、三、四个点

new(5)=old(1)+3*diff;%画出第二幅图中第一、二、三、四、五个点

subplot(3,3,2);

plot(new);%画出第二幅图

axis equal

n=[new(1),new(2),new(3),new(4),new(5)];%建立一个新数组

for k=3:9;%第k幅图

o=n;

a=length(o)-1;%a为线段的数目

for j=0:a-1;%对每条线段进行计算

d=(o(j+2)-o(j+1))/3;%线段1/3的长度

n(4*j+1)=o(j+1);%起始点的坐标

n(4*j+2)=o(j+1)+d;%线段的1/3处

n(4*j+3)=n(4*j+2)+d*((1+sqrt(3)*i)/2);%移动后的点的坐标 n(4*j+4)=o(j+1)+2*d;

end

n(4*a+1)=[i];%终点的坐标

subplot(3,3,k);

plot(n);

axis equal

End

运行结果：