Matlab笔记

对x(t)进行插0，即每2个数据点中插入一个值，
clear all;
%%%%%单周期和多周期方波信号频谱
T = 0.01;%采样频率100Hz
t = 0:T:10-T;%信号持续时间10秒
inc = find(t<=5);%占空比为5/10
f = zeros(size(t));f(inc) = 1;%生成单周期方波脉冲信号
subplot(2,2,3);plot(abs(Fx1));title('AR过程功率谱（方法论1）');
h=ifft(H);%AR系统的冲突响应
x=conv(h,Nn);%系统工程输出信号
Rx=conv(x,fliplr(x));%输出信号的自相关
Fx2=fft(Rx);%AR过程功率谱
subplot(2,2,4);plot(abs(Fx2));title('AR过程功率谱（方法2）');
Matlab实践笔记
高某调频信号为u(t)=[1-0.7cos(t)][cos(10 t)+sin(9 t)]
试画出现1-20秒内的信号曲线，画出该曲线的包络线。
clear all;
t = 0:0.01:20;
y = (1-0.7+cos(t)).*(cos(10*pi*t)+sin(9*pi*t));
系数a1、a2为实数，取值如下
=0.1
=-0.8
画出x（n）的功率谱密度
clear all;
%%%%%AR过程功谱%%%%%%%
Nn=randn(1,4096);%产生高斯白噪声
Rn=conv(Nn,fliplr(Nn));%通过卷积
Fn=fft(Rn);%求噪声功谱
subplot(2,2,1); plot(abs(Fn));title('噪声功率谱');%
z1 = interp(x,2);%内插信号，每2个点内插一个值，保持信号包络不变，不是直接插0
Fz1 = fft(z1);%内插信号z1频谱
figure(3);
subplot(3,1,1);stem(x,'.');title('原信号x');
subplot(3,1,2);stem(y1,'.');title('抽样信号y');
figure %用于指定窗口的名称
freqs(b,a);
b = [1 0 5];
a = [3 4 1];
h = freqs(b,a);
freqs(b,a);
例2.4
目的：掌握离散序列频谱的求解和绘制方法。
信号：该例中考查两个正弦序列的叠加后的幅频特性
x(n)=sin(w1*n)+sin(w2*n)
f1 = 20; f2 = 50; %两个连续正弦信号的频率分别为10Hz和50Hz
y = conv(x,h);
subplot(5,1,1);stem(n,u1);title('u1(n)');
subplot(5,,2);stem(n,u2);title('u2(n)');
subplot(5,1,3);stem(n,x);title('x(n)=e(n)-e(n-5)');
subplot(5,1,4);stem(n,h);title('h(n)=(0.5)^n * e(n)');
subplot(3,1,3);stem(z1,'.');title('插值信号z');
figure(4);
subplot(3,1,1);plot(abs(Fx));title('原信号x幅频特性');
subplot(3,1,2);plot(abs(Fy1));title('抽样信号y1幅频特性');
% fs = 0.1; duty = 50;
% f = square(2*pi*fs*t,duty)/2+1/2; %直接生成单周期方波，持续时间t，频率fs，占空比duty
% f = circshift(f,[0,50]); %改变脉冲位置
F = fft(f);%单脉冲方波脉冲频谱
omega = linspace(0,2,length(t));%数字频率0-2PI
%plot(t,y,t,up,'-.',t,down,'-.')
plot(t,y);hold on
plot(up_x，up);hold on;
plot(down_x，down)
例题：假设一个AR过程x(n)由下式产生
x(n)=- x(n-1)- x(n-2)+v(n)
其中v（n）是零均值高斯白噪声，方差
subplot(3,1,1);
stem(n,x); %绘制时域信号
title('时域信号');
subplot(3,1,2);
stem(mFx); %绘制[0,2pi]区间内的频谱
title('[0,2pi]区间内的频谱');
subplot(3,1,3);
stem(n(1:32),mFx(1:32));%绘制[0,pi]区间内的频谱
y(inc) = x(inc);%抽取位置赋原信号x值
Fy = fft(y);%抽样信号频谱
z = kron(x,[0,1]);%内插0，每2个数据点之间内插1个0
% z = upsample(x,2); %内插0，每2个数据点之间内插1个0
Fz = fft(z);
figure(1)
subplot(3,1,1);stem(x,'.');title('原信号x');
a1=0.1; a2=-0.8;
a=[1,a1,a2];
b=1;%AR系统函数的系数
[H,W]=freqz(b,a,length(Fn),'whole');%求AR系统的频率响应，0-2P
subplot(2,2,2);plot(W,abs(H));title('AR系统幅频特性');
Fx1=abs(H).^2.*(Fn.');%直接通过随机过程通过学习LTI系统的性质求输出AE过程的功率谱
up1_inc = find (diff(sign(diff(y)))==-2)+1;
up1 = y(up1_inc);
up1_x = x(up1_inc);
downl_inc = fid (diff(sign(diff(y)))==2)+1;
down1 = y(down1_inc);
down_x = x(down1_inc);
T = 0.01;%采样频率100Hz
t = 0:T:10-T;%信号持续时间10s
f = 0.3;%调频斜率参数
x = sin(pi*f*t.^2);%线性调频信号
Fx = fft(x);% x信号频谱
inc = find(mod(t,5*T)==0);%每2个点抽取一次，确定抽取点的位置
y = 0*t;%其余位置赋0值
fs = 400; %采样频率取为400Hz
w1 = 2*pi*f1/fs; %w1，w2分别对应的数字角频率
w2 = 2*pi*f2/fs;
n = [-31:32]; %数据段范围
x = sin(w1*n)+sin(w2*n);
Fx = fft(x,64); %对信号x作64点FFT
mFx = abs(Fx); %求频谱的模
11月24日
将单脉冲信号周期延拓即得到脉冲串，考虑一个10周期矩形脉冲串，占空比为50%，画出其幅频特性，并与单脉冲的进行比较
将占空比改为5%，画出此时的幅频特性，与前面的有何不同？
数据抽取和插值对频谱的影响
考虑线性调频信号时画出其幅频响应
对x(t)每5个点抽取一次，得到y(t);此时y(t)数据与x(t)一样长，但是有4/5人位置为0，画出y(t)的幅频响应，并与x(t)的比较
subplot(3,1,3);plot(abs(Fz));title('插值信号z1幅频特性');
y1 = downsample(x,5); %第5个点抽样1次，其余位置不赋0值
% y1 = decimate(x,5); %第5个点抽样1次，其余位置不赋0值
Fy1 = fft(y1);%抽样信号y1频谱
subplot(5,1,5);stem(y);title('卷积后的y');
n = -5:34;
u1 = [n>=0];
N = 10;%延拓10个周期
t1 = 0:T:10*N-T;%多周期信号持续时间100s采样频率不变
f1 = kron(ones(1,10),f);%利用矩阵直积进行周期延拓，得到10个周期的就问信号
F1 = fft(f1);%多周期信号频谱
% f1 = square(2*pi*fs*t1,duty)/2+1/2; %直接生成我周期方波，持续t1，
up = up1;
down = dwon1;
clear all;
t = 0:0.01:20;
y = (1-0.7+cos(t)).*(cos(10*pi*t)+sin(9*pi*t));
[up,down] = envelope(t,y);
up1_inc = find (diff(sign(diff(y)))==-2)+1;
title('[0,pi]区间内的频谱');
例2.5
目的：掌握离散序列响应卷积和求解的Matlab实现。
信号：h(n)=(0.5)^n * e(n)