Where is the problem with my code. it never ends !!

Question

geometry geometry 2018-5-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/401426-where-is-the-problem-with-my-code-it-never-ends

回答： Geoff Hayes 2018-5-18

在 MATLAB Online 中打开

What is the problem with the following code?

It calculates the exponential coefficients fo fourier series of periodic function f(t).

when I run it the compilation never ends !! and I must close matlab and run it again.

function [an]=fc(f,T)
syms t n
n=-20:1:20;
an=(1/T)*int(f*exp(-1i*n*2*pi*t/T),t,0,T);
end
syms t  n
 f=rectangularPulse(t/4)*(abs(t)-1);
%f=rectangularPulse(t);
% f=sin(t);
 Y=fc(f,20);
subplot(211)
stem(-20:1:20,abs(Y));
ylabel('|an|')
grid;
subplot(212)
stem(-20:1:20,angle(Y));
ylabel('<an')
grid;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Geoff Hayes 2018-5-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/401426-where-is-the-problem-with-my-code-it-never-ends#answer_320870

在 MATLAB Online 中打开

geometry - your function appears to be recursive

function [an]=fc(f,T)
syms t n
n=-20:1:20;
an=(1/T)*int(f*exp(-1i*n*2*pi*t/T),t,0,T);
end
syms t  n
f=rectangularPulse(t/4)*(abs(t)-1);
%f=rectangularPulse(t);
% f=sin(t);
Y=fc(f,20);    % <------------------