对锯齿波信号频谱分析，并用两个非原始信号周期整倍数的正弦波信号进行双边带调制
close all clear Fs=100000; %采样频率 t=0:1/Fs:0.1; x1=sawtooth(1000/7*2*pi*t); %信号周期2pi/(t的系数) figure plot(t,x1,'LineWidth',1.2) axis([0 0.025 -1 1]) xlabel('时间 (ms)') ylabel('幅度') title('原始信号(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on %傅立叶变换 x1=x1-mean(x1); %消去直流分量，使频谱更能体现有效信息 N=10001; %即采样个数。其实就是length(x1); z=fft(x1); %频谱分析 f=(0:N-1)*Fs/N; Mag=2*abs(z)/N; %幅值，单位同被测变量x1 Pyy=Mag.^2; %能量；对实数系列X，有 X.*X=X.*conj(X)=abs(X).^2=X.^2，故这里有很多表达方式 figure plot(f,Pyy,'LineWidth',1.2) %显示频谱图 axis([0 2000 0 0.35]) xlabel('频率 (Hz)') ylabel('幅度') title('频谱图(频域)') set(gca,'FontSize',12,'FontWeight','bold') grid on; x2=cos(700/7*2*pi*t); figure plot(t,x2,'LineWidth',1.2) axis([0 0.025 -1 1]) xlabel('时间 (ms)') ylabel('幅度') title('余弦信号a(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on x3=cos(1100/7*2*pi*t); figure plot(t,x3,'LineWidth',1.2) axis([0 0.025 -1 1]) xlabel('时间 (ms)') ylabel('幅度') title('余弦信号b(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on x12=x1.*x2; figure plot(t,x12,'LineWidth',1.2) xlabel('时间 (ms)') ylabel('幅度') title('已调信号a(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on %傅立叶变换 x12=x12-mean(x12); %消去直流分量，使频谱更能体现有效信息 z12=fft(x12); %频谱分析 f12=(0:N-1)*Fs/N; Mag12=2*abs(z12)/N; %幅值，单位同被测变量x1 Pyy12=Mag12.^2; %能量；对实数系列X，有 X.*X=X.*conj(X)=abs(X).^2=X.^2，故这里有很多表达方式 figure plot(f12,Pyy12,'LineWidth',1.2) %显示频谱图 axis([0 2000 0 0.1]) xlabel('频率 (Hz)') ylabel('幅度') title('已调信号a频谱图(频域)') set(gca,'FontSize',12,'FontWeight','bold') grid on; x13=x1.*x3; figure plot(t,x13,'LineWidth',1.2) xlabel('时间 (ms)') ylabel('幅度') title('已调信号b(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on %傅立叶变换 x13=x13-mean(x13); %消去直流分量，使频谱更能体现有效信息 z13=fft(x13); %频谱分析 f13=(0:N-1)*Fs/N; Mag13=2*abs(z13)/N; %幅值，单位同被测变量x1 Pyy13=Mag13.^2; %能量；对实数系列X，有 X.*X=X.*conj(X)=abs(X).^2=X.^2，故这里有很多表达方式 figure plot(f13,Pyy13,'LineWidth',1.2) %显示频谱图 axis([0 2000 0 0.12]) xlabel('频率 (Hz)') ylabel('幅度') title('已调信号b频谱图(频域)') set(gca,'FontSize',12,'FontWeight','bold') grid on x23=x12+x13; figure plot(t,x23,'LineWidth',1.2) xlabel('时间 (ms)') ylabel('幅度') title('合成信号(时域)') set(gca,'FontSize',12,'FontWeight','bold') set(gca,'xticklabel',get(gca,'xtick')*1e3) grid on %傅立叶变换 x23=x23-mean(x23); %消去直流分量，使频谱更能体现有效信息 z23=fft(x23); %频谱分析 f23=(0:N-1)*Fs/N; Mag23=2*abs(z23)/N; %幅值，单位同被测变量x1 Pyy23=Mag23.^2; %能量；对实数系列X，有 X.*X=X.*conj(X)=abs(X).^2=X.^2，故这里有很多表达方式 figure plot(f23,Pyy23,'LineWidth',1.2) %显示频谱图 axis([0 2000 0 0.12]) xlabel('频率 (Hz)') ylabel('幅度') title('合成信号频谱图(频域)') set(gca,'FontSize',12,'FontWeight','bold') grid on
