Ряд Фурье на 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;

Может кто-нибудь объяснить, почему я получаю те разрывы в осциллограмме, когда график обнуляется и как это исправить?

0 ответов

Другие вопросы по тегам