%% comparison of PM and FM (time domain)
f_i=1e3; u_s=2; % information
k_PM=1;k_FM=k_PM*2*pi*f_i; % modulator constants
% time spans of functions
t_end=t_f1+t_f2+t_f3+t_f4;
% partial functions (symbolic and numeric)
f1=@(t) zeros(size(t));s1=f1(t1);
f2=@(t) u_s*(t)/t_f2; s2=f2(t2);
f3=@(t) u_s*ones(size(t));s3=f3(t3);
f4=@(t) u_s*cos(2*pi*f_i*t); s4=f4(t4);
% information signal (numeric)
% integrated functions (symbolic and numeric)
g1=int(f1,t);g1=double(subs(g1,t,t1));
g2=int(f2,t);g2=double(subs(g2,t,t2))+g1(end);
g3=int(f3,t);g3=double(subs(g3,t,t3))+g2(end);
g4=int(f4,t);g4=double(subs(g4,t,t4))+g3(end);
% integrated signal (numeric)
y_T=u_c*cos(2*pi*f_c*tt); % carrier
y_PM=u_c*cos(2*pi*f_c*tt+k_PM*s); % PM-signal
y_FM=u_c*cos(2*pi*f_c*tt+k_FM*g); % FM-signal
disp('comparison of FM and PM');
plot(tt,y_T,'LineWidth',1,'LineStyle',':','Color',[0.5,0.5,0.5])
plot(tt,s,'LineWidth',1,'Color',[0.9,0.3,0])
plot(tt,y_PM,'LineWidth',1,'Color',[0,0.3,0.6])
legend('carrier','information','PM-signal','location','southeast')
plot(tt,y_T,'LineWidth',1,'LineStyle',':','Color',[0.5,0.5,0.5])
plot(tt,s,'LineWidth',1,'Color',[0.9,0.3,0])
plot(tt,g,'LineWidth',1,'Color',[0.9,0.3,1])
plot(tt,y_FM,'LineWidth',1,'Color',[0,0.3,0.6])
legend('carrier','information','information (integrated)','FM-signal','location','southeast')