MATLAB相关问题

MATLAB问题可以问我

function st=mstg
%产生信号序列向量 st,并显示 st 的时域波形和频谱
%st=mstg 返回三路调幅信号相加形成的混合信号，长度 N=1600
N=1600 %N 为信号 st 的长度。
Fs=10000;T=1/Fs;Tp=N*T; %采样频率 Fs=10kHz，Tp 为采样时间
t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;
fc1=Fs/10; %第 1 路调幅信号的载波频率 fc1=1000Hz,
fm1=fc1/10; %第 1 路调幅信号的调制信号频率 fm1=100Hz
fc2=Fs/20; %第 2 路调幅信号的载波频率 fc2=500Hz
fm2=fc2/10; %第 2 路调幅信号的调制信号频率 fm2=50Hz
fc3=Fs/40; %第 3 路调幅信号的载波频率 fc3=250Hz,
fm3=fc3/10; %第 3 路调幅信号的调制信号频率 fm3=25Hz
xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第 1 路调幅信号
xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第 2 路调幅信号
xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第 3 路调幅信号
st=xt1+xt2+xt3; %三路调幅信号相加
fxt=fft(st,N); %计算信号 st 的频谱
%====以下为绘图部分，绘制 st 的时域波形和幅频特性曲线
subplot(3,1,1)
plot(t,st);grid;xlabel('t/s');ylabel('s(t)');
axis([0,Tp/8,min(st),max(st)]);title('(a) s(t)的波形')
subplot(3,1,2)
stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b) s(t)的频谱')
axis([0,Fs/5,0,1.2]);
xlabel('f/Hz');ylabel('幅度')
fp=280;fs=450;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;
[N, wp]=ellipord(wp, ws, rp, rs);
[B, A]=ellip(N, rp, rs, wp);
y1t=filter(B,A, st);
figure (2) ;subplot(2, 1, 1);
p1ot (B, A);
yt='y_ 1(t)';
subplot (2, 1, 2);tplot(y1t,T,yt);
他就提示我未定义函数或变量，我初学

大佬，我想问一下如何完成c之后的步骤

x = [1 1 0 1 0 1 1 1 0 1 ]
m=10
P=0.5
% x=2*round(rand(1,m)+P-0.5)-1;
% % x = round(rand(1,m))
figure(5)
stem(x)
t0=300;
t=0:1/t0:length(x);
for i = 1:length(x)
if(x(i) ==1)
for j = 1 : t0/2;
y( t0/2 * (2 * i - 2) + j ) = 1;
y( t0/2 * (2 * i - 1) + j ) = 0;
end
else
for j = 1 : t0/2;
y( t0/2 * (2 * i - 2) + j ) = -1;
y( t0/2 * (2 * i - 1) + j ) = 0;
end
end
end
y = [y , x(i)];
M = max(y);
m = min(y);
figure(1)
plot(t,y);grid on;
axis([0,i,m-0.1,M+0.1])
Ts=1;
N=i;
N_sample=150;
dt = Ts/N_sample;
df = 1.0/(20.0*Ts);
t = -10*Ts:dt:10*Ts;
f = -2/Ts:df:2/Ts;
alpha = 0.2;
for n = 1:length(alpha)
for k = 1 : length(f)
if abs(f(k))>0.5*(1+alpha(n))/Ts
Xf(n,k) = 0;
elseif abs(f(k))<0.5*(1-alpha(n))/Ts
Xf(n,k) = Ts;
else
Xf(n,k) = 0.5*Ts*(1+cos(pi*Ts/(alpha(n))+eps)*(abs(f(k))-0.5*(1-alpha((n))/Ts)));
end
end
xt(n,:) = sinc(t/Ts).*(cos(alpha(n)*pi*t/Ts))./(1-4*alpha(n)^2*t.^2/Ts^2+eps);
end
figure(2);
plot(t,xt)
axis([-10 10 -0.5 1.1])
%
r = conv(y,xt)
figure(3)
plot(r)
这是写的截至c波形的代码