matlab结课论文2_可编辑.docxVIP

山西大同大学matlab课程结课作业

MATLAB程序应用

姓名：

课程序号：

2

班级：

学号：

2021年12月

1.实验内容：，编写一个程序求满足的最大的值以及此时的值。

functionn

n=2;m=1;

whilem=10^100

m=m.*n;n=n+1;

end

m=m/(n-1)；n=n-2;

mn

m=1.7112e+098n=69

2.设，试根据公式编出计算pi的Matlab主程序文件，pi的精度为0.00001。

程序：

k=0;n=1;b=0;a=0;

whileabs((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.1416;a-pi=9.918493839577991e-006

3.有两个矩阵A和B如下：，，

将A中所有等于-1的元素改为-2，将B中所有小于0的元素改为1，然后将B中等于0的元素的值改为A的相应位置元素的值。请用Matlab函数文件实现上述运算。

clear;

clc;

A=[1-168;23-10;-2457;1-177];

B=[-3-76-5;-213-1;0013;26-10];

C=A;A(A==-1)=-2;U=A;

D=B;B(B0)=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周。

cleaclose;clc;r;

axis([-66-66],equal,manual);holdon;

ezplot(x^2+y^2-9);

h=ezplot(x^2+y^2-1);

x=get(h,xdata);

y=get(h,ydata);

fort=1:7200

set(h,xdata,x+4*cosd(t),ydata,y+4*sind(t));

drawnow;

end

5.matlab画图命令，以为心，绘制一个正七边形，并填涂成红色。写出执行语句并画出图形。

提示：利用fill命令绘图。

theta=linspace(0,2*pi,8);

a=1+sin(theta);

b=cos(theta);

plot(a,b,r);

axisequal;

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=

516120

2260400

2152320

0000

Fa=

1.0e+002*

2.6000-0.3600-1.2800i0.0400-0.3600+1.2800i

-0.7200-1.2200i-0.5600+0.5400i0-0.0200i0.6400-0.1800i

0.16000-0.0800i00+0.0800i

-0.7200+1.2200i0.6400+0.1800i0+0.0200i-0.5600-0.5400i

Fb=

1.0e+002*

2.6000-