clc clear close all % Plot of Some Signals t = -2:0.01:2; figure; plot(t, ustep(t), 'k', 'LineWidth', 2); axis([-2 2 -0.25 1.25]); xlabel('t', 'FontSize', 13, 'FontWeight', 'Bold'); ylabel('ustep(t)', 'FontSize', 13, 'FontWeight', 'Bold'); figure; plot(t, rect(t), 'k', 'LineWidth', 2); axis([-2 2 -0.25 1.25]); xlabel('t', 'FontSize', 13); ylabel('rect(t)', 'FontSize', 13); figure; plot(t, ramp(t), 'k', 'LineWidth', 2); axis([-2 2 -0.25 2.25]); xlabel('t', 'FontSize', 13); ylabel('ramp(t)', 'FontSize', 13); figure; plot(t, sawtooth(t), 'k', 'LineWidth', 2); axis([-2 2 -0.25 1.25]); xlabel('t', 'FontSize', 13); ylabel('sawtooth(t)', 'FontSize', 13); figure; plot(t, triang(t), 'k', 'LineWidth', 2); axis([-2 2 -0.25 1.25]); xlabel('t', 'FontSize', 13); ylabel('triang(t)', 'FontSize', 13); t = -2:0.01:4; figure; plot(t, expfun(t), 'k', 'LineWidth', 2); axis([-2 4 -0.25 1.25]); xlabel('t','FontSize', 13); ylabel('expfun(t)', 'FontSize', 13); t = -4:0.01:4; figure; plot(t,sinc(t), 'k', 'LineWidth', 2); axis([-4 4 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('sinc(t)', 'FontSize', 13); t = -4:0.01:4; figure; plot(t,hsinc(t), 'k', 'LineWidth', 2); axis([-4 4 -1 1]); xlabel('t', 'FontSize', 13); ylabel('hsinc(t)', 'FontSize', 13); figure; plot(t,dmpsin(t), 'k', 'LineWidth', 2); axis([-4 4 -1 1]); xlabel('t', 'FontSize', 13); ylabel('dampsin(t)', 'FontSize', 13); figure; plot(t,dmpcos(t), 'k', 'LineWidth', 2); axis([-4 4 -1 1.25]); xlabel('t', 'FontSize', 13); ylabel('dampcos(t)', 'FontSize', 13); figure; plot(t,funa(t), 'k', 'LineWidth', 2); axis([-4 4 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(t)', 'FontSize', 13); % Even Part and Odd Part fe = @(t) (funa(t) + funa(-t))/2; fo = @(t) (funa(t) - funa(-t))/2; figure; plot(t, fe(t), 'k', 'LineWidth', 2); axis([-4 4 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('fe(t)', 'FontSize', 13); figure; plot(t, fo(t), 'k', 'LineWidth', 2); axis([-4 4 -0.5 0.5]); xlabel('t', 'FontSize', 13); ylabel('fo(t)', 'FontSize', 13); % Signal Transformations t = -6:0.01:6; figure; plot(t, funa(2*t), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(2*t)', 'FontSize', 13); figure; plot(t, funa(t/2), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(t/2)', 'FontSize', 13); figure; plot(t, funa(-t), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(-t)', 'FontSize', 13); figure; plot(t, funa(t-2), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(t-2)', 'FontSize', 13); figure; plot(t, funa(t+2), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(t+2)', 'FontSize', 13); figure; plot(t, funa(-2*t+6), 'k', 'LineWidth', 2); axis([-6 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(-2*t+6)', 'FontSize', 13); t = -10:0.01:6; figure; plot(t, funa(-t/2-2), 'k', 'LineWidth', 2); axis([-10 6 -0.5 1.25]); xlabel('t', 'FontSize', 13); ylabel('funa(-t/2-2)', 'FontSize', 13); % Zero Input Response yo = dsolve('D2y + 3*Dy + 2*y = 0', 'Dy(0) = 1', 'y(0) = 1'); yo = dsolve('D2y + 4*Dy + 13*y = 0', 'Dy(0) = 12', 'y(0) = 6'); % Unit Impulse Response h = dsolve('D2y + 3*Dy + 2*y = 0', 'Dy(0) = 1', 'y(0) = 0'); h = dsolve('D2y + 4*Dy + 13*y = 0', 'Dy(0) = 1', 'y(0) = 0'); h = dsolve('Dy + 6*y = 0', 'y(0) = 1');