山西大同大学matlab课程结课作业
MATLAB程序应用
姓 名：
课程序号： 班 级： 学 号： 2
2013年12月
1.实验内容：已知n!123n，编写一个程序求满足n!10最大的n值以及此时n!的值。
function n
n=2;m=1;
while m
m=m.*n;n=n+1;
end
m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 100的
2.设22(11
31511111)79111315，试根据公式编出计算pi的Mat
lab主程序文件，pi的精度为0.00001。
程序：
k=0;n=1;b=0;a=0;
while abs((pi-a))>0.00001
a=2*sqrt(2)*k;
k=( bcos( *pi/2)+sin(b*pi/2))/n+k;
n=n+2;
b=b+1;
end
a
输出a=3.[**************] ; a-pi= 9.[**************]e-006
1-16231A2453.有两个矩阵A和B如下：1178077，37652131B00132610，
将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
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、W分别存放A、B中间变换结果。
;
4.用matlab主程序文件产生动画：呈现一小圆（半径为
1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。
clea close;clc;r;
axis([-6 6 -6 6],'equal','manual');hold on;
ezplot('x^2+y^2-9');
h=ezplot('x^2+y^2-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
5.matlab画图命令，以1,0
执行语句并画出图形。
提示：利用fill命令绘图。
theta=linspace(0,2*pi,8);
a=1+sin(theta);
b=cos(theta);
plot(a,b,'r');
axis equal;
fill(a,b,'r')
为心，绘制一个正七边形，并填涂成红色。写出
6、（自选题）根据自己专业的信号与系统分析以及个人兴趣，完成一个matlab的综合应用
【题目】
证明 时域卷积定理
【思路和方法】
用矩阵形式表示两个函数
算出两个函数的卷积
扩充两个函数的范围以及卷积的范围
写出卷积的傅里叶变换形式
写出两个函数的傅里叶变换形式
比较卷积的傅里叶变换形式和两个函数傅里叶变换的乘积是否相等
【源程序】
f1= [1,2;3,4];
f2= [5,6;7,8];
f3=conv2(f1, f2); %两个函数的卷积
f3(4,4)=0 %矩阵扩充
Fa=fft2(f3) %卷积的傅里叶变换
f1(4,4) = 0; %矩阵扩充
f2(4,4) = 0;%矩阵扩充
F1= fft2(f1);
F2 = fft2(f2);
Fb= F1 .* F2 %两个函数傅里叶变换的乘积
Fa-Fb
【运行结果】
f3 =
5 16 12 0
22 60 40 0
21 52 32 0
0 0 0 0
Fa =
1.0e+002 *
2.6000 -0.3600 - 1.2800i 0.0400 -0.3600 + 1.2800i -0.7200 - 1.2200i -0.5600 + 0.5400i 0 - 0.0200i 0.6400 - 0.1800i 0.1600 0 - 0.0800i 0 0 + 0.0800i
-0.7200 + 1.2200i 0.6400 + 0.1800i 0 + 0.0200i -0.5600 - 0.5400i
Fb =
1.0e+002 *
2.6000 -0.3600 - 1.2800i 0.0400 -0.3600 + 1.2800i -0.7200 - 1.2200i -0.5600 + 0.5400i 0 - 0.0200i 0.6400 - 0.1800i 0.1600 0 - 0.0800i 0 0 + 0.0800i
-0.7200 + 1.2200i 0.6400 + 0.1800i 0 + 0.0200i -0.5600 - 0.5400i
ans =
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
7.意见和建议
尊敬的老师您好：
matlab的学习已经就要结束了，感觉时间飞快，虽然也经历了好几周，
但我总感觉时间过的太快，这么说的原因很简单，因为这是我上大学以来唯一一门感觉学着有意思而且用处颇多的课程，还有很多课能用matlab，这让我再一次认识到这个软件的强大与简单易懂，以往枯燥无味而且复杂的自控原理也头一次在这个软件中看起来不那么讨厌了，我没有上网查一大堆的有关MATLAB
的资料，简单的几句话就是我学习MATLAB的最真实感想，我喜欢这个软件，我也会花更多的时间去学习这个软件，相信它会给我更大的裨益，虽然实验课上的不理想，使之难以消化，但我想通过我的认真学习，也不算什么，我相信借助matlab它以后会让我们的学习更方便，更简单
山西大同大学matlab课程结课作业
MATLAB程序应用
姓 名：
课程序号： 班 级： 学 号： 2
2013年12月
1.实验内容：已知n!123n，编写一个程序求满足n!10最大的n值以及此时n!的值。
function n
n=2;m=1;
while m
m=m.*n;n=n+1;
end
m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 100的
2.设22(11
31511111)79111315，试根据公式编出计算pi的Mat
lab主程序文件，pi的精度为0.00001。
程序：
k=0;n=1;b=0;a=0;
while abs((pi-a))>0.00001
a=2*sqrt(2)*k;
k=( bcos( *pi/2)+sin(b*pi/2))/n+k;
n=n+2;
b=b+1;
end
a
输出a=3.[**************] ; a-pi= 9.[**************]e-006
1-16231A2453.有两个矩阵A和B如下：1178077，37652131B00132610，
将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
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、W分别存放A、B中间变换结果。
;
4.用matlab主程序文件产生动画：呈现一小圆（半径为
1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。
clea close;clc;r;
axis([-6 6 -6 6],'equal','manual');hold on;
ezplot('x^2+y^2-9');
h=ezplot('x^2+y^2-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
5.matlab画图命令，以1,0
执行语句并画出图形。
提示：利用fill命令绘图。
theta=linspace(0,2*pi,8);
a=1+sin(theta);
b=cos(theta);
plot(a,b,'r');
axis equal;
fill(a,b,'r')
为心，绘制一个正七边形，并填涂成红色。写出
6、（自选题）根据自己专业的信号与系统分析以及个人兴趣，完成一个matlab的综合应用
【题目】
证明 时域卷积定理
【思路和方法】
用矩阵形式表示两个函数
算出两个函数的卷积
扩充两个函数的范围以及卷积的范围
写出卷积的傅里叶变换形式
写出两个函数的傅里叶变换形式
比较卷积的傅里叶变换形式和两个函数傅里叶变换的乘积是否相等
【源程序】
f1= [1,2;3,4];
f2= [5,6;7,8];
f3=conv2(f1, f2); %两个函数的卷积
f3(4,4)=0 %矩阵扩充
Fa=fft2(f3) %卷积的傅里叶变换
f1(4,4) = 0; %矩阵扩充
f2(4,4) = 0;%矩阵扩充
F1= fft2(f1);
F2 = fft2(f2);
Fb= F1 .* F2 %两个函数傅里叶变换的乘积
Fa-Fb
【运行结果】
f3 =
5 16 12 0
22 60 40 0
21 52 32 0
0 0 0 0
Fa =
1.0e+002 *
2.6000 -0.3600 - 1.2800i 0.0400 -0.3600 + 1.2800i -0.7200 - 1.2200i -0.5600 + 0.5400i 0 - 0.0200i 0.6400 - 0.1800i 0.1600 0 - 0.0800i 0 0 + 0.0800i
-0.7200 + 1.2200i 0.6400 + 0.1800i 0 + 0.0200i -0.5600 - 0.5400i
Fb =
1.0e+002 *
2.6000 -0.3600 - 1.2800i 0.0400 -0.3600 + 1.2800i -0.7200 - 1.2200i -0.5600 + 0.5400i 0 - 0.0200i 0.6400 - 0.1800i 0.1600 0 - 0.0800i 0 0 + 0.0800i
-0.7200 + 1.2200i 0.6400 + 0.1800i 0 + 0.0200i -0.5600 - 0.5400i
ans =
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
7.意见和建议
尊敬的老师您好：
matlab的学习已经就要结束了，感觉时间飞快，虽然也经历了好几周，
但我总感觉时间过的太快，这么说的原因很简单，因为这是我上大学以来唯一一门感觉学着有意思而且用处颇多的课程，还有很多课能用matlab，这让我再一次认识到这个软件的强大与简单易懂，以往枯燥无味而且复杂的自控原理也头一次在这个软件中看起来不那么讨厌了，我没有上网查一大堆的有关MATLAB
的资料，简单的几句话就是我学习MATLAB的最真实感想，我喜欢这个软件，我也会花更多的时间去学习这个软件，相信它会给我更大的裨益，虽然实验课上的不理想，使之难以消化，但我想通过我的认真学习，也不算什么，我相信借助matlab它以后会让我们的学习更方便，更简单