![](https://img-blog.csdnimg.cn/20190918140037908.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
数字信号处理(DSP)
数字信号处理
ReCclay
CSDN认证博客专家、CSDN万粉博主、CSDN嵌入式领域优质创作者、CSDN-2020博客之星年度总评选Top25。自入站以来,凭借坚持与热爱,以博文的方式分享所学,累计博文数量达1千余篇,受益人次达300w+次,涉及领域包括但不限于前/后端软件开发、嵌入式软件开发、Linux驱动开发、自动驾驶算法研究、视觉SLAM算法研究、ROS机器人开发、FPGA开发、IC设计与验证等。欢迎您的关注,我们一起进步!
展开
-
数字信号中的上采样和下采样
数字信号处理中的上采样和下采样~原创 2022-01-07 14:53:47 · 13741 阅读 · 0 评论 -
巴特沃兹低通滤波器设计之已知阶数及3dB截止频率,求增益响应曲线
代码如下% 4阶模拟巴特沃斯低通滤波器设计clc;clear;% 15位近似定点数format long;% 确定零点/极点/增益因子[z,p,k] = buttap(4);disp('极点位于');disp(p);% 传递函数系数的确定[pz, pp] = zp2tf(z, p, k);% S幂递减的列印系数disp('分子多项式系数'); disp(pz);disp...原创 2020-01-17 20:41:52 · 5260 阅读 · 0 评论 -
Matlab中的DSP应用中的bilinear函数
利用双线性变换法,把下列S域的模拟滤波器转换为数字滤波器clc;clear;bs=[1 1];as=[1 5 6];T=0.1;Fs=1/T;[bz,az]=bilinear(bs,as,Fs)过程及结果bz = 0.0415 0.0040 -0.0375az = 1.0000 -1.5573 0.6047...原创 2020-01-17 20:39:39 · 17018 阅读 · 0 评论 -
Matlab中的DSP应用中的impinvar函数
利用脉冲响应不变法,把下列S域的模拟滤波器转换为数字滤波器clc;clear;b=[1,1];a=[1 5 6];T=0.1;[bz,az]=impinvar(b,a,1/T)输出bz = 0.1000 -0.0897az = 1.0000 -1.5595 0.6065...原创 2020-01-17 20:39:55 · 9970 阅读 · 0 评论 -
巴特沃斯滤波器设计中的常用函数
1、[z p k] = buttap(N)计算3dB截止频率为1,阶数为N的归一化巴特沃斯低通滤波器,输出为零点、极点、系统增益。要想得到对应的传输函数,可以[b, a] = zp2tf(z, p, k),b和a分别对应传输函数的分子和分母系数。2、[num den] = butter(N, Wn, 's')由3dB截止频率和阶数N,设计一个低通滤波器。若Wn是一个含有两个元素的向量[W...原创 2020-01-17 20:40:05 · 8428 阅读 · 0 评论 -
DSP实验报告(一)之离散时间信号和系统的时域分析
一、产生如下因果线性时不变系统的冲激响应的前45个样本并绘图:clf;N = 45;num = [0.9 -0.45 0.35 0.002];den = [1 0.71 -0.46 -0.62];y = impz(num,den,N);% Plot the impulse responsestem(y);xlabel('Time index n'); ylabel('Amplit...原创 2020-01-17 20:40:15 · 2644 阅读 · 0 评论 -
DSP实验报告(二)之离散时间信号和线性时不变离散时间系统的频域分析
长度为6的滑动平均滤波器的增益响应clf;M = 6;%设定滤波器的长度num = ones(1,M)/M;[g,w] = gain(num,1);%求增益函数plot(w/pi,g);gridaxis([0 1 -50 0.5])xlabel('\omega /\pi');ylabel('Gain in dB');title(['M = ', num2str(M)])gain...原创 2020-01-17 20:40:22 · 894 阅读 · 0 评论 -
DSP实验报告(三)之连续时间信号的数字处理
一、对连续时间信号进行离散抽样,样例程序clf;t = 0:0.0005:1;f = 3;xa = cos(2*pi*f*t);subplot(2,1,1)plot(t,xa);gridxlabel('Time, msec');ylabel('Amplitude');title('Continuous-time signal x_{a}(t)');axis([0 1 -1.2 1...原创 2020-01-17 20:40:28 · 1572 阅读 · 0 评论 -
DSP实验报告(四)之数字滤波器结构
一、生成如下因果无限冲激响应传输函数的级联实现,画出级联实现的结构图。程序clc;clear;format rat;num = input('分子系数向量 = ');den = input('分母系数向量 = ');[z,p,k] = tf2zp(num,den);%由传递函数转换为零极点增益sos = zp2sos(z,p,k)%由零极点增益转换为二次分式过程及结果分子...原创 2020-01-17 20:40:41 · 2449 阅读 · 0 评论 -
DSP实验报告(五)之离散时间信号和系统的时域分析
一、用双线性变换法及Matlab编程设计一个IIR数字带通滤波器,其中通带波纹为1dB,阻带最小衰减为40dB,通带截止频率分别为500Hz和700Hz,阻带截止频率分别为400Hz和800Hz,抽样频率为2000Hz。画出该滤波器的幅频特性图。No.1巴特沃斯带通滤波器答案版本的…clc;clear;Rp = 1; % 单位为dB的最大通带衰减Rs = 40; % 单位为dB的最小...原创 2020-01-17 20:39:37 · 2845 阅读 · 1 评论 -
Matlab中DSP的系统函数的表示法及其转换
系统函数表示一、传递函数法a=[1 a(2) a(3)… a(N)]b=[b(1) b(2)… b(M)]二、零极点增益法零点向量 Z=[z1 z2 zM-1];极点向量 P=[z1,z2,…,zN-1]k为系统增益。三、部分分式法分母极点向量 p=[p(1) p(2) … p(n)]分子留数向量r=[r(1) r(2) … r(n)]余数多项式系数向量 k=[k(1...原创 2020-01-17 20:40:49 · 3486 阅读 · 0 评论 -
Matlab中的DSP应用中的freqs函数
注意这是个S域分析…返回量H则包含了离散系统频响在0——pi范围内N个频率等分点的值,向量w则包含范围内N个频率等分点。调用中若N默认,默认值为512。clc;clear;a = [1 0.4 1];b = [0.2 0.3 1];w = logspace(-1, 1);%生成从10的a次方到10的b次方之间按对数等分的n个元素的行向量。n如果省略,则默认值为50。自动挑选200个...原创 2020-01-17 20:41:03 · 4654 阅读 · 0 评论 -
Matlab中的DSP应用中的tf2zp函数及zplane(z,p)
[z,p,K]=tf2zp(num,den)求得有理分式形式的系统转移函数的零、极点zplane(z,p)绘出零、极点分布图原创 2020-01-17 20:41:14 · 10207 阅读 · 0 评论 -
Matlab对图像的傅里叶变换初探
原图像:clc;clear;I=imread('1.jpg'); % 读取图像I=rgb2gray(I); % RGB图像转换为灰度图 - (可试试去掉然后做对比)I=im2double(I); % 图像转换为double,此时图像的数值范围由原来的【0,255】转化为【0,1】 -(可试试去掉然后做对比)F=fft2(I); ...原创 2020-01-17 20:41:57 · 12941 阅读 · 3 评论 -
获取微信中的音频方法
针对我的chrome浏览器step1、把微信对应的文章用PC浏览器打开step2、F12 -> 进入审查元素界面step3、依次按照下图操作。注意如果name中没有东西,可以先在网页上播放一下音频。step4、复制对应的链接到新的标签页,然后另存为即可。...原创 2020-01-17 20:41:32 · 2567 阅读 · 0 评论 -
MATLAB中的信号短时傅里叶变换获得信号的时频联合分析
短时傅里叶变换 -> 时频联合分析双音多频DTMF对于一个声音,如果能够检测出其中高频和低频的频率,就可以指出对应的按键编码。现在可以通过专用的数字处理芯片直接获得双音频解码,也可以通过计算软件进行解码,后者就是这篇博文将要演示的。...原创 2020-01-17 20:41:23 · 1855 阅读 · 0 评论 -
巴特沃兹低通滤波器设计之阶数确定
文章目录题目代码输出题目已知一个具有最大平坦低通特性的传输函数Ha(s)H_{a}(s)Ha(s)的1dB截止频率为1kHz,它在5kHz处有一个40dB的最小衰减,试确定它的最低阶数。代码clcclearWp = 2*pi*1000; Ws = 2*pi*5000; Rp = 1; Rs = 40;[N, Wn] = buttord(Wp, Ws, Rp, Rs, 's');%...原创 2020-01-17 20:41:43 · 11819 阅读 · 0 评论 -
DSP课程设计中的一点记录
关于单位频率 f 的单位是:Hz模拟角频率 Ω 的单位是:rad/s数字角频率 w 的单位是:rad模拟角频率和数字角频率之间的关系是:w = Ω*T = 2πfT = 2πf/fs关于FFTN点FFT,其中N的含义是把频率分成N等份(0~N-1)对应关系为:0~N-1 ===> -2/fx ~ 2/fs(-π ~ π) ,所以在取坐标的时候要特别注意。...原创 2020-01-17 20:39:28 · 701 阅读 · 0 评论