如何用matlab模拟简单的信号调制，初学求代码

[X,FS,NBITS]=wavread('春天里')

N= length(X)

sound (X,FS)

W=2\N*[0:N-1]%连续谱的数字角频率

FY=fft(X,N)

figure

subplot(211)plot(X)title('原音乐信号波形')

subplot(212)plot(W,abs(FY))title('原音乐信号频谱')

%减抽样

D=12j=0

for i=1:D:length(X/20),

j=j+1

X1(j)=X(i)

end

sound(X1,FS/12)

N1=1024

W1=2\N1*[0:N1-1]

F1Y=fft(X1,N1)

figure

subplot(211)plot(X1)title('减抽样后信号波形')

subplot(212)plot(W1,abs(F1Y))title('减抽样后信号频谱')

%AM调制

n=0:N-1

x=cos(n*pi*0.8)%调制信号

Nt=length(x)

FtY=fft(x,Nt)

Wt=2/Nt*[0:Nt-1]

y=X.*x'%对信号进行调制

N2=length(y)

F2Y=fft(y,N2)

W2=2/N2*[0:N2-1]

sound(y,FS)

figure

plot(Wt,abs(FtY))

grid on

title ('调制信号cos的频谱图')

figure

subplot(2,1,1)plot(W,abs(FY))

grid on

title ('原音乐信号的频谱')

subplot(2,1,2)plot(W2,abs(F2Y))

grid on

title ('已调信号频谱')

最近正在学习这个东西，希望有所帮助

#include<reg51.h>

sbit P1_0 = P1^0 //定义P10脚名称为P1_0

sbit P1_1 = P1^1

sbit P1_2 = P1^2 //定义P11脚名称为P1_1

delay(unsigned char ms) //毫秒级延时函数,ms<255

{

unsigned char i,j

for(i=msi>0i--)

for(j=121j>0j--)

}

usdelay(unsigned char us)//微秒级延时函数,us<255

{

for(us>0us--)

}

main()

{

P1_2 = 0

while(1)

{

P1_0 = 0 delay(1) P1_0 = 1

P1_1 = 0 usdelay(10) P1_1 = 1

}

}

//P12是用来对比的，你可以在P10，P11，P12各用一个LED来检测效果

u=(1+m).*c%调制信号

%先产生高斯噪声，方差为1

n = randn(size(u))

% 确定噪声标准差，控制噪声的大小，如0.01，0.1等

sigmaU = 0.1

u = u + n*sigmaU

%=== === == ==消息信号的频谱

blabla

---