курсовая работа / kurs_oulsh

clear();
%Описание сигнала по Уолшу
Taun=27*pi; %Нижняя граница ИО
Tauv=28*pi; %Верхняя граница ИО
dt=0.0015;
%Масштабирование
tau = Taun:dt:Tauv;
tn=0.00001:dt:1;
T0 = (tn ) * (Tauv - Taun) + Taun;
T0_alt = 2 * (tau - Taun) / (Tauv - Taun) - 1;
X = (2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)));
%Описание
r1=(sign(sin(2*pi*tn)));
r2=(sign(sin(4*pi*tn)));
r3=(sign(sin(8*pi*tn)));
walt0=1;
walt1=r1;
walt2=r1.*r2;
walt3=r2;
walt4=r2.*r3;
walt5=r1.*r2.*r3;
walt6=r1.*r3;
walt7=r3;
%Коэфы обэмы
C0fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun))))*walt0);
C0 =integral(C0fun,0,1);

C1fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(2*pi*tn))));
C1 =integral(C1fun,0,1);

C2fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(2*pi*tn))).*(sign(sin(4*pi*tn))));
C2 =integral(C2fun,0,1);

C3fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(4*pi*tn))));
C3 =integral(C3fun,0,1);

C4fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(4*pi*tn))).*(sign(sin(8*pi*tn))));
C4 =integral(C4fun,0,1);

C5fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(2*pi*tn))).*(sign(sin(4*pi*tn))).*(sign(sin(8*pi*tn))));
C5 =integral(C5fun,0,1);

C6fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(2*pi*tn))).*(sign(sin(8*pi*tn))));
C6 =integral(C6fun,0,1);

C7fun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).*(sign(sin(8*pi*tn))));
C7 =integral(C7fun,0,1);

F=(C0*walt0+C1*walt1+C2*walt2+C3*walt3+C4*walt4+C5*walt5+C6*walt6+C7*walt7);
%Погрешности
e=X-F;
g=max(abs(e));
sumC=C0*C0+C1*C1+C2*C2+C3*C3+C4*C4+C5*C5+C6*C6+C7*C7;
Ffun=@(tn) (((2 * cos (2*((tn ) * (Tauv - Taun) + Taun)) + 2 * sin(sqrt(2)*((tn ) * (Tauv - Taun) + Taun)))).^2);
Fail=integral(Ffun,0,1);
Ekv=Fail-sumC;

sqrtEkv=sqrt(Ekv);

figure(1);
plot(tn, X, 'g');
hold on;
plot(tn, F, '--k');
legend('Сигнал','Описание по Оулшу');
hold on;
figure(2);
plot(tn, e, 'g');
hold on;
plot([0 1], [g g], 'r');
hold on;
plot([0 1], [sqrtEkv sqrtEkv], 'b');
hold on;
legend('Эмперическая погрешность','max эмп. погрешности','Корень СКВ');
hold on;