matlab结课论文.doc

: .
页眉内容
山西大同大学matlab课程结课作业
MATLAB程序应用
姓 名：
课程序号： 2
班 级：
学 号：
2013年12月
页眉内容
1. 实验内容：已知n! 1 2 3 L n，编写一个程序求满足n! 10100的 最大的n值以及此时n!的值。
fun cti on n
n=2;m=1;
while m<=10A100
m=m.* n;n=n+1;
end
m=m/(n-1) ； n=n-2;
m n
m = +098 n =69
2^(1 1111丄丄丄)
2. 设 3 5 7 9 11 13 15 ，试根据公式编出计算pi的Mat
lab主程序文件，。
程序：
k=0 ;n=1;b=0;a=0;
while abs((pi-a))>
a=2*sqrt(2)*k;
k=( bcos( *pi/2)+si n(b*pi/2))/n+k;
n=n+2;
b=b+1;
end
a
输出 a= ; a-pi= -006
1
-1
6
8
3
7
6
5
2 A
3
1
0
B
2
1
3
1
2
4
5
7
0
0
1
3
： 1
1
7
7
2
6
1
0
将A中所有等于-1的元素改为-2，将B中所有小于0的元素改为1，然后将B中等于 0的元素的值改为A的相应位置元素的值。请用Matlab函数文件实现上述运算。
clear;
clc;
A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0];
C=A;A(A==-1)=-2;U=A;
D=B;B(B<0)=1;V=B;
A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D;
A,B,W,U,V %用函数文件实现矩阵中元素的变换。
%A B为输入变量。
%U V、V分别存放A B中间变换结果。
4. 用matlab主程序文件产生动画：呈现一小圆(半径为 1)在一大圆(半径为3
)的圆周外部滚动的动画，要求连续滚动 20周。
clea close;clc;r;
axis([-6 6 -6 6],'equal','manual');hold on;
ezplot('xA2+yA2-9'); h=ezplot('xA2+yA2-1');
x=get(h,'xdata');
y=get(h,'ydata');
for t=1:7200
set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t));
drawnow;
end
1,0
5. matlab 画图命令，以 为心，绘制一个正七边形，并填涂成红色。写出 执行语句并画出图形。
提示：利用 fill 命令绘图。
thet