c31.m

来自「matlab例程MATLAB 6.x版的内建数据类型（Built-in data」· M 代码 · 共 24 行

% (c31.m)
%M is the number of coefficients to be computed
T0=pi;N0=256;T=T0/N0;M=10;
t=0:T:T*(N0-1);t=t';
f=exp(-t/2);f(1)=0.604;
Dn=fft(f/N0)
[Dnangle,Dnmag]=cart2pol(real(Dn),imag(Dn));
k=0:length(Dn)-1;k=k';
subplot(211),stem(k,Dnmag)
subplot(212), stem(k,Dnangle)

%Computation of trigonometric Fourier coefficients:

%c31;clg
C0=Dnmag(1); Cn=2*Dnmag(2:M);
Amplitudes=[C0;Cn]
Angles=Dnangle(1:M);
Angles=Angles*(180/pi);
disp('Amplitudes Angles')
[Amplitudes Angles]
k=0:length(Amplitudes)-1; k=k';
subplot(211),stem(k,Amplitudes)
subplot(212), stem(k,Angles)