实验二、傅里叶级数与傅里叶变换

  • 格式:doc
  • 大小:149.50 KB
  • 文档页数:6

实用标准文案

精彩文档 实验二、傅里叶级数与傅里叶变换

上机作业:

1、 在(0,T)周期上,锯齿波可表为f(x)=x/3,锯齿波可展开为傅里叶级数:

112sin63kTTkxkT

试画出锯齿波的幅度频谱图。

解:

%Fig2d2.m

T=1;x=0.2;H=1;

t=-0.5*T:0.01:0.5*T;

k=0:10;

wk=2*k*pi/T;

Ak=abs(sin(wk*x)./(2*wk/3));

Ak(1)=T/6-pi*x/3;

plot(k,Ak,'b--')

hold on

stem(k,Ak,'o')

xlabel('k')

ylabel('A_k')

title('幅频响应曲线')

set(gca,'xtick',[0:10])

实用标准文案

精彩文档

2、求函数sin2**15*sin2**40*xtt的傅里叶变换,画出其幅频响应及相频响应曲线图。

解:

%Fig2d3.m

x=0:1/100:10-1/100;

y=sin(2*pi*15*x)+sin(2*pi*40*x);

Y=fft(y); %求傅里叶变换

n=fix(length(Y)/2);

p=unwrap(angle(Y));

freq=[0:n-1]./length(Y); 实用标准文案

精彩文档 f=(0:length(Y)-1)'*100/length(Y);

subplot(121)

plot(freq,abs(Y(1:n)),'b-')

xlabel('f')

ylabel('|F(2\pif)|')

title('幅频响应曲线')

subplot(122)

plot(f,p*180/pi)

xlabel('p')

ylabel('|F(2\pip)|')

title('相频响应曲线')

3、已知平面波展开为柱面波的公式是:

cos012cosikmmmeJkiJkm

取上式的实部,令cosx,得 实用标准文案

精彩文档 2021cos2cos2nnnkxJkiJkn

试画出上式左边和右边(n=10和n=20)的柱面波图,并加以比较。

解:

%Fig2d12.m

close all

clear all

[X,Z]=meshgrid(0.05:0.1:10);

g=cos(X);

subplot(221)

contour(g)

meshc(X,Z,g)

xlabel('z')

ylabel('x')

title('向Z方向传播的平面波')

[Q,R]=cart2pol(X,Z);

sqrtR=sqrt(pi/2./R);

Bes0=sqrtR.*besselj(0,R);

qiu=Bes0

for k=2:2:10

Bes=sqrtR.*besselj(k,R); %求k阶贝塞尔函数

qiu=qiu+2*i^k*Bes.*cos(2*k*Q); 实用标准文案

精彩文档 end

subplot(223)

meshc(X,Z,qiu)

xlabel('z')

ylabel('x')

title('n=10时的柱面波图')

for k=12:2:20

Bes=sqrtR.*besselj(k,R); %求k阶贝塞尔函数

qiu=qiu+2*i^k*Bes.*cos(2*k*Q);

end

subplot(224)

meshc(X,Z,qiu)

xlabel('z')

ylabel('x')

title('n=20时的柱面波图') 实用标准文案

精彩文档