Question

哪位大神会MATLAB的，给下面的程序每个语句注释一下呗

clear all;
close all;
clc;
t0=0.15;
dt=0.0001;
t=[0:dt:1];
fc=250;
phi0=0;
ct=cos(2*pi*fc*t+phi0); %载波信号
fm=10;
fun = @(t)1.*(t>=0&t<=t0/3)+(-2).*(t>=t0/3&t<=2*t0/3)+0.*(t>=2*t0/3);
mt=fun(t);
mh=imag(hilbert(mt));
b=sin(2*pi*fc.*t)
ussb=mt.*ct-imag(hilbert(mt)).*b;
lssb=mt.*ct+imag(hilbert(mt)).*b;
dt=t(2)-t(1); % 采样周期
f=1/dt; % 采样频率(Hz)
X1=fft(ussb); % 计算x的快速傅立叶变换X
X2=fft(lssb);
N=1/dt;
F1=X1(1:N/2+1); % F(k)=X(k)(k=1:N/2+1)
F2=X2(1:N/2+1);
f=f*(0:N/2)/N; % 使频率轴f从零开始
figure
subplot(2,1,1),plot(f,abs(F2)),hold on,title('lssb已调信号频谱图');plot(-f,abs(F2));xlim([-fc*2,fc*2]);
xlabel('Frequency');
ylabel('|F(k)|')

ut=lssb
[m,n]=size(ut);
ni=0.05*randn(m,n);
u0=ut+ni;
w1=2*dt*(fc-2*fm);
w2=2*dt*(fc+2*fm);
[b,a]=butter(4,[w1,w2],'bandpass');
u1=filter(b,a,u0);
figure(2)
u2=u1.*ct;
B=2*fm;
wn3=2*dt*B;
[b,a]=butter(4,wn3,'low');
u3=filter(b,a,u2);
subplot(3,1,1),plot(t,mt),title('mt----原信号')
subplot(3,1,2),plot(t,-u3),title('u3----解调后的信号')
subplot(3,1,3),plot(t,-u3),hold on,plot(t,mt),title('原信号与解调信号的比较')

Answer 1

就是载波信号的调制和解调啊 这类的课程设计 在百度文库有不少 建议你找找

每句都注释 还是等等其他高手把

追问

关键的几句注释就好，特别是第二个程序

Answer 2

真想发给我们老师让他帮你。。

追问

在做课设。。。第二段程序不大明白。。。

追答

我们上学期交的。我老师微博不是经常在。

Answer 3

兄弟，好好看书，这么多句还得靠自己

追问

结果弄出来的，第一个程序基本懂，主要是第二个不怎么懂解释