?? dm09301.m
字號:
% dm09301
% 周期方波信號頻譜分析
function CTFS_SQ
% 繪制并觀察周期方波信號頻譜
% Nf:傅里葉級數展開的項數
% an:各次諧波余弦項系數
display('Please input the value of Nf ');
Nf = input('Nf = ');
an = zeros(Nf+1,1);
an(1) = 0;
for i = 1:Nf
an(i+1) = 2/(i*pi)*sin(i*pi/2); %計算系數an
cn(i+1) = abs(2/(i*pi)*sin(i*pi/2)); %計算幅度譜
end
t = -5:0.01:5;
x = square(pi*(t+0.5)); %用square函數求方波信號
subplot(211);
plot(t,x); %繪制方波信號
axis([-5 5 -1.5 1.5]);
title('周期方波信號','Fontsize',8);
xlabel('t (單位:s)', 'Fontsize',8);
subplot(212);
k = 0:Nf;
stem(k,cn);
hold on;
plot(k,cn);
title('幅度頻譜','Fontsize',8);
xlabel('諧波次數', 'Fontsize',8);
% End
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -