n = 0:1:99; % 产生长度为100的序列 x1 = rand(1,100); % 产生随机序列 x2 = linspace(0,1,100); % 产生线性序列 % 进行相加、乘及卷积运算 y1 = x1 + x2; % 相加 y2 = x1 .* x2; % 相乘 y3 = conv(x1,x2); % 卷积 % 利用fft函数求离散傅里叶变换 N = 1024; % 采样点数 X1 = fft(x1,N); % 对x1进行FFT变换 X2 = fft(x2,N); % 对x2进行FFT变换 % 绘制信号和变换结果的图像 subplot(2,2,1);plot(n,x1);title('x1'); subplot(2,2,2);plot(n,x2);title('x2'); subplot(2,2,3);plot(n,y1);title('x1+x2'); subplot(2,2,4);plot(n,y2);title('x1x2');能把卷积的图像也体现出来吗
时间: 2023-08-05 19:05:31 浏览: 82
当进行卷积运算时,得到的结果序列长度为 N + M - 1,其中 N 和 M 分别为输入序列的长度。因此,在上述代码中,进行卷积运算的结果序列长度应该为 199,而不是 100。为了能够正确地绘制卷积的图像,需要对代码进行如下修改:
```matlab
n_conv = 0:1:198; % 产生长度为199的序列
y3 = conv(x1,x2); % 卷积
% 绘制信号和变换结果的图像
subplot(2,2,1);plot(n,x1);title('x1');
subplot(2,2,2);plot(n,x2);title('x2');
subplot(2,2,3);plot(n,y1);title('x1+x2');
subplot(2,2,4);plot(n_conv,y3);title('x1*x2');
```
这样就可以绘制出卷积的图像了。
相关问题
帮我逐句解释一下这段代码%生成m序列 周期为2^N-1 x1=0;x2=0;x3=1; m=350;%伪随机码的周期 for i=1:m y3=x3;y2=x2;y1=x1; x3=y2;x2=y1; x1=xor(y3,y1); %xor逻辑异运算 l(i)=y1; end for i=1:m M(i)=1-2*l(i); %将单极性m序列变为双极性m序列 end k=1:1:m; ym=fft(M,4096); %傅里叶变换 magm=abs(ym); fm=(1:2048)*200/2048; [c,d]=xcorr(M,'unbiased'); %xcorr计算互相关 %随机生成50位二进制比特序列,进行扩频编码 n=50;a=0; x_rand=rand(1,n); for i=1:n if x_rand(i)>=0.5 x(i)=1;a=a+1; else x(i)=0; end %大于等于0.5取1,小于0.5取0 end t=0:n-1; tt=0:349; L=1:7*n; y(L)=0; y=rectpulse(x,7); %利用矩形脉冲将序列扩展为350 s(L)=0; for i=1:350 %扩频后,码率变为100/7*7=100 s(i)=xor(L(i),y(i)); end tt=0:7*n-1; %对扩频前后信号进行BPSK调制,观察其时域波形 %BPSK调制采用2Khz信号cos(2*2000*t)作为载波 fs=2000; %载频频率 fc=100000 %采样率 T=1/fs; ts=0:0.00001:3.5-0.00001;%为了使信号看起来更光滑,作图时采样频率为100kHz,ps=cos(2*pi*fs*ts) s_b=rectpulse(s,1000); %将冲激信号补成矩形信号 s_bpsk=(1-2.*s_b).*cos(2*pi*fs*ts);%扩频后信号BPSK调制时域波形,(1-2.*s_b)是1,-1序列 s_bb=rectpulse(x,7000); s_bpskb=(1-2.*s_bb).*cos(2*pi*fs*ts);%无扩频信号BPSK调制时域波形
这段代码的作用是生成一个周期为2^N-1的伪随机码,并进行扩频编码和BPSK调制,最终观察其时域波形。
首先,定义了三个变量x1、x2、x3,初始值分别为0、0、1,以及周期m为350。然后通过循环生成伪随机码,其中使用了xor逻辑异运算,将生成的单极性m序列变为双极性m序列。
接着,通过FFT傅里叶变换和xcorr函数计算互相关,得到了伪随机码的频域特性和自相关函数。
然后,使用rand函数随机生成50位二进制比特序列,并进行扩频编码,将序列扩展为350个,再将扩频后的信号进行BPSK调制,最终得到了扩频后信号和无扩频信号的时域波形。其中,采用cos(2*2000*t)作为载波,fs=2000为载波频率,fc=100000为采样率,ts设置了一个时间范围,使用rectpulse函数将信号补成矩形信号,并使用(1-2.*s_b)对信号进行BPSK调制,得到扩频后信号的时域波形。对于无扩频信号,直接使用rectpulse和BPSK调制函数得到其时域波形。
clear all; clc; %% 生成双极性7位M序列 X1 = 0; X2 = 0; X3 = 1; m = 350; % 重复50遍的7位单极性m序列 for i = 1:m Y1 = X1; Y2 = X2; Y3 = X3; X3 = Y2; X2 = Y1; X1 = xor(Y3, Y1); L(i) = Y1; end for i = 1:m M(i) = 1 - 2 * L(i); % 将单极性m序列变为双极性m序列 end %% 画出双极性7位M序列频谱 k = 1:1:m; figure(1); subplot(2, 1, 1); stem(k-1, M); axis([0, 7, -1, 1]); xlabel('k'); ylabel('M序列'); title('双极性7位M序列'); subplot(2, 1, 2); ym = fft(M, 4096); magm = abs(ym); % 求双极性m序列频谱 fm = (1:2048) * 200 / 2048; plot(fm, magm(1:2048) * 2 / 4096); title('双极性7位M序列的频谱'); xlabel('Hz'); %% 生成扩频前待发送二进制信息序列和扩频后序列码 N = 50; a = 0; x_rand = rand(1, N); % 产生50个0与1之间随机数 for i = 1:N if x_rand(i) >= 0.5 % 大于等于0.5的取1,小于0.5的取0 x(i) = 1; a = a + 1; else x(i) = 0; end end t = 0:N-1; figure(2); subplot(2, 1, 1); stem(t, x); title('扩频前待发送二进制信息序列'); tt = 0:349; subplot(2, 1, 2); L = 1:7*N; y = rectpulse(x, 7); s(L) = 0; for i = 1:350 % 扩频后,码率变为100/7*7=100Hz s(i) = xor(L(i), y(i)); end tt = 0:7*N-1; stem(tt, s); axis([0, 350, 0, 1]); title('扩频后的待发送序列码');
% 清空变量并清空命令窗口
clear all;
clc;
%% 生成双极性7位M序列
X1 = 0; X2 = 0; X3 = 1; % 初始化三个寄存器
m = 350; % 重复50遍的7位单极性m序列
% 生成双极性7位M序列
for i = 1:m
Y1 = X1; Y2 = X2; Y3 = X3; % 保存当前寄存器值
X3 = Y2; X2 = Y1; % 移位操作
X1 = xor(Y3, Y1); % 异或操作
L(i) = Y1; % 保存当前寄存器最低位的值
end
% 将单极性m序列变为双极性m序列
for i = 1:m
M(i) = 1 - 2 * L(i);
end
%% 画出双极性7位M序列频谱
k = 1:1:m;
figure(1); % 新建图像窗口
subplot(2, 1, 1); % 画图1
stem(k-1, M); axis([0, 7, -1, 1]); % 设置坐标轴范围
xlabel('k'); ylabel('M序列'); % 设置坐标轴标签
title('双极性7位M序列'); % 设置图像标题
subplot(2, 1, 2); % 画图2
ym = fft(M, 4096); % 对双极性m序列进行FFT
magm = abs(ym); % 求双极性m序列频谱
fm = (1:2048) * 200 / 2048; % 计算频率
plot(fm, magm(1:2048) * 2 / 4096); % 画图
title('双极性7位M序列的频谱'); xlabel('Hz');
%% 生成扩频前待发送二进制信息序列和扩频后序列码
N = 50; % 信息序列长度
a = 0;
x_rand = rand(1, N); % 产生50个0与1之间随机数
% 生成待发送二进制信息序列
for i = 1:N
if x_rand(i) >= 0.5 % 大于等于0.5的取1,小于0.5的取0
x(i) = 1;
a = a + 1;
else
x(i) = 0;
end
end
t = 0:N-1; % 生成时间序列
figure(2); % 新建图像窗口
subplot(2, 1, 1); % 画图1
stem(t, x); % 画图
title('扩频前待发送二进制信息序列');
tt = 0:349; % 生成时间序列
subplot(2, 1, 2); % 画图2
L = 1:7*N; y = rectpulse(x, 7); % 7倍扩频
s(L) = 0;
for i = 1:350 % 扩频后,码率变为100/7*7=100Hz
s(i) = xor(L(i), y(i)); % 异或操作
end
tt = 0:7*N-1; % 生成时间序列
stem(tt, s); axis([0, 350, 0, 1]); % 设置坐标轴范围
title('扩频后的待发送序列码');
阅读全文