汉宁窗 在这里插入图片描述 123456789101112clearclose allclcN=45;wn=[0.3,0.5];%截止频率window=hanning(N+1); %计算长度为N的汉宁窗b=fir1(N,wn,window); %N是阶数,wn是截止频率[H,w]=freqz(b); subplot(2,1,1);plot(w/pi,20*log10(abs(H)));...

符号含义 # 例题 脉冲响应不变法 123456789101112131415161718clearclose allclcfs=1000;%采样频率fc=200;%通带截止频率fr=300;%阻带截止频率T=0.001; %采样周期%%%%%%%脉冲响应不变法wp1=2*pi*fc;%通带截止频率wr1=2*pi*fr;%阻带截止频率[N1,wn1]=buttord(wp1,wr1,1...

1234567891011121314151617181920clearclose allclcfs=1000;fc=200;fr=300;T=0.001; wp1=2*pi*fc;wr1=2*pi*fr; [N1,wn1]=buttord(wp1,wr1,1,25,'s'); [B1,A1]=butter(N1,wn1,'s'); [num1,de...

这题原本画出来就四条谱线，除了频率0的直流分量，就三根了，太少了，为了看着方便，我自己增加了序列 12345678910111213clearclose allclcn=0:8;x=[4 3 2 6 7 8 3 3 5];subplot(1,2,1);stem(n,x)xlabel('n');ylabel('幅度');title('原序列&...

这里我们直接展示补0的核心代码 生成三角波和反三角波 123456789101112N=8;n=0:7;x1=zeros(size(n));x2=zeros(size(n));xa=(n>=0)&(n<=3);x1(xa)=n(xa);x2(xa)=4-n(xa);xa=(n>=4)&(n<=7);x1(xa)=8-n(xa);x2(xa)=n(xa...

xa(n)是高斯序列 xb(n)是衰减正弦序列 线性卷积很简单，利用conv函数就可以 那么圆周卷积要怎么实现呢 12345678910111213141516171819202122clc; n1=0:1:15; %16个点n2=0:1:30; %线性卷积，序列长度变成了16+16-1p=8;q=2;a=0.1;f=0.0625;xa=exp(-((n1-p).^2)/q);xb=exp...

在这里插入图片描述 12345678910111213141516171819clearclose allclcfs=8000;%采样频率n=0:99;%采样点100个y=sin(0.00625*2*pi*n)+sin(0.0625*pi*2*n)+sin(0.125*2*pi*n);%采样后的信号 0.00625=50/8000subplot(2,1,1);stem(n,y)titl...

原函数 假设我们有这么一个函数 x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π) n为0到15的整数 这是他的函数图像，如何将其进行周期化 我在MATLAB中文论坛看到了一个很巧妙的方法 1b = mod(a,m) b = mod(a,m) 返回 a 除以 m 后的余数，其中 a 是被除数，m 是除数。此函数通常称为取模运算，表达式为 b = a-m...

今天做作业的时候发现要对一个信号进行FFT变换，在网上找了半天也没找到个能看懂的（因为我太菜了），后来自己研究了一下，感觉一知半解的 起因是这道作业题 ## 例题-满足奈奎斯特 我画了两个图，一个是原信号经过采样后的离散图，一个就是此信号经过FFT后的频谱图 因为是8kHZ采样，所以信号不会失真，频谱也是正确的 解答如下： 12345678910111213141516171819clea...

因为之前MATLAB没有好好学，而且DSP也没有好好学，做起作业来十分困难，为了好好写作业，我又重新研究了DSP和MATLAB，学了一下信号的抽样和插值 先看一看x(n)原本的图像 代码如下 123456789101112%定义基本参数f=1;%原信号频率fs=16;%抽样信号频率N=52;%采样点数n=0:N-1;%画原序列图像xn=cos(2*pi*n*(f/fs));stem(n,...