Ряд Фурье на MATLAB
Я искал повсюду, а также пробовал несколько вещей, чтобы найти решение проблемы, с которой я столкнулся. Я не могу построить ряд Фурье, который должен дать мне трапециевидную форму волны. Но это прерывается между ними. Вот код:
Br = 0.245;
Mu0 = 4*pi*10^-7;
r2 = 13.25e-3;
r4 = 16.75e-3;
W0 = 1;
r = 17.75e-3;
R = r2/r;
theta = 60;
psi = -2*pi:0.01:2*pi;
phi = 0;
syms x k;
clear M;
Mk = (4/pi)*(Br/Mu0)*(int(cos(k*x),0,pi/10)+int((5-40*x/pi).*cos(k*x),pi/10,3*pi/20)-int(cos(k*x),3*pi/20,7*pi/20)+int((40*x/pi-15).*cos(k*x),7*pi/20,2*pi/5)+int(cos(k*x),2*pi/5,pi/2));
Mr = symsum(Mk*R/Mu0*cos(k*psi),k,1,100);
figure;
plot(psi,Mr);
grid on;
Может кто-нибудь объяснить, почему я получаю те разрывы в осциллограмме, когда график обнуляется и как это исправить?