
Добавил:
PRO100TIT
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:курсовая работа / kurs_polinom_N_1_poryadka
.m Taun=27*pi;
deltaTau=0.015;
Tauv=28*pi;
Tau=Taun:deltaTau:Tauv;
x=2*cos(2*Tau)+2*sin(sqrt(2)*Tau);%исходный сигнал от Tau
figure(1);
plot(Tau,x,'k');
hold on;
grid on;
t=-1:deltaTau:1;
xt=2*cos(pi*t+55*pi)+2*sin(sqrt(2)*(pi*t/2+27.5*pi));%сигнал от t
figure(2);
plot(t,xt,'k');
hold on;
grid on;
t0=-1;%узловые точки
t1=1;
f0=3.0916;
f1=0.094;
p=(f0*(t-t1)/(t0-t1))+(f1*(t-t0)/(t1-t0));
plot(t,p,'g.');
hold on;
title ('График полинома Лагранжа 1 порядка');
legend ('Сигнал','Полином 1 порядка','Location','Best');
ep1=xt-p;%эмпирическая погрешность
dx2=2*pi^2*cos(pi*t)-pi^2*sin(sqrt(2)*(pi*t/2+55/2*pi));%оценка сверху максимальной погрешности
m2=max(abs(dx2));
q2=(t-t0).*(t-t1);
es1=m2*max(abs(q2))/2;
figure(3);
plot([-1 1], [es1 es1], 'k--');
hold on;
plot(t, ep1, 'k--');
title ('График эмпирической погрешности и график оценки сверху максимальной погрешности');
legend('1 порядок');
deltaTau=0.015;
Tauv=28*pi;
Tau=Taun:deltaTau:Tauv;
x=2*cos(2*Tau)+2*sin(sqrt(2)*Tau);%исходный сигнал от Tau
figure(1);
plot(Tau,x,'k');
hold on;
grid on;
t=-1:deltaTau:1;
xt=2*cos(pi*t+55*pi)+2*sin(sqrt(2)*(pi*t/2+27.5*pi));%сигнал от t
figure(2);
plot(t,xt,'k');
hold on;
grid on;
t0=-1;%узловые точки
t1=1;
f0=3.0916;
f1=0.094;
p=(f0*(t-t1)/(t0-t1))+(f1*(t-t0)/(t1-t0));
plot(t,p,'g.');
hold on;
title ('График полинома Лагранжа 1 порядка');
legend ('Сигнал','Полином 1 порядка','Location','Best');
ep1=xt-p;%эмпирическая погрешность
dx2=2*pi^2*cos(pi*t)-pi^2*sin(sqrt(2)*(pi*t/2+55/2*pi));%оценка сверху максимальной погрешности
m2=max(abs(dx2));
q2=(t-t0).*(t-t1);
es1=m2*max(abs(q2))/2;
figure(3);
plot([-1 1], [es1 es1], 'k--');
hold on;
plot(t, ep1, 'k--');
title ('График эмпирической погрешности и график оценки сверху максимальной погрешности');
legend('1 порядок');
Соседние файлы в папке курсовая работа