matlab命令大全

绘制正弦曲线和余弦曲线。

x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));求方程3x4+7x3 +9x2-23=0的全部根。

p=[3,7,9,0,-23]; %建立多项式系数向量x=roots(p) %求根用求特征值的方法解方程。

3x5-7x4+5x2+2x-18=0p=[3,-7,0,5,2,-18];A=compan(p); %A的伴随矩阵x1=eig(A) %求A的特征值x2=roots(p) %直接求多项式p的零点求积分quad('x.*log(1+x)',0,1)求解线性方程组。

a=[2,-3,1;8,3,2;45,1,-9];b=[4;2;17];x=inv(a)*babs函数可以求实数的绝对值、复数的模、字符串的ASCII码值。

用于取整的函数有fix、floor、ceil、round，要注意它们的区别。

rem与mod函数的区别。

rem(x,y)和mod(x,y)要求x,y必须为相同大小的实矩阵或为标量。

在MATLAB中，还可以用linspace函数产生行向量。

其调用格式为：linspace(a,b,n)其中a和b是生成向量的第一个和最后一个元素，n是元素总数。

显然，linspace(a,b,n)与a:(b-a)/(n-1):b等价。

常用的产生通用特殊矩阵的函数有：zeros：产生全0矩阵(零矩阵)。

ones：产生全1矩阵(幺矩阵)。

eye：产生单位矩阵。

rand：产生0～1间均匀分布的随机矩阵。

randn：产生均值为0，方差为1的标准正态分布随机矩阵。

分别建立3×3、3×2和与矩阵A同样大小的零矩阵。

(1) 建立一个3×3零矩阵。

zeros(3)(2) 建立一个3×2零矩阵。

zeros(3,2)(3) 设A为2×3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。

A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵Azeros(size(A)) %产生一个与矩阵A同样大小的零矩阵(1) 在区间[20,50]内均匀分布的5阶随机矩阵。

(2) 均值为0.6、方差为0.1的5阶正态分布随机矩阵。

命令如下：x=20+(50-20)*rand(5)y=0.6+sqrt(0.1)*randn(5)此外，常用的函数还有reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。

求(x+y)5的展开式。

在MATLAB命令窗口，输入命令：pascal(6)矩阵次对角线上的元素1,5,10,10,5,1即为展开式的系数。

点运算符有.*、./、.\和.^。

两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。

1．矩阵平方根sqrtmsqrtm(A)计算矩阵A的平方根。

2．矩阵对数logmlogm(A)计算矩阵A的自然对数。

此函数输入参数的条件与输出结果间的关系和函数sqrtm(A)完全一样求一元二次方程ax2 +bx+c=0的根。

程序如下：a=input('a=?');b=input('b=?');c=input('c=?');d=b*b-4*a*c;x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);计算分段函数的值。

x=input('请输入x的值:');if x<=0y= (x+sqrt(pi))/exp(2);elsey=log(x+sqrt(1+x*x))/2;endY先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。

程序如下：A=[1,2,3;4,5,6]; B=[7,8,9;10,11,12];tryC=A*B;catchC=A.*B;endClasterr %显示出错原因一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。

输出全部水仙花数。

for m=100:999m1=fix(m/100); %求m的百位数字m2=rem(fix(m/10),10); %求m的十位数字m3=rem(m,10); %求m的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend例3-14 编写函数文件求半径为r的圆的面积和周长。

function [s,p]=fcircle(r)%CIRCLE calculate the area and perimeter of a circle of radii r%r 圆半径%s 圆面积%p 圆周长%2004年7月30日编s=pi*r*r;p=2*pi*r;例3-15 利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。

function [rho,theta]=tran(x,y)rho=sqrt(x*x+y*y);theta=atan(y/x);调用tran.m的命令文件main1.m：x=input('Please input x=:');y=input('Please input y=:');[rho,the]=tran(x,y);rhothe利用MATLAB 产生指数序列 x [k ]=Kaku [k ]a = input('输入指数 a = ');K = input('输入常数K = ');N = input ('输入序列长度N = ');k = 0:N-1;x = K*a.^k;stem(k,x);xlabel('时间');ylabel('幅度');title(['\alpha = ',num2str(a)]);利用MATLAB 函数conv 计算两个序列的离散卷积。

x=[-0.5,0,0.5,1]; kx=-1:2;h=[1,1,1]; kh=-2:0;y = conv(x, h);k=kx(1)+kh(1):kx(end)+kh(end);stem(k,y);xlabel('k'); ylabel('y');的幅度频谱画出Ω-Ω-=j j e 11)e ( αX%Computation of Sequence DTFTb=[1];a1=[1 -0.9]; a2=[1 0.9];w=linspace(0,2*pi,512);X1=freqz(b,a1,w);X2=freqz(b,a2,w);plot(w/pi,abs(X1),w/pi,abs(X2),':');legend('\alpha=0.9','\alpha=-0.9');已知输入信号为角频率分别为0.1和0.4的离散正弦信号。

设计简单FIR 高通滤波器, 滤除低频分量，保留高频分量。

%信号通过简单FIR DF 的响应% Set up the filter coefficientsa=1;b=[-6.76195 13.456335 -6.76195];% Generate the two sinusoidal sequencesN=100; k = 0:N-1;x1 = cos(0.1*k);x2 = cos(0.4*k);%Generate the filter output sequencey = filter(b, a, x1+x2);plot(k,y,'r',k,x2,'b--',k,x1+x2,'k:');axis([0 100 -2 4]);ylabel('Amplitude'); xlabel('Time index k');legend('y[k]', ' x2[k]', ' x1[k]+x2[k]');试求下面系统函数的零极点形式二阶因子形式。

128.016.08.004.0)(233-+-+=z z z zz z H%Determination of the factored form and%the second order section form of a% rational z-transformb =[1 0 0.04 0];a =[1 -0.8 0.16 -0.128];[z,p,k]=tf2zp(b,a);disp('Zeros are at'); disp(z);disp('Poles are at'); disp(p);disp('Gain constant');disp(k);sos=zp2sos(z,p,k);disp('Second-order sections');disp(real(sos));已知一离散因果的LTI 系统的系统函数为3.0005.05.02)(2323+--++=z z z zz z z H 画出z 平面的零极点分布图。

%zeros and poles of the transfer functionb =[1 2 1 0];a =[1 -0.5 -0.005 0.3];zplane(b,a);已知一段语音信号中混入了一频率fn=1200Hz 正弦型噪音信号。

抽样频率fsam= 22050Hz 。

试用二阶带阻滤波器滤除语音信号中的噪音信号。

Fn=1200;W3dB=0.06;%读入语音信号[x,Fs] = wavread('sample2');%播放语音信号sound(x,Fs);pause;%设计滤波器W0=2*pi*Fn/Fs;beta=cos(W0);alpha=min(roots([1 -2/cos(W3dB) 1]));a=[1,-beta*(1+alpha),alpha];b=[1 -2*beta 1]*(1+alpha)/2;%信号滤波y=filter(b,a,x);%播放滤波后语音信号sound(y,Fs);4 ,16 ),/ π2cos(][===r N N rk k x 利用MATLAB 计算16点序列x [k ]的512点DFT N = 16;k = 0:N-1;L = 0:511; x = cos(2*pi*k*4./16);X = fft(x);plot(k/16,abs(X),'o');hold onXE = fft(x,512);plot(L/512,abs(XE))利用MATLAB 由DFT 计算x [k ]* h [k ]。