沧海一粟

摘要： 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 阅读全文

摘要： 代码： 运行结果： 运用Parks-McClellen算法，滤波器长度M=27，比较合适 脉冲响应、幅度谱和误差函数 振幅谱 下面是P7.9的结果，可看出，P-M法得到的长度（M=27)比窗函数（M=43）得到的小得多。 第7章终于弄完了，内心小欢喜一下，放一张田地的照片放松放松，明天开始第8章！ 阅读全文

摘要： 代码： 运行结果： 先用firpmord求出N、f、m、weights，这些参数作为firpm函数参数，再求出脉冲响应序列 求出脉冲响应序列的幅度谱、误差响应函数 振幅谱，按照题目要求画了网格线 依前面经验，Parks-McClellan方法要比窗函数法、频率采样法设计的滤波器长度小，节约硬件资源， 阅读全文

摘要： FFTW ( the Faster Fourier Transform in the West) 是一个快速计算离散傅里叶变换的标准C语言程序集，其由MIT的M.Frigo 和S. Johnson 开发。可计算一维或多维实和复数据以及任意规模的DFT。 我们需要的是fftw单精度版本 解压fftw软 阅读全文

摘要： GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。 OGR是GDAL项目的一个分支，功能与GDAL类似，只不过它提供对矢量数 阅读全文

摘要： 代码： 运行结果： 窗函数法，使用了矩形窗、三角窗、Hann窗、Hamming窗、Blackman窗、Kaiser窗， 1、Rectangular窗 2、Bartlett三角窗 3、Hann、Hamming窗、Blackman窗的图这里不放了，直接放Kaiser窗的结果 4、频率采样方法 5、PM法 阅读全文

摘要： 代码： 运行结果： 30阶时，稍有欠缺 增大滤波器长度，M=35时， 幅度谱、误差函数如下 振幅谱如下 阅读全文

摘要： 以前安过，忘了，今天记笔记这里 阅读全文

摘要： 代码： 运行结果： 用P-M法设计出来滤波器长度M=23 幅度谱如下 振幅谱如下： 第2个输入，及对应的输出 第2个输入输出的谱（DTFT）对比： 从图中看，输入0.5π分量幅度有所增大，几乎达到原来的2倍。 阅读全文

摘要： 代码： 运行结果： 设计一个50阶（即长度M=51）的线性相位FIR，通带宽度不超过0.02π，阻带衰减达到30dB， 最后要把输入中的高斯噪声过滤掉。 As=33dB，满足设计要求。 用P-M方法设计的脉冲响应幅度谱， 振幅谱 输入信号 滤波前后，输入数出对比 输入输出的谱： 右下图可见，随即噪声 阅读全文