1.      采用fir1函数设计,fir1函数可以设计低通.带通.高通.带阻等多种类型的具有严格线性相位特性的FIR滤波器.语法形式: b = fir1(n, wn) b = fir1(n, wn, ‘ftype’) b = fir1(n, wn, ‘ftype’, window) b = fir1(n, wn, ‘ftype’, window, ‘noscale’) 参数的意义及作用: b:返回的FIR滤波器单位脉冲响应,脉冲响应为偶对称,长度为n+1: n:滤波器的介数: wn:滤波器的…

CIC滤波器是滑动平均滤波器的非常高效的迭代实现,只需要一个减法和一个加法,而滑动平均需要N-1个加法. cic滤波器相当于一个梳状滤波器y(n)=x(n)-x(n-D),H(z)=1-z-D,和一个积分滤波器y(n)=x(n)+y(n-1),H(z)=(1-z-1)-1的级联,两个级联后y(n)=x(n)-x(n-D)+y(n-1),H(z)=(1-z-D)/(1-z-1),和求和表达式y(n)=x(n)+x(n-1)+x(n-2)...+x(n-D),H(z)=(1-z-D)/(1-z-1)…

一.设计模糊控制器1.1 创建项目文件夹在此路径如图 1.2 打开MATLAB打开MATLAB R2012a切换当前目录为上一步路径,如图 1.3 设计模糊控制器打开模糊控制器设计对话框 根据模糊控制器的输入输出设计模糊控制器,在此以二输入一输出为例. 完成后如图(左)所示,然后对每个输入输出变量设置隶属函数,如图(右). 添加论域数量 设置隶属函数 完成后如图所示 设计模糊规则 保存刚刚设计的模糊控制器,如下图所示 加载模糊控制器到MATLAB中 二.设置控制系统打开simulink仿真器 设…

前言:本程序是我去年实现论文算法时所做.主要功能为标记切割肝脏区域.时间有点久,很多细节已经模糊加上代码做了很多注释,因此在博客中不再详述. NOTE: 程序分几大段功能模块,仔细阅读,对解决医学图像还是有一定的借鉴意义 想借鉴本文的一定要仔细阅读代码和注释,中间有人机交互部分,空跑会抛异常 .dcm数据,我放到了我的百度云盘,有兴趣的可以下载,实测一下代码.dcm数据连接 clc,clear img_1=dicomread('10011.dcm');%读取dcm文件 (所谓的灰度值) meta…

simulink 使用的动态系统仿真,仿真需要求状态方程和输出方程,关键是求状态方程,而状态方程的求解有多种算法,可变步长和定步长,所以仿真时对求解器的选择和步长的设置就比较重要. 所谓基於数据流的仿真方式,就是指在整个仿真时间段上,根据信号流动的先后顺序逻辑,从信号输入端开始,对一个子系统进行仿真得出在整个仿真时间段上的状态数据和输出,然后基於这些数据再对下一个子系统进行仿真,直到整个系统的信号输出端.在基於数据流的仿真方式下,只有当上一个阶段的计算完成后,才开始进入下一阶段的计算,因此在求解…

第一份工作时做啸叫抑制的仿真,调大0.3可以有大的啸叫产生,下图的SIMULINK仿真模型 实现移相有多种方法: 1.iir实现 2.FFT实现 3.使用FIR实现 所有信号均可以由正弦信号叠加而成.…

1.只需要在运行仿真后,在命令窗口内输入: ,'ShowHiddenHandle','on'); set(gcf,'menubar','figure'); scope最上方会出现一个菜单栏,选择Tools->Edit Plot,即可修改图像所有属性. 2.双击Scope->Parameters->Data History->勾选Save data to workspace 在命令窗口输入: >> tt = ScopeData.time; >> xx = Sc…

1.spatialgabor.m描述gabor函数 % SPATIALGABOR - applies single oriented gabor filter to an image%% Usage:%  [Eim, Oim, Aim] =  spatialgabor(im, wavelength, angle, kx, ky, showfilter)%% Arguments:%         im         - Image to be processed.%         wavel…

clc; Rs=10*10^6; %码速率 es=0.707; %阻尼系数 fs=8*10^6; %采样频率.系统时钟频率 Bit_Loop_out=27;%环路滤波器输入输出位宽 N=31; %NCO频率字位宽 Tdds=8; %NCO频率字更新所需的系统时钟周期个数 K=2*pi*Tdds*2^(Bit_Loop_out-2)/2^N %环路总增益 K=K*90 Wn=73.867*10^3; %环路固有谐振频率(rad/s) c1=((8*es*Wn/fs)+4*Wn*Wn/fs/fs)/…

注:本设计的参数为:D=2,R=5,N=3:时钟频率为50mhz,输入信号为有符号8位,根据公式bmax=bin+N*log(2,R*D):可以得到bmax=18: 1,cic抽取滤波器原理 网上资料一大堆,不说了.重点在于传递函数,以及各个部分的结构. 2,simulink仿真 模型图 频谱仪显示结果 3,cic滤波器verilog 代码 module cic_dec(clk,rst_n,datain,dataout);input clk,rst_n;input [7:0] datain;ou…