把值付给四个变量
wp=2pi5000;ws=2pi12000;Rp=2;As=30;
1 buttord
(1)[N,wc]=buttord(wp,ws,αp,αs)
调用参数wp,ws分别为数字滤波器的通带、阻带截止频率的归一化值,要求:0≤wp≤1,0≤ws≤1。1表示数字频率pi。
αp,αs分别为通带最大衰减和组带最小衰减(dB)。
当ws≤wp时,为高通滤波器;
当wp和ws为二元矢量时,为带通或带阻滤波器,这时wc也是二元向量。
N,wc作为butter函数的调用参数。
(2)[N,Ωc]=buttord(Ωp,Ωs,αp,αs,‘s’)
用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。
Ωp,Ωs,Ωc均为实际模拟角频率。
说明:buttord函数使用阻带指标计算3dB截止频率,这样阻带会刚好满足要求,而通带会有富余。
2buttap(N)
[z0,p0,k0]=buttap(N)
用于计算N阶巴特沃斯归一化(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子。
说明:如果要从零、极点模型得到系统函数的分子、分母多项式系数向量ba、aa,可调用
[B,A]=zp2tf(z0,p0,k0)
3butter
(1)[b,a]=butter(N,wc,‘ftype’)
计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。
调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值(关于pi归一化),一般是调用buttord(1)格式计算N和wc。
系数b、a是按照z-1的升幂排列。
(2)[B,A]=butter(N,Ωc,‘ftype’,‘s’)
计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量ba、aa。
调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率),可调用buttord(2)格式计算N和Ωc。
系数B、A按s的正降幂排列。
tfype为滤波器的类型:
◇ftype=high时,高通;Ωc只有1个值。
◇ftype=stop时,带阻阻;此时Ωc=[Ωcl,Ωcu],分别为带阻滤波器的通带3dB下截止频率和上截止频率。
◇ ftype缺省时:
若Ωc只有1个值,则默认为低通;
若Ωc有2个值,则默认为带通;其通带频率区间Ωcl < Ω < Ωcu。
注意:所设计的带通和带阻滤波器系统函数是2N阶。因为带通滤波器相当于N阶低通滤波器与N阶高通滤波器级联。
[N,wc]=buttord(wp,ws,Rp,As,'s');
[B,A]=butter(N,wc,'s');
k=0:511;fk=0:14000/512:14000;wk=2pifk;
Hk=freqs(B,A,wk);
画图,两行两列
subplot(2,2,1);
plot(fk/1000,20log10(abs(Hk)));grid on
坐标轴意义
xlabel('频率(kHz)');ylabel('幅度(dB)')
坐标的取值范围
axis([0,14,-40,5])
imread:读图像存入矩阵中
rgb2gray:RGB彩图转灰度图
fun是一个函数,作用是将返回一个与输入矩阵x一样尺寸的矩阵,其中每个元素都是输入矩阵的x标准差
blkproc:对图像分块处理,本程序中表示分成32×32的块,每个块应用fun函数
通俗点讲,最后效果是把输入彩图转灰度图并马赛克化。
如果n是定值的话,比如说100
程序如下:
xn=[];
for x=1:1:100
xn(x)=x+1;
end
xn
然后xn就可以显示了。
或者你如果要n是可以输入的话你就新建一个shuliem文件如下
function xn=shulie(x)
for x=1:1:100
xn(x)=x+1;
end
然后主程序调用就是在command window 输入shulie(100)就可以 了。
不懂欢迎追问~有帮助请采纳~
直接调用conv函数就可以实现,具体程序如下:
>> x=[1,2,3,4];
>> y=[1,1,1,1,1];
>> z=conv(x,y)
z =
1 3 6 10 10 9 7 4
plot(z)
希望我的解答能让你满意!
以上就是关于matlab一段简单的程序,麻烦注释每一句的意思.全部的内容,包括:matlab一段简单的程序,麻烦注释每一句的意思.、matlab简单程序、简单小程序求解:::matlab中计算Xn=n+1;对于不同的n的值,将X1至Xn的值都计算出来,求大神等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)