查看: 8217|回复: 3|关注: 0
Matlab如何让显示的有效位数由4为到8位
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
&& x=[6.173,8.231,9.62,10.288,12.347,14.919];
&& y=[-00,-45];
&& t=polyfit(x,y,4)
&&1.0e+005 *
& &-0.0013& & 0.0217& &-0.1787& &-1.6351& & 7.3405
问题：我想让t得到的数据有效位数显示的多一些，比如由目前的4位变为8位有效数字，！（注意是显示，因为我想把系数记录下来用在另一个语言中），多谢各位大侠！
<h1 style="color:#1 麦片财富积分
关注者： 42
<h1 style="color:# 麦片财富积分
就可以了。
关注者： 186
1、可以通过format来设置数据精度；
2、可以通过vpa自定义精度或者digits设置精度；
站长推荐 /2
Powered by5543人阅读
mat_lab（50）
x=[20,22,24,26,28,30,32,34,36,38,40,42,44];y=[62.9,68.8,71.2,82.5,84.1,88.6,88.4,88.4,88.0,88.0,88.0,88.0,88.0];plot(x,y,'-r*');xlabel('Number of Gabor features');ylabel('classification recognition rate (%)');legend('pain vs.non-pain','location','northwest');xlim([20, 45]);ylim([60, 90]);
set(gca,'XTick',[20:2:45]) %改变x轴坐标间隔显示 这里间隔为2axis([10,50,50,100])& % axis([xmin,xmax,ymin,ymax])用这个语句调整坐标轴数字显示
参考知识库
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：1730863次
积分：20681
积分：20681
排名：第265名
原创：365篇
转载：251篇
评论：292条
(1)(1)(2)(1)(7)(4)(3)(14)(3)(8)(2)(2)(11)(25)(29)(1)(10)(26)(153)(85)(15)(15)(13)(1)(1)(11)(9)(8)(4)(59)(48)(8)(14)(24)(3)(2)(2)摘要《数字信号处理》课程是一门理论性和实践性都很强，它具备高等代数、数值分析、概率统计、随机过程等计算学科的知识;要求我们学生掌握扎实的基础知识和理论基础。又是跟其他学科密切相关，即与通信理论、计算机、微电子技术不可分，又是人工智能、模式识别、神经网络等新兴学科的理论基础之一。本次数字滤波器设计方法是基于MATLAB的数字滤波器的设计。此次设计的主要内容为：IIR数字滤波器和FIR数字滤波器的设计关键词：IIR、FIR、低通、高通、带阻、带通Abstract"DigitalSignalProcessing"isatheoreticalandpracticalnaturearestrong,andithasadvancedalgebraandnumericalanalysis,probabilityandstatistics,randomprocesssuchascalculationorequiresstudentstoacquirebasicknowledgeandasolidtheoreticalbasis.Iscloselyrelatedwithothersubjects,namely,andcommunicationtheory,computers,microelectronicscannotbeseparated,butalsoinartificialintelligence,patternrecognition,neuralnetworktheoryoneoftheemergingdiscipline.ThedigitalfilterdesignmethodisbasedonMATLABfordigitalfilterdesign.Themainelementsofdesign:IIRandFIRdigitalfilterdesignofdigitalfilterKeyWords:IIR,FIR,lowpass,highpass,bandstop,bandpass目录一、前言3二、课程设计的目的3三、数字信号处理课程设计说明及要求3四、滤波器的设计原理44.1数字滤波器简介44.2IIR滤波器的设计原理44.3FIR滤波器的设计原理54.4FIR滤波器的窗函数设计法6五、设计内容65.1设计题目：65.2设计程序代码及结果：7六、结束语15七、参考文献16一、前言数字信号处理（DigitalSignalProcessing，简称DSP）是一门涉及许多学科而又广泛应用于许多领域的新兴学科。随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。二、课程设计的目的1)三、数字信号处理课程设计说明及要求所需硬件：PC机四、滤波器的设计原理4.1数字滤波器简介数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应，让它满足一定的要求，从而对通过该系统的信号的某些特定的频率成分进行过滤，这就是滤波器的基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统是一个离散系统，则滤波器称为数字滤波器。信号通过线性系统后，其输出就是输入信号和系统冲激响应的卷积。除了外，的波形将不同于输入波形。从频域分析来看，信号通过线性系统后，输出信号的频谱将是输入信号的频谱与系统传递函数的乘积。除非为常数，否则输出信号的频谱将不同于输入信号的频谱，某些频率成分较大的模，因此，中这些频率成分将得到加强，而另外一些频率成分的模很小甚至为零，中这部分频率分量将被削弱或消失。因此，系统的作用相当于对输入信号的频谱进行加权。4.2IIR滤波器的设计原理IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计，通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。IIR数字滤波器的设计步骤：（1）按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标；（2）根据模拟滤波器技术指标设计为响应的模拟低通滤波器；（3）很据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器；（4）如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。4.3FIR滤波器的设计原理FIR滤波器通常采用窗函数方法来设计。窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后街区（加窗）它的脉冲响应得到线性相位和因果FIR滤波器。我们用Hd（e^jw）表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽wc&pi的低通滤波器由下式给定：为了从hd(n)得到一个FIR滤波器，必须同时在两边截取hd（n）。而要得到一个因果的线性相位滤波器，它的h(n)长度为N，必须有：这种操作叫做加窗，h(n)可以看做是hd(n)与窗函数w（n）的乘积：h(n)=hd(n)w(n)其中根据w（n）的不同定义，可以得到不同的窗结构。在频域中，因果FIR滤波器响应H(e^jw)由Hd(e^jw)和窗响应W(e^jw)的周期卷积得到，即常用的窗函数有矩形窗、巴特利特（BARTLETT）窗、汉宁（HANNING）窗、海明（HAMMING）窗、布莱克曼（BLACKMAN）窗、凯泽（KAISER）窗等。4.4FIR滤波器的窗函数设计法FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列，主要设计步骤为：(1)通过傅里叶逆变换获得理想滤波器的单位脉冲响应hd(n)。(2)由性能指标确定窗函数W(n)和窗口长度N。(3)求得实际滤波器的单位脉冲响应h(n)，h(n)即为所设计FIR滤波器系数向量b(n)。五、设计内容5.1设计题目：1-1.试用MATLAB设计一巴特沃斯低通数字滤波器，要求通带截至频率Wp=30HZ，主带截至频率为Ws=35HZ，通带衰减不大于0.5DB，主带衰减不小于40DB，抽样频Fs=100HZ。1-2．基于Butterworth模拟滤波器原型，使用双线性状换设计数字滤波器：各参数值为：通带截止频率Omega=0.2*pi,阻带截止频率Omega=0.3*pi,通带波动值Rp=1dB,阻带波动值Rs=15dB,设Fs=20000Hz。1-3设计一巴特沃斯高通数字滤波器，要求通带截止频率0.6*pi，通带衰减不大于1dB，阻带衰减15DB，抽样T=1。1-4.设计一巴特沃斯带阻数字滤波器，要求通带上下截至频率为0.8*PI、0.2*PI，通带衰减不大于1DB，阻带上下截至频率0.7*PI、0.4*PI阻带衰减不小于30DB，2-1.用窗函数法设计一个线性相位FIR低通滤波器，并满足性能指标:通带边界频率Wp=0.5*pi,阻带边界频率Ws=0.66*pi,阻带衰减不小于40dB,通带波纹不大于3dB。选择汉宁窗。2-4．用海明窗设计一个FIR滤波器，其中Wp=0.2*pi，Ws=0.3*pi,通带衰减不大于0.25dB,阻带衰减不小于50dB。5.2设计程序代码及结果：1-1一.试用MATLAB设计一巴特沃斯低通数字滤波器，要求通带截至频率Wp=30HZ，阻带截至频率为Ws=35HZ，通带衰减不大于0.5DB，阻带衰减不小于40DB，抽样频Fs=100HZ。代码为：fp=30;fs=35;Fs=100;wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;wp=tan(wp/2);ws=tan(ws/2);%通带最大衰减为0.5dB，阻带最小衰减为40dB[N,wn]=buttord(wp,ws,0.5,40,'s');%模拟低通滤波器极零点[z,p,k]=buttap(N);%由极零点获得转移函数参数[b,a]=zp2tf(z,p,k);%由原型滤波器获得实际低通滤波器[B,A]=lp2lp(b,a,wp);[bz,az]=bilinear(B,A,.5);[h,w]=freqz(bz,az,256,Fs);figureplot(w,abs(h))gridon图1巴特沃斯数字低通滤波器1-2基于Butterworth模拟滤波器原型，使用双线性状换设计数字滤波器：各参数值为：通带截止频率Omega=0.2*pi,阻带截止频率Omega=0.3*pi,通带波动值Rp=1dB,阻带波动值Rs=15dB,设Fs=4000Hz。代码：wp=0.2*ws=0.3*Fs=4000;T=1/Fs;OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);rp=1;rs=15;as=15;ripple=10^(-rp/20);attn=10^(-rs/20);[n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');[z,p,k]=Buttap(n);[b,a]=zp2tf(z,p,k);[bt,at]=lp2lp(b,a,wn);[b,a]=bilinear(bt,at,Fs);[db,mag,pha,grd,w]=freqz_m(b,a);%%下面绘出各条曲线subplot(2,2,1);plot(w/pi,mag);title('MagnitudeFrequency幅频特性');xlabel('w(/pi)');ylabel('|H(jw)|');axis([0,1,0,1.1]);set(gca,'XTickMode','manual','XTick',[00.20.31]);set(gca,'YTickMode','manual','YTick',[0attnripple1]);gridsubplot(2,2,2);plot(w/pi,db);title('MagnitudeFrequency幅频特性(db)');xlabel('w(/pi)');ylabel('dB');axis([0,1,-30,5]);set(gca,'XTickMode','manual','XTick',[00.20.31]);set(gca,'YTickMode','manual','YTick',[-60-as-rp0]);gridsubplot(2,2,3);plot(w/pi,pha/pi);title('PhaseFrequency相频特性');xlabel('w(/pi)');ylabel('pha(/pi)');axis([0,1,-1,1]);subplot(2,2,4);plot(w/pi,grd);title('GroupDelay群延时');xlabel('w(/pi)');ylabel('Sample');axis([0,1,0,15]);set(gca,'XTickMode','manual','XTick',[00.20.31]);grid运行结果：图2巴特沃思数字低通滤波器幅频-相频特性1-3设计一巴特沃斯高通数字滤波器，要求通带截止频率0.6*pi，通带衰减不大于1dB，阻带衰减15DB，抽样T=1。Wp=0.6*Ws=0.4*Ap=1;As=15;[N,wn]=buttord(Wp/pi,Ws/pi,Ap,As);%计算巴特沃斯滤波器阶次和截止频率%频率变换法设计巴特沃斯高通滤波器[db,mag,pha,grd,w]=freqz_m(b,a);%数字滤波器响应plot(w,mag);title('数字滤波器幅频响应|H(ej\Omega)|')图3巴特沃斯数字高通滤波器2-1用窗函数法设计一个线性相位FIR低通滤波器，并满足性能指标:通带边界频率Wp=0.5*pi,阻带边界频率Ws=0.66*pi,阻带衰减不小于40dB,通带波纹不大于3dB。选择汉宁窗。代码：wp=0.5*ws=0.66*wdelta=ws-N=ceil(8*pi/wdelta)ifrem(N,2)==0N=N+1;end);运行结果：给分就给你个全的!图6低通FIR滤波器六、结束语本次数字滤波器设计方法是基于MATLAB的数字滤波器的设计，是用学过的数字信号理论为依据，用MATLAB代码来实现。课程设计过程中，通过IIR数字滤波器和FIR数字滤波器的设计实例，说明如何利用MATLAB来完成数字滤波器的设计。窗函数法中相位响应有严格的线性,不存在稳定性问题,设计简单。双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，但会产生频率混碟现象，使数字滤波器的频响偏移模拟滤波器的频响。由滤波器的频谱图和滤波前后的语音信号的频谱图对比可知本设计选用双线性变换法设计的IIR滤波器比较好。在同样的技术指标的要求下，IIR滤波器所要求的阶数N也比较小，实现起来比较容易。通过综合运用数字信号处理的理论知识进行滤波器设计，通过理论推导得出相应结论，再利用MATLAB作为编程工具进行计算机实现，从而加深了对所学知识的理解，建立概念。对以前在课本上所学的东西有了更深入的理解和掌握。最后，无论做什么课程设计，都需要要有一定的理论知识作为基础，同时通过这次课程设计，我对于以前所学的数字信号处理知识有了更深的理解。七、参考文献1.程佩青《数字信号处理教程》北京清华大学出版社2007年2月.2.赵知劲、刘顺兰《数字信号处理实验》.浙江大学出版社.3.S.K.Mitra．DigitalSignalProcessing:AComputer-BasedApproach.NewYork,NewYork:McGraw-Hill,thirded,2006．4.肖伟、刘忠等《MATLAB程序设计与应用》清华大学出版社、北京交通大学出版社.5.胡良剑、孙晓君《MATLAB数学实验》.高等教育出版社.
略懂社热议
去略懂看看
问问青春版，什么都略懂一点
去略懂看看