05 2019 档案

《DSP using MATLAB》Problem 8.3
摘要：代码： 运行结果： 系统函数部分分式展开，查表求逆z变换就可得到h(n) 零极点的模和幅角 将脉冲序列当成输入得到h_chk(n)，系统函数求逆z变换得到h(n)， 二者幅度谱、相位谱、群延迟对比如下，可见，幅度谱一样，相位谱和群延迟有所不同。 阅读全文
posted @ 2019-05-29 20:09 跑啊跑 阅读(189) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 8.2
摘要：代码： 运行结果： 这里的G为增益系数，使得幅度谱在共振频率处最大，等于1 系统函数部分分式展开 系统函数的零极点图 h_chk是由将脉冲序列当成系统输入而得到的，h是由系统函数部分分时展开后查表求逆z变换得到的， 二者幅度谱一致，但是相位谱和群延迟稍有不同。 r=0.9和0.99的结果这里就不放了 阅读全文
posted @ 2019-05-28 20:38 跑啊跑 阅读(189) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 8.1
摘要：代码： 运行结果： 系统函数部分分式展开， 零极点的模和幅角： 用脉冲序列当输入得到脉冲响应序列h_chk(n)，系统函数H(z)取逆z变换得h(n)，二者如下图 h_chk(n)的幅度谱、相位谱、群延迟 h(n)的幅度谱、相位谱、群延迟 r=0.9、0.99的图这里就不放了。 阅读全文
posted @ 2019-05-26 20:10 跑啊跑 阅读(242) 评论(0) 推荐(0)
linux查看系统内容
摘要： 阅读全文
posted @ 2019-05-25 19:56 跑啊跑 阅读(200) 评论(0) 推荐(0)
redhat4.4下安装GMT4.5.11
摘要：GMT是地学界常用的开源软件，不仅是因为其开源的特性，还有着独特的魅力。 所需要的软件如下 安装步骤： 1. Put the soft packages in one folder, i.e. /home/gmt/soft 2. tar zxvf gdal-1.10.1.tar.gz cd gdal 阅读全文
posted @ 2019-05-24 22:01 跑啊跑 阅读(530) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.38
摘要：代码： 运行结果： 运用Parks-McClellen算法，滤波器长度M=27，比较合适 脉冲响应、幅度谱和误差函数 振幅谱 下面是P7.9的结果，可看出，P-M法得到的长度（M=27)比窗函数（M=43）得到的小得多。 第7章终于弄完了，内心小欢喜一下，放一张田地的照片放松放松，明天开始第8章！ 阅读全文
posted @ 2019-05-23 20:33 跑啊跑 阅读(167) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.37
摘要：代码： 运行结果： 先用firpmord求出N、f、m、weights，这些参数作为firpm函数参数，再求出脉冲响应序列 求出脉冲响应序列的幅度谱、误差响应函数 振幅谱，按照题目要求画了网格线 依前面经验，Parks-McClellan方法要比窗函数法、频率采样法设计的滤波器长度小，节约硬件资源， 阅读全文
posted @ 2019-05-23 07:31 跑啊跑 阅读(250) 评论(0) 推荐(0)
fftw3.3.3在redhat4.4下安装
摘要：FFTW ( the Faster Fourier Transform in the West) 是一个快速计算离散傅里叶变换的标准C语言程序集，其由MIT的M.Frigo 和S. Johnson 开发。可计算一维或多维实和复数据以及任意规模的DFT。 我们需要的是fftw单精度版本 解压fftw软 阅读全文
posted @ 2019-05-22 06:55 跑啊跑 阅读(858) 评论(0) 推荐(0)
gdal在redhat4.4下安装
摘要：GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。 OGR是GDAL项目的一个分支，功能与GDAL类似，只不过它提供对矢量数 阅读全文
posted @ 2019-05-21 20:39 跑啊跑 阅读(306) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.36
摘要：代码： 运行结果： 窗函数法，使用了矩形窗、三角窗、Hann窗、Hamming窗、Blackman窗、Kaiser窗， 1、Rectangular窗 2、Bartlett三角窗 3、Hann、Hamming窗、Blackman窗的图这里不放了，直接放Kaiser窗的结果 4、频率采样方法 5、PM法 阅读全文
posted @ 2019-05-21 17:01 跑啊跑 阅读(303) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.35
摘要：代码： 运行结果： 30阶时，稍有欠缺 增大滤波器长度，M=35时， 幅度谱、误差函数如下 振幅谱如下 阅读全文
posted @ 2019-05-19 17:03 跑啊跑 阅读(157) 评论(0) 推荐(0)
cmake 2.8.12在redhat 4.4下安装
摘要：以前安过，忘了，今天记笔记这里 阅读全文
posted @ 2019-05-19 12:49 跑啊跑 阅读(195) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.34
摘要：代码： 运行结果： 用P-M法设计出来滤波器长度M=23 幅度谱如下 振幅谱如下： 第2个输入，及对应的输出 第2个输入输出的谱（DTFT）对比： 从图中看，输入0.5π分量幅度有所增大，几乎达到原来的2倍。 阅读全文
posted @ 2019-05-18 20:25 跑啊跑 阅读(349) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.33
摘要：代码： 运行结果： 设计一个50阶（即长度M=51）的线性相位FIR，通带宽度不超过0.02π，阻带衰减达到30dB， 最后要把输入中的高斯噪声过滤掉。 As=33dB，满足设计要求。 用P-M方法设计的脉冲响应幅度谱， 振幅谱 输入信号 滤波前后，输入数出对比 输入输出的谱： 右下图可见，随即噪声 阅读全文
posted @ 2019-05-16 20:57 跑啊跑 阅读(235) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.32
摘要：代码： 运行结果： 阅读全文
posted @ 2019-05-14 22:06 跑啊跑 阅读(219) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.31
摘要：参照Example7.27，因为0.1π=2πf1 f1=0.05，0.9π=2πf2 f2=0.45 所以0.1π≤ω≤0.9π，0.05≤|H|≤0.45 代码： 运行结果： 根据线性相位FIR性质，differentiator为第3类线性相位FIR，下图为脉冲响应、幅度谱和振幅谱。 脉冲响应和 阅读全文
posted @ 2019-05-11 16:43 跑啊跑 阅读(193) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.30
摘要：代码： 运行结果： 滤波器长度M=47，阻带衰减满足设计指标。 幅度谱和相位谱 振幅谱，把阻带、通带放大，数数极值点的个数。 下图，9个极值点 下图，8个极值点 下图，9个极值点 总共有9+8+9=26个极值点，M=47，L=(M-1)/2=23，0到π上，最多L+3=26个极值点。 输入输出序列 阅读全文
posted @ 2019-05-08 23:04 跑啊跑 阅读(193) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.29
摘要：代码： 运行结果： 阶数N=22时，阻带衰减43dB，不满足要求。 所以增大N，当N=24时，As=52dB，滤波器长度M=25 阅读全文
posted @ 2019-05-07 20:38 跑啊跑 阅读(190) 评论(0) 推荐(0)
《DSP using MATLAB》Problem 7.28
摘要：又是一年五一节，朋友圈都是晒名山大川的，晒脑袋的，我这没钱的待在家里上网转转吧 频率采样法设计带通滤波器，过渡带中有一个样点 代码： 运行结果： 频率采样法，得到的脉冲响应、振幅响应和幅度响应： 幅度响应、相位响应和群延迟 采用fir2函数求脉冲响应序列 matlab自带的fir2函数比频率采样法得 阅读全文
posted @ 2019-05-04 13:48 跑啊跑 阅读(418) 评论(0) 推荐(1)