
Добавил:
PRO100TIT
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:курсовая работа / kurs_polinom_N_4_poryadka
.m Tn=27*pi;
deltaT=0.015;
Tv=28*pi;
t=Tn:deltaT:Tv;
x=2*cos(2*t)+2*sin(sqrt(2)*t);
figure(1);
plot(t,x,'k');
hold on;
grid on;
t0=Tn;
dt=pi/4;
for j=1:1:5
f(j)=2*cos(2*(dt*(j-1)+27*pi))+2*sin(sqrt(2)*(dt*(j-1)+27*pi));
plot(27*pi+pi.*(j-1)/4,f(j),'gh');
end
t1=27*pi;
t2=27*pi+pi/4;
t3=27*pi+2*pi/4;
t4=27*pi+3*pi/4;
t5=27*pi+4*pi/4;
p=(f(1)*(t-t2).*(t-t3).*(t-t4).*(t-t5))/((t1-t2).*(t1-t3).*(t1-t4).*(t1-t5))+(f(2)*(t-t1).*(t-t3).*(t-t4).*(t-t5))/((t2-t1).*(t2-t3).*(t2-t4).*(t2-t5))+(f(3)*(t-t1).*(t-t2).*(t-t4).*(t-t5))/((t3-t1).*(t3-t2).*(t3-t4).*(t3-t5))+(f(4)*(t-t1).*(t-t2).*(t-t3).*(t-t5))/((t4-t1).*(t4-t2).*(t4-t3).*(t4-t5))+(f(5)*(t-t1).*(t-t2).*(t-t3).*(t-t4))/((t5-t1).*(t5-t2).*(t5-t3).*(t5-t4));
plot(t,p,'g.');
hold on;
title ('График полинома Лагранжа 4 порядка');
legend ('Сигнал','Полином 4 порядка','Location','Best');
hold on;
ep3=x-p;%эмпирическая погрешность
figure(2);
plot(t, ep3, 'k--');
title ('График эмпирической погрешности');
legend('4 порядок');
deltaT=0.015;
Tv=28*pi;
t=Tn:deltaT:Tv;
x=2*cos(2*t)+2*sin(sqrt(2)*t);
figure(1);
plot(t,x,'k');
hold on;
grid on;
t0=Tn;
dt=pi/4;
for j=1:1:5
f(j)=2*cos(2*(dt*(j-1)+27*pi))+2*sin(sqrt(2)*(dt*(j-1)+27*pi));
plot(27*pi+pi.*(j-1)/4,f(j),'gh');
end
t1=27*pi;
t2=27*pi+pi/4;
t3=27*pi+2*pi/4;
t4=27*pi+3*pi/4;
t5=27*pi+4*pi/4;
p=(f(1)*(t-t2).*(t-t3).*(t-t4).*(t-t5))/((t1-t2).*(t1-t3).*(t1-t4).*(t1-t5))+(f(2)*(t-t1).*(t-t3).*(t-t4).*(t-t5))/((t2-t1).*(t2-t3).*(t2-t4).*(t2-t5))+(f(3)*(t-t1).*(t-t2).*(t-t4).*(t-t5))/((t3-t1).*(t3-t2).*(t3-t4).*(t3-t5))+(f(4)*(t-t1).*(t-t2).*(t-t3).*(t-t5))/((t4-t1).*(t4-t2).*(t4-t3).*(t4-t5))+(f(5)*(t-t1).*(t-t2).*(t-t3).*(t-t4))/((t5-t1).*(t5-t2).*(t5-t3).*(t5-t4));
plot(t,p,'g.');
hold on;
title ('График полинома Лагранжа 4 порядка');
legend ('Сигнал','Полином 4 порядка','Location','Best');
hold on;
ep3=x-p;%эмпирическая погрешность
figure(2);
plot(t, ep3, 'k--');
title ('График эмпирической погрешности');
legend('4 порядок');
Соседние файлы в папке курсовая работа