matlab

da=polyval(den,j*wc)。 g=na/da。 g1=abs(g)。 h=20*log10(g1)。 a=10^(h/10)。 wm=wc。 T=1/(wm*(a)^(1/2))。 alphat=a*T。 Gc=tf([T 1],[alphat 1])。 elseif key==3 gama=vars(1)。 wc=vars(2)。 gama1=gama+5。 num={1}。

save s1 A B %将 A、 B 保存到文件 中 clear %清除当前工作空间中的所有变量 load s1 %装入文件 的内容到数组 A,B 中 MATLAB 程序设计与优化函数的应用 _讲稿 11 A %显示 A的内容 A = 1 4 7 2 5 8 B %显示 B 的内容 B = 2 4 6 8 【 例 】 通过 office 建立一个文件 F:\，内容为： 1 2 3 45 6 7

是在低高强度和饱和，这增加了输出图像 J 的对比度值。 14 J = imadjust(I,[low_in。 high_in],[low_out。 high_out]) ， 将图像 I中的亮度值映射到 J中的新值，即将 low_in 至 high_in 之间的值映射到 low_out至 high_out 之间的值。 low_in 以下与 high_in 以上的 值被剪切掉了，也就是说，

图 21a svd变换后运行结果 从左至右依次为：原图像转换为灰度图后的图像，矩阵 U等价的图像，矩阵 S等价的图像 图 21b svd变换后运行结果 从左至右依次为：矩阵 V等价的图像，矩阵 SM等价的图像，矩阵 SN等价的图像 武汉理工大学《 matlab 课程设计》报告 9 图像的正交分解 正交分解理论知识 实数矩阵 A 的 QR 分解 是把 A 分解为 QRA （公式 24） 这里的

击响应长度9，即滤波器阶数为8b=fir1(N1,wc/pi,39。 high39。 ,boxcar(N))。 n=0:1:8。 %调用fir1函数subplot(1,2,1)。 stem(n,b,39。 .39。 ) %输出单位脉冲冲击响应axis([0,8,1])。 grid。 xlabel(39。 n39。 )。 ylabel(39。 h(n)39。 )。 title(39。 冲激响应39

27） 比较方程（ 217）与（ 218），可知（ 219）中的位移 x 是（ 220）中复数 x 的虚部，因此（ 225）的虚部就是方程（ 212）的特解，即有 )sin(   tBx （ 228） 其中 B 为振幅，  为相位差。 由式（ 226）、 223）及（ 224）得出稳态强迫振动有如下的基本特点： 辽宁工程技术大学毕业设计（论文） 13

易出错，结果不直观。 数字高通滤波器 的设计  ，  ，通带最大衰减为 1dB，阻带最小衰减为 100dB。 基于 MATLAB 的 IIR 数字滤波器的设计 13 用 脉冲响应不变法 设计 巴特沃斯高通滤波器 ， 程序 如下： Wp=*pi/T。 Ws=*pi/T。 %设置归一化通带和阻带截止频率 T=1。 %设置采样周期为 1 fs=1/T。 %采样频率为周期倒数 Ap=1。

一个滤波器的截止频率 Wn 和最小阶数 n。 调用形式为： R s )R p ,W s ,p,c h e b2 o r d ( WW n ][ n ,  （ 7） 通过上式求得滤波器的 n 和 Wn 后，在通过 MATLAB 软件中的 cheby2 函数的调用得到滤波器 传递函数的 分母系数 a 和 滤波器传递函数的分子 系数 b。 调用形式 为： W n)R s ,c he by2( n

入到科学研究及工程计算各个领域的重要原因。 （ 3）强大的科学计算机数据处理能力 Matlab 是一个包含大量计算算法的集合。 其拥有 600 多个工程中要用到的数学运算函数，可以方便的实现 用户所需的各种计算功能。 函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。 在通常情况下，可以用它来代替底层编程语言，如 C 和 C++。 在计算要求相同的情况下，使用

index1=find(x1==0)。 %找出 x1数组中为 0的元素 x1(index1)=ones(1,length(index1))。 %返回一个 index1 x2=rand(1,length(t))。 x2=round(x2)。 index1=find(x2==0)。 x2(index1)=ones(1,length(index1))。 x3=rand(1,length(t))。