
Добавил:
PRO100TIT
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:курсовая работа / kurs_BestBribl
.m clear();
%Полином наилучшего приблежения
Taun=27*pi; %Нижняя граница ИО
Tauv=28*pi; %Верхняя граница ИО
dt=0.0015;
%Масштабирование
tau = Taun:dt:Tauv;
T0 = 2 * (tau - Taun) / (Tauv - Taun) - 1;
T0kv=(2 * (tau - Taun) / (Tauv - Taun) - 1).^2;
X = 2 * cos (2 * ((T0+ 1) * (pi) / 2 + Taun)) + 2 * sin(sqrt(2) * ((T0 + 1) * (pi) / 2 + Taun));
XT0=X.*(2 * (tau - Taun) / (Tauv - Taun) - 1);
P0 = (max(X)-abs(min(X)))/2;
coef=polyfit(T0,X,2);
Xpol=7.195*T0.*T0+1.6552*T0-2.9516;
Y=1.0391119*T0-0.3235143;
%Построение n=0 n=1 n=2
X1=2*cos(2*T0)+2*sin(T0*sqrt(2));
figure(1);
plot(T0, X, 'r');
hold on;
plot([-1 1], [P0 P0], 'g');
hold on;
plot(T0, Y, 'b');
hold on;
plot(T0, Xpol, 'k');
hold on;
legend('Сигнал','Полином 0-го порядка','Полином 1-го порядка','Полином 2-го порядка');
%Эмперическая
e0=X-P0;
e1=X-Y;
e2=X-Xpol;
g1=max(abs(e2));
g2=max(abs(e1));
g3=max(abs(e0));
%Графики погрешностей
figure(2);
plot(T0,e0,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot([-1 1], [P0 P0], 'b');
hold on;
plot([-1 1], [g3 g3], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 0-го порядка','max эмп. погрешности 0-го порядка');
grid on;
hold on;
%
figure(3);
plot(T0,e1,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot(T0, Y, 'b');
hold on;
plot([-1 1], [g2 g2], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 1-го порядка','max эмп. погрешности 1-го порядка');
grid on;
hold on;
%
figure(4);
plot(T0,e2,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot(T0, Xpol, 'b');
hold on;
plot([-1 1], [g1 g1], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 2-го порядка','max эмп. погрешности 2-го порядка');
grid on;
hold on;
%ГрафоАнал Метод
em=max(X)-max(Y);
%Полином наилучшего приблежения
Taun=27*pi; %Нижняя граница ИО
Tauv=28*pi; %Верхняя граница ИО
dt=0.0015;
%Масштабирование
tau = Taun:dt:Tauv;
T0 = 2 * (tau - Taun) / (Tauv - Taun) - 1;
T0kv=(2 * (tau - Taun) / (Tauv - Taun) - 1).^2;
X = 2 * cos (2 * ((T0+ 1) * (pi) / 2 + Taun)) + 2 * sin(sqrt(2) * ((T0 + 1) * (pi) / 2 + Taun));
XT0=X.*(2 * (tau - Taun) / (Tauv - Taun) - 1);
P0 = (max(X)-abs(min(X)))/2;
coef=polyfit(T0,X,2);
Xpol=7.195*T0.*T0+1.6552*T0-2.9516;
Y=1.0391119*T0-0.3235143;
%Построение n=0 n=1 n=2
X1=2*cos(2*T0)+2*sin(T0*sqrt(2));
figure(1);
plot(T0, X, 'r');
hold on;
plot([-1 1], [P0 P0], 'g');
hold on;
plot(T0, Y, 'b');
hold on;
plot(T0, Xpol, 'k');
hold on;
legend('Сигнал','Полином 0-го порядка','Полином 1-го порядка','Полином 2-го порядка');
%Эмперическая
e0=X-P0;
e1=X-Y;
e2=X-Xpol;
g1=max(abs(e2));
g2=max(abs(e1));
g3=max(abs(e0));
%Графики погрешностей
figure(2);
plot(T0,e0,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot([-1 1], [P0 P0], 'b');
hold on;
plot([-1 1], [g3 g3], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 0-го порядка','max эмп. погрешности 0-го порядка');
grid on;
hold on;
%
figure(3);
plot(T0,e1,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot(T0, Y, 'b');
hold on;
plot([-1 1], [g2 g2], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 1-го порядка','max эмп. погрешности 1-го порядка');
grid on;
hold on;
%
figure(4);
plot(T0,e2,'g');
hold on;
plot(T0, X, 'r');
hold on;
plot(T0, Xpol, 'b');
hold on;
plot([-1 1], [g1 g1], 'k');
hold on;
legend('Эмперическая погрешность','Сигнал','Полином 2-го порядка','max эмп. погрешности 2-го порядка');
grid on;
hold on;
%ГрафоАнал Метод
em=max(X)-max(Y);
Соседние файлы в папке курсовая работа