matlab上机实验4,CAD上机实验报告4

CAD作业4-35

一、题目

利用Matlab程序和尖顶余弦脉冲的分解公式：

sin cos ( ) 0 (1 cos ) 1( ) sin cos

(1 cos ) 2(sinn cos ncosn sin ) n( ) n (n2 1)(1 cos )

画出尖顶余弦脉冲分解系数图。

二、matlab程序

t=0:0.001:1;

rad=t*pi;

n=6;

A=zeros(6,1001);

A(1,:)=(sin(rad)-rad.*cos(rad))/pi./(1-cos(rad));

A(2,:)=(rad-sin(rad).*cos(rad))/pi./(1-cos(rad));

for i=3:n

A(i,:)=2*(sin(i*rad).*cos(rad)-i*cos(i*rad).*sin(rad))/i/pi/(i*i-1)./(1-cos(rad)); end;

t=t*180;

plot(t,A(1,:),'r');hold on;

plot(t,A(2,:),'b');hold on;

plot(t,A(3,:),'g');hold on;

plot(t,A(4,:),'m');hold on;

plot(t,A(5,:), 'k');hold on;

grid on;

xlabel('θ°');

ylabel('αn');

title('余弦脉冲的谐波分解系数');

三、程序截图

实验四 第1页