matlab结课论文

山西大同大学matlab课程结课作业MATLAB程序应用
姓名：
课程序号： 2
班级：
学号：
2013年12月
1.实验内容：已知!123n n =⨯⨯⨯⋅⨯ ，编写一个程序求满足100!10n ≤的
最大的n 值以及此时!n 的值。

function n
n=2;m=1;
while m<=10^100
m=m.*n;n=n+1;
end
m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131
1(22 +--++--+=π，试根据公式编出计算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.141602572083633 ; a-pi= 9.918493839577991e-006
3.有两个矩阵A 和B 如下：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=771175420132861-1A ，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=0162310013125673B ，
将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、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、（自选题）根据自己专业的信号与系统分析以及个人兴趣，完成一个matla b的综合应用
【题目】
证明时域卷积定理
【思路和方法】
用矩阵形式表示两个函数
算出两个函数的卷积
扩充两个函数的范围以及卷积的范围
写出卷积的傅里叶变换形式
写出两个函数的傅里叶变换形式
比较卷积的傅里叶变换形式和两个函数傅里叶变换的乘积是否相等
【源程序】
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 1
6 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它以后会让我们的学习更方便，更简单。