%=========================================================================
% 信號的時域運算  程序名：product.m
% 本程序給出兩個信號之間的加法、乘法運算，并畫出相應的波形圖
% 設兩個信號分別為：f1(t)=(-t+4)[u(t)-u(t-4)],f2(t)=sin(2*pi*t)
% 畫出信號f3(t)=f1(t)+f1(-t)、f5(t)=f2(t)*f3(t)、f6=f1(t)*f2(t)的波形圖
%=========================================================================
syms t                            %定義符號變量t
f1=sym('(-t+4)*(u(t)-u(t-4))');  %定義符號表達式（方程）f1
subplot(2,3,1)
ezplot(f1,[-5,5])                 %繪制符號方程f1的波形
f2=sym('sin(2*pi*t)');            %定義符號表達式（方程）f2
subplot(2,3,2);
ezplot(f2,[-5,5])
subplot(2,3,5);
ezplot(f2,[-5,5])

y1=subs(f1,t,-t);                  %將信號f1反褶得到信號y1
f3=f1+y1;                         %將信號f1和y1相加得到信號f3，其波形為三角波
subplot(2,3,4);
ezplot(f3,[-5,5])                 %繪制三角波的波形圖
f5=f2*f3;
subplot(2,3,6);
ezplot(f5,[-5,5])
f6=f1*f2;
subplot(2,3,3),ezplot(f6,[-5,5])
%==========================================================================