Function for fourier series

14 ビュー (過去 30 日間)
Rayan  Maye
Rayan Maye 2022 年 2 月 24 日
コメント済み: Star Strider 2022 年 2 月 24 日
hello
I have written a function for fourier series but when i try to plot it ,only the consine one shows up did I forget something ? Danke schön
function ft = signal(S,t,A,M,T,N)
for n=1:N
if S=="square"
ft=M+sum(((4*A)./(n.*pi)).*sin((n.*pi)/2).*cos((2*n.*pi.*t)./T));
elseif S=="triangle"
ft=M+sum(((8*A)./(n.^2*pi^2)).*cos((2*n*pi.*t)./T));
elseif S=="saw"
ft=M+sum(-1^(n+1).*(2*A./n.^pi).*sin((2*n.*pi*t)./T));
elseif S=="cosine"
ft=M+A*cos((2*pi.*t)./T);
end
end
%%plot part
M=0;
A=1;
T=10;
t=0:0.1:60;
N=10;
grid on
ft=signal("square",t,A,M,T,N);
plot(t,ft,'r')
hold on
ft=signal("triangle",t,A,M,T,N);
plot(t,ft,'b')
hold on
ft=signal("saw",t,A,M,T,N);
plot(t,ft,'g')
hold on
ft=signal("cosine",t,A,M,T,N);
plot(t,ft,'k')
%%

採用された回答

Star Strider
Star Strider 2022 年 2 月 24 日
編集済み: Star Strider 2022 年 2 月 24 日
Try this —
%%plot part
M=0;
A=1;
T=10;
t=0:0.1:60;
N=10;
figure
subplot(4,1,1)
ft=signal("square",t,A,M,T,N);
plot(t,ft,'r')
subplot(4,1,2)
ft=signal("triangle",t,A,M,T,N);
plot(t,ft,'b')
subplot(4,1,3)
ft=signal("saw",t,A,M,T,N);
plot(t,ft,'g')
subplot(4,1,4)
ft=signal("cosine",t,A,M,T,N);
plot(t,ft,'k')
%%
function ft = signal(S,t,A,M,T,N)
n=(1:N).';
if S=="square"
ft=M+sum(((4*A)./(n.*pi)).*sin((n.*pi)/2).*cos((2*n.*pi.*t)./T));
elseif S=="triangle"
ft=M+sum(((8*A)./(n.^2*pi^2)).*cos((2*n*pi.*t)./T));
elseif S=="saw"
ft=M+sum(-1.^(n+1).*(2*A./n.^pi).*sin((2*n.*pi*t)./T));
elseif S=="cosine"
ft=M+A*cos((2*pi.*t)./T);
end
end
The loop was not executing correctly. I eliminated the for loop and changed the code to use ‘n’ as a column vector to take advantage of MATLAB’s vectorisation abilities. (I plotted them as subplots because I wanted to see the individual results more clearly.) Otherwise (except for a minor change in the ‘saw’ calculation) the code is unchanged.
EDIT — (24 Feb 2022 at 15:06)
Changed original answer text to one that actually addresses the problem.
.
  2 件のコメント
Rayan  Maye
Rayan Maye 2022 年 2 月 24 日
How dumb of me of not vectorizing n,thanks for your help mate
Star Strider
Star Strider 2022 年 2 月 24 日
As always, my pleasure!

サインインしてコメントする。

その他の回答 (1 件)

Rayan  Maye
Rayan Maye 2022 年 2 月 24 日
Tried it ,did'nt work but thanks.I'm suspecting the problem is maybe a mistake that I did in the transcription of mathematical formulas.

カテゴリ

Help Center および File ExchangeAnnotations についてさらに検索

タグ

製品


リリース

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by