matlab结课论文

文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. 山西大同大学matlab课程结课作业

MATLAB程序应用

姓名：

课程序号： 2

班级：

学号：

2013年12月

1.实验内容：已知!123n n =⨯⨯⨯⋅⨯L ，编写一个程序求满足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.2083633 ; a-pi= 9.9577991e-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它以后会让我们的学习更方便，更简单