Математичне моделювання

Міністерство освіти і науки України

Національний технічний університет України

"Київський політехнічний інститут"

 

Приладобудівний факультет

Кафедра приладів і систем орієнтації та навігації

 

Математичне моделювання

Фізичних процесів

Звіт з комп’ютерного практикуму

 

Виконав: студент групи ПГ-11, ПБФ,

Попов О. С.

Перевірив:

Сопілка Ю. В.

 

 

Київ, 2013

Домашнє завдання №1

Fs=1000; t=-10:1/Fs:10; A1=((-1)^i)*(2+0.1*i), A2=3+((-1)^i)*(0.1*i),

Домашнє завдання №2

Fs=100; t=0:1/Fs:8; A1=2+0.1*i, A2=3-0.05*i, A3=1+0.1*i, x1=1+0.1*i, x2=3+0.1*i, x3=5+0.1*i,

Лабораторна робота №1

Fs=1000; %частота дискретизації

t=-0.04:1/Fs:0.04;

T=0.02;

A=5; %амплітуда

s=-A*rectpuls(t+T/2,T)+A*rectpuls(t-T/2,T);

plot (t,s), grid;

title('Приклад застосування функції RECTPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

ylim([-6 6])

 

 

 

 

Лабораторна робота №2

t=0:0.01:10;

y=0.75*rectpuls(t-3,2)+0.5*rectpuls(t-8,0.4)+1.35*rectpuls(t-5,0.8);

plot(t,y), grid;

title('Приклад застосування функції RECTPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

ylim([-0.05 1.4])

 

 

 

Лабораторна робота №3

Fs=1000;

t=-0.05:1/Fs:0.05;

T1=0.02;

T2=0.06;

A=10;

s=A*(T2*tripuls(t,T2)-T1*tripuls(t,T1))/(T2-T1);

plot (t,s), grid;

title('Приклад застосування функції TRIPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

ylim([-0.5 10.5])

 

 

 

Лабораторна робота №4

t=0:0.01:6;

y=0.75*tripuls(t-1,0.5)+0.5*tripuls(t-5,0.5,-1)+1.35*tripuls(t-3,0.8,1);

plot(t,y), grid;

title('Приклад застосування функції TRIPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

ylim([-0.05 1.4])

 

 

Лабораторна робота №5

Fs=100, t=0:1/Fs:10,

y=0.75*gauspuls(t-3,1,0.5),

plot (t,y), grid;

title('Приклад застосування функції GAUSPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

Лабораторна робота №6

Fs=16000; %частота дискретизації

t=-0.006:1/Fs:0.006;

Fc=4000;

bw=0.1;

bwr=-20;

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

plot (t,S), grid;

title('Приклад застосування функції GAUSPULS');

xlabel('Час (с)'); ylabel('Процес Y(t)');

 

 

Лабораторна робота №7

t=0:0.01:50;

y1=0.7*sinc(pi*(t-25)/5);

plot (t,y1), grid;

title('Функція y(t)=0.7*sinc(pi*(t-25)/5)');

xlabel('Час (с)'); ylabel('Процес Y(t)');

 

 

 

Лабораторна робота №8

t=0:0.01:50;

d=[0:50/5:50];

y=0.75*pulstran(t,d,'rectpuls',3);

plot (t,y), grid;

title('Функція y(t)=0.75*pulstran (t,d,''rectpuls'',3)');

xlabel('Час (с)'); ylabel('Вихідний процес Y(t)');

 

 

 

Лабораторна робота №9

t=0:0.01:50;

d=[0:50/5:50];

y=0.7*pulstran(t,d,'tripuls',5,-0.5);

plot (t,y), grid;

title('Функція y(t)=0.75*pulstran (t,d,''tripuls'',5,-0.5)');

xlabel('Час (с)'); ylabel('Вихідний процес Y(t)');

 

 

Лабораторна робота №10

Fs=1000;

t=-0.1:1/Fs:0.5;

tau=0.02; %тривалість імпульсу

d=[20 60 160 260 360]*0.001;%затримка імпульсів

y=pulstran(t,d,'tripuls',tau);

plot (t,y), grid;

title('Функція y=pulstran (t,d,''tripuls'',tau)');

xlabel('Час (с)'); ylabel('Вихідний процес Y(t)');

 

 

Лабораторна робота №11

t=0:0.01:50;

d=[0:50/5:50'];

y=0.7*pulstran(t,d,'gauspuls',1,0.5);

plot (t,y),grid;

title('Y(t)=0.7*pulstran(t,d,''gauspuls'',1,0.5)');

xlabel('Час (с)'); ylabel('Вихідний процес Y(t)');

 

 

Лабораторна робота №12

tau=0.06; %продовжність імпульса t0=0:1/Fs0:tau; %дискретний час для імпульса s0=sin(pi*t0/tau).^2; %вектор сигналів імпульсу

Лабораторна робота №13

t=0:0.01:50;

y1=0.7*sin(pi*t/5);

plot (t,y1), grid;

title('Гармонічні коливання Y(t)');

xlabel('Час (с)'); ylabel('Вихідний Процес Y(t)');

 

Лабораторна робота №14

t=0:0.01:50;

y=0.7*square(pi*t/5,40);

plot (t,y), grid;

title('Прямокутні хвилі Y(t)=0.7*square(pi*t/5,40) ');

xlabel('Час (с)'); ylabel('Процес Y(t)');

 

 

Лабораторна робота №15

Fs=1000;%частота дискретизації

t=-0.01:1/Fs:0.05;

A=3; f0=50; tau=0.05;

S=(square(2*pi*t*f0,f0*tau*100+1)*A/2);

plot(t,S),grid;

title('Приклад застосування функції sguare');

xlabel('Час (с)'); ylabel('Вихідний процес Y(t)');

xlim([-0.025 0.125])