Matlab code for square wave generation from fourier series (symsum fails for upperlimit greater than 1000 in my case)

32 次查看（过去 30 天）
Rishav Raj
Rishav Raj 2020-1-15
评论： Rishav Raj 2020-3-7
function [a0,An,Bn]=mysquare(T,N)
syms t n
w0=2*pi/T;
a0=1/T*(int(1,t,0,T/2));
an=2/T*(int(1*cos(n*w0*t),t,0,T/2));
bn=2/T*(int(1*sin(w0*n*t),t,0,T/2));
An=symsum(an.*cos(n*w0*t),n,1,N);
Bn=symsum(bn.*sin(n*w0*t),n,1,N);
%main script
T=4;
[a0,an,bn]=mysquare(T,10000);
ft=a0+an+bn;
fplot(ft);
This works fine for N<=1000 but for larger values it takes alot of time and no result.

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

采纳的回答

Charan Jadigam
Charan Jadigam 2020-3-4
Hi,
The function 'symsum' is trying to find a closed form expression as an answer and it is unable to do so.
Instead, as a workaround, use ‘sum’ function to solve the equation. The usage in your example would be,
An=sum(subs(an.*cos(n*w0*t),n,1:10000));
Bn=sum(subs(bn.*sin(n*w0*t),n,1:10000));

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Discrete Fourier and Cosine Transforms 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by