Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Цифровая обработка сигналов. Методическая разработка к лабораторным занятиям для студентов-бакалавров дневного и заочного обучения
.pdf
005 Y3(n)=4*(T(n)/P)*(D-(T(n)/P)^2);
006 end
007 endfunction;
001 function Y5=polinom5(T)
002 D1=0.745966;
003 D3=-0.305500;
004 D5=0.138149;
005 for n=1:nmax
006 Y5(n)=4*(T(n)/P)*(D1+4*(T(n)/P)^2*(D3+D5*(T(n)/P)^2));
007 end
008 endfunction;
058 //Трапецеидальное колебание
059 clf()
060 n=[1:nmax];
061 z=pila(A,P,0,nmax);
062 T=(abs(z)-P/2)*B;
063 for n=1:nmax
064 if T(n)>P
065 Tp(n)=P;
066 elseif T(n)<-P
067 Tp(n)=-P;
068 else Tp(n)=T(n);
069 end
070 end
071 n=[1:nmax];
072 subplot(3,2,1)
073 plot2d3(n,Tp,style=[color("green")])
074 xgrid
075 xtitle('Трапецеидальное колебание')
076 //Определение спектра трапецеидального колебания
077 i0=ceil(2*P/A);
078 zs=zeros(1:i0);
079 for i=1:i0
080 Tps(i)=Tp(i);
081 end
082 STp=dpf(Tps,6,i0);
083 STpo=abs(STp)/abs(STp(2));
084 subplot(3,2,2)
085 i=[1:6];
086 k=i-1;
087 plot2d3(k,STpo,style=[color("green")])
088 xgrid
089 xtitle('Спектр трапецеидального колебания')
090 AmpTp=zeros(1:6);
091 for m=1:6
092 if STpo(m)<0.00001 then
093 AmpTp(m)=0;
094 else
095 AmpTp(m)=STpo(m)*100;
096 end
097 end
098 // Косинусно-синусный генератор
111
099 // Пилообразные колебания, сдвинутые друг
100 // относительно друга на четверть периода
101 zC=pila(A,P,0,nmax);
102 zS=pilaS(A,P,0,nmax);
103 subplot(3,2,3)
104 plot2d(n,zC,style=[color("blue")])
105 plot2d(n,zS,style=[color("red")])
106 xgrid
107 xtitle('Пилообразные колебания zC и zS')
108 // Формирование треугольных колебаний Tc и Ts
109 Tc=P/2-abs(zC);
110 Ts=abs(zS)-P/2;
111 subplot(3,2,4)
112 plot2d(n,Tc,style=[color("blue")])
113 plot2d(n,Ts,style=[color("red")])
114 xgrid
115 xtitle('Треугольные колебания Tc и Ts')
116 // Определение косинусной и синусной компонент
117 // с использованием полиномов третьей и пятой степени
118 C3=polinom3(Tc);
119 S3=polinom3(Ts);
120 C5=polinom5(Tc);
121 S5=polinom5(Ts);
122 C=C3;
123 S=S3;
124 subplot(3,2,5)
125 plot2d(n,C,style=[color("blue")])
126 plot2d(n,S,style=[color("red")])
127 xgrid
128 xtitle('Синусная и косинусная компоненты КСГ')
129 //Определение спектра генерируемых колебаний
130 for i=1:i0
131 Cs(i)=C(i);
132 end
133 SCs=dpf(Cs,6,i0);
134 SCso=abs(SCs)/abs(SCs(2));
135 subplot(3,2,6)
136 i=[1:6];
137 k=i-1;
138 plot2d3(k,SCso,style=[color("blue")])
139 xgrid
140 xtitle('Спектр выходных колебаний КСГ')
141 Amp3=zeros(1:6);
142 for m=1:6
143 if SCso(m)<0.00001 then
144 Amp3(m)=0;
145 else
146 Amp3(m)=SCso(m)*100;
147 end
148 end
149 C=C5;
150 S=S5;
112
151 subplot(3,2,5)
152 plot2d(n,C,style=[color("blue")])
153 plot2d(n,S,style=[color("red")])
154 xgrid
155 //Определение спектра выходных колебаний]
156 //при использовании полинома пятой степени
157 for i=1:i0
158 Cs(i)=C(i);
159 end
160 SCs=dpf(Cs,6,i0);
161 SCso=abs(SCs)/abs(SCs(2));
162 subplot(3,2,6)
163 i=[1:6];
164 k=i-1;
165 plot2d3(k,SCso,style=[color("blue")])
166 xgrid
167 Amp5=zeros(1:6);
168 for m=1:6
169 if SCso(m)<0.00001 then
170 Amp5(m)=0;
171 else
172 Amp5(m)=SCso(m)*100;
173 end
174 end
113
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]