数字信号处理：用matlab编程实现音乐信号的AM调制

clear all; close all;clc

[x,fs,bit]=wavread('F:\费玉清-一剪梅00_01_23-00_01_28.wav'); size(x);

x=x(:,1);

N=length(x);

sound(x,fs);

figure

subplot(2,1,1);

plot(x)

F1=fft(x,N);

grid on

title('音乐信号时域波型')

xlabel('Time')

ylabel('Magnitude')

w1=2/N*[0:N-1];

subplot(2,1,2);plot(w1,abs(F1))

grid on

title('音乐信号频域波型')

xlabel('Frequency/Hz')

ylabel('Magnitude')

n=0:N-1;

x1=cos(0.5*n);

figure

subplot(2,1,1);

plot(x1)

title('余弦信号时域波型')

xlabel('Time')

ylabel('Magnitude')

F2=fft(x1,N);

w1=2/N*[0:N-1];

subplot(2,1,2);

plot(w1,abs(F2))

grid on

title('余弦信号频域波型')

xlabel('Frequency/Hz')

ylabel('Magnitude')

x2=x.*x1';

figure

subplot(2,1,1);

plot(x2)

F3=fft(x2,N);

grid on

title('高频调制后音乐信号时域波型')

xlabel('Time')

subplot(2,1,2);

plot(w1,abs(F3))

grid on

title('高频调制后音乐信号频域波型') xlabel('Frequency/Hz')

ylabel('Magnitude')