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 日

1 投票

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 日

0 投票

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.

カテゴリ

ヘルプ センター および File Exchange で MATLAB についてさらに検索

製品

リリース

R2021b

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by