37. Радиоимпульс пен оның спектрінің графигін алыңыз.

Гаустық орайжанауыш және бірлік амплитудалы жалғыз ралиоимпульсті құру үшін gauspuls функцияы қолданылады: y=gauspuls(t,fc,bw,bwr) мұндағы fc тасымалдаушы жиілік, bw спектрдің салыстырмалы ені, bwr спектрдің ені б/ша өлшеу жүргізілетін децибелмен өлшенетін деңгей. У векторы келесі формуламен есептеледі: у=exp(at²)cos(2πf_ct)

а коэффиценті импульс ұзақтығын, енін де анықтайды.

Gauspuls функциясын шақырғанда, бірден үшке дейін шығыс параметрлерін қолдануға болады: [y,yq,ye]=gauspuls(t)

Yq анықталған у радиоимпульстің квадратуралық қосымшасын береді. Yq векторы у векторынан 90 градусқа тасымалдаушы тербелістің ығысуымен ажыратылады. Ye пайда болған радиоимпульстің орайжанауышын береді.

Мысал. 4 кГц тасымалдаушы жиілігімен, -20 дБ деңгейімен, өлшенген спектрдің салыстырмалы ені 10% радиоимпульстің және оның спектрінің графиктерін салайық. Дискреттеу жиілігі = 16 кГц.

Fs=16e3; t=-10e-3:1\Fs:10e-3; Fc= 4e3; bw=0.1; bwr=-20;

s=gauspuls(t,Fc,bw,bwr);

N=2^nextpow2(length(s));

Spectra=fft(s,N);

spectradB=20*log10(abs(spectr));

f=(0:N-1)\N*Fs;

subplot(2,1,1);

plot(t,s,`k`);

title(`gaupuls`);

subplot(2,1,2); plot(f(1:N\2)

spectradB(1:N\2),`k`);

maxspectra=20*log10(max(abs(spectra)));

freq=Fc*[1-bw\2 1+bw\2];

hold on, plot (freq, maxspectra ([1 1]) +bwr, `k+`);

title(`Spectra of gaupuls`);

38. Матрицалық ойындар есебіне мысал келтіріңіз.

	Kizil tuz	Kara tuz	Kizil ekilik
Kizil tuz	1	-1	-2
Kara tuz	-1	1	1
Kizil ekilik	2	-1	0

>>A=[1 -1 2;-1 1 1;2 -1 0]; %tikelei esepke habarlasu;

>>[x,fval]=linprog(ones(3,1),-A',-ones(3,1),[],[],zeros(3,1))

x = 0.0000

3.0000

2.0000

fval = 5.0000

>>% keri esepke habarlasu:

>>[y,gval]=linprog(-ones(3,1),A',ones(3,1),[],[],zeros(3,1))

y =

0.0000

1.0000

1.0000

gval =-2.0000

>> v=1./fval

v =0.2000

>> p=v*x

p =0.0000

0.6000

0.4000

>> q=v*y

q = 0.0000

0.2000

0.2000

39. Математикалық тербелістердің еркін тербелістерін алыңыз. Еріксіз тербелістер.

Тербелістер туралы есептер физиканың барлық салаларында кездеседі. Көп жағдайларда мүлдем өзгеше физикалық обьектілердің тербелістері бір біріне ұқсас. Қарапайым мысалдар маятниктің шағын тербелістері және конденсатор мен шарғыдан тұратын тізбектегі электр тербелістері.

Ұзындығы L қатты сырықтың бір шетіне қадалған, екінші шеті ілу нүктесінде бекітілген, массасы m жүктің қозғалысын қарастырайық. Жүк центрі о нүктесіндегі L радиусты шеңбердің доғасы бойынша қозғалады, жүктің орны сырықтың вертикальдан ауытқу бұрышымен сипатталады. Программаның мәтінімен салыстыруға ыңғайлы болу үшін маятниктің вертикальды бағыттан ауытқу бұрышын х пен белгілейміз.

Маятниктің аз ауытқу қозғалысы белгілібұл гармониялық тербелістер. Маятниктің қозғалыс заңы келесі түрде жазылады: х=acos(ω_ot+φ_o).

Егер ауытқулар аз болмаса, тербелістерді сипаттайтын теңдеу ẍ +ω²sinx=o. Соңғы теңдеу бейсызық, оның шешімі, басқа көп бейсызық теңдеулердің шешімдері сияқты, элементар функциялар жазылмайды, сондықтан оның сандық шешімін алу қажет. Мейлінше аз ауытқу бұрышы кезінде теңдеу сызықты болады ẍ +ω²sinx=o.

function F=mayat(t,y);

global omega;

F=[y(2); - omega^2*sin(y(1))];

global omega;

x0=[2 3];

for omega=1:0.3:4

[T,Y]=ode45('mayat',[0:0.05:4*pi],x0);

hold on

plot(Y(:,1),Y(:,2),'k');

grid on

axis([-5 20 -8 8]);

title('mayatnik')

end

Еріксіз т/стер. Маятникке горизонталь бағытталған F(t)=Bcoswt сырткы айнымалы күш әсер еткенде жүйе мәжбүр т/с д.а. x(2)+w²sinx=Bcoswt.

Егер маятник бастапкы уакытта тыныштыкта турса, оган периодты куш әсер етсе, онда ол біраз уақыт теңселіп тұрады, регулярлы емес козгалыста б/ды. Бул кезде амплитуданың тек өсуі ғана емес, периодты үдеуі не азаюы мүмкін. Біраз уакыт өткеннен кейін козгалыс калыптасады, сырткы куш w жиілікпен периодты б/ды.