Факультет энергетики

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

Новосибирск, 2013

Задача№1

Построить график функции y(x) заданной вариантом (приложение 1). На графике подписать оси и нанести сетку.

Текст программы:

clear all % очистить workspace

clc % очистить command window

x=-10: 0.0001: 10;

y=x.^2-0, 75*x-19; %функция

plot(x, y, 'g') % построение графика

xlim([-10 10]) % лимит по х

ylim([-20 20]) % лимит по у

grid % сетка

xlabel('Ось х') % подписать ось х

ylabel('Ось у') % подписать ось у

title('Парабола') % название графика

legend(' y= x.^2-0, 75*x-19 ') % обозначение графика

График:

Задача№2

Построить график кусочной функции заданной вариантом (приложение 2). Все кривые на графике должны быть одного цвета. Подписать оси и нанести сетку.

Текст программы с find:

clear all % очистить workspace

clc % очистить command window

t=-5: 0.001: 5; % массив с шагом

n=find(t< =-pi); % выбор значений по условию

y(n)=cos(t(n)); % функция графика

n=find(t> -pi & t< =0); % выбор значений по условию

y(n)= t(n)/pi; % функция графика

n=find(t> 0 & t< =1); % выбор значений по условию

y(n)=0; % функция графика

n=find(t> 1 & t< =2); % выбор значений по условию

y(n)= -(t(n)-1).^1/pi; % функция графика

n=find(t> pi); % выбор значений по условию

y(n)= cos(t(n)-1); % функция графика

plot(t, y) % построение графика

grid % сетка

xlabel('Ось х') % подписать ось х

ylabel('Ось у') % подписать ось у

Текст программы с hold:

clc %Очистить workspace

clear all %очистка command window

t=-5: 0.00001: -pi;

y= cos(t);

plot (t, y)

hold on

t=-pi: 0.00001: 0;

y=t/pi;

plot (t, y)

t=0: 0.00001: 1;

y=0;

plot (t, y)

t=1: 0.00001: 2;

y=-((t-1)/pi);

plot (t, y)

t=pi: 0.00001: 5;

y= cos(t-1);;

plot (t, y)

grid

xlabel ('Ось Х') %подпись х

ylabel ('Ось У') %подпись у

title ('График кусочной функции') %название графика

hold off

Задача № 3

Используя команды plot3, mesh и surf построить график функции заданной вариантом (приложение 3).

clc %Очистить workspace

clear all %очистка command window

[x, y]=meshgrid(-2:.1: 2, -2:.1: 2); % массив с шагом

z= cos(3*x).*cos(3*y).*exp(-x^2-y^2); % функция

plot3(x, y, z)

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

clc %Очистить workspace

clear all %очистка command window

[x, y]=meshgrid(-2:.1: 2, -2:.1: 2); % массив с шагом

z= cos(3*x).*cos(3*y).*exp(-x^2-y^2);; % функция

surf(x, y, z)

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

clc %Очистить workspace

clear all %очистка command window

[x, y]=meshgrid(-2:.1: 2, -2:.1: 2); % массив с шагом

z= cos(3*x).*cos(3*y).*exp(-x^2-y^2); % функция

mesh(x, y, z)

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

Задача№4

Используя команду sublot построить на отдельных графиках кривые:

-синусоидальный сигнал частотой 50Гц;

-затухающий синусоидальный сигнал частотой кратной 50Гц (гармонику). Частота и амплитуда сигнала, а также коэффициент затухания заданы вариантом (приложение 4).

clc %Очистить workspace

clear all %очистка command window

subplot(3, 1, 1)

x=0: 1e-4: 0.1;

y=sin(2*pi*50*x);

plot(x, y)

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

legend('синусоидальный сигнал')

subplot(3, 1, 2)

x1=0: 1e-4: 0.1;

y1=1.5*sin(2*pi*5*50*x1).*exp(-40*x1);

plot(x1, y1, 'k')

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

legend('затухание сигнал')

subplot(3, 1, 3)

x2=0: 1e-4: 0.1;

y2=sin(2*50*pi*x1)+1.5*sin(2*pi*5*50*x2).*exp(-40*x2);

plot(x2, y2, 'k')

grid on % сетка

xlabel('x') % подписать ось х

ylabel('y') % подписать ось у

zlabel('z') % подписать ось z

legend('затухание синусоидальный сигнал')