matlab的一些函数.doc

Matlab基本算术命令： 数组的计算： x=[3,5,11,17,8,9,31,40,12]; y=x+16 //数组里每个数加16 y=sqrt(x) //数组每个数开平方 n=length(x) //数组里偶数加2 for i=1:n if mod(x(i),2)==0 x(i)=x(i)+2 end end sum(x(1:2:length(x))) //数组里系数为奇数的数的和 ======================================== 数组的阶乘运算须知： 当x是一组数时，如：x=1:0.1:4 此时：y=x^3、y=exp(x*x) 会出错 应为：y=x.^3、y=exp(x.*x) 用“点乘” ======================================== 数组操作命令： x=[3,1,5,7,9,2,6]; 定义一个数组 x(3) 该数组的第三个数是？ x(1:7) （显示该数组第1到第7个数） x(1:end-1) 示该数组1~（7-1）个数 x(6:-2:1) 第六个数开始，从后往前隔着取数 sum(x) (数组内所有数求和) ======================================== 矩阵命令： A = [ 1 2 3; 3 4 5; 6 7 8] 按行写3X3矩阵 B = [ [1 2 3] [2 4 7] [3 5 8]] 按列写3X3矩阵 B = zeros (m,n) m X n零矩阵 E = ones (m, n) m X n单位矩阵 N = rand(m, n) 随即矩阵 A 矩阵转置 其他矩阵操作命令见第一次作业 ===================================== 控制系统模型相关命令： 用ode45 解常微分方程： 例题：求解dy/dt = -2y/t + 4t 在区间[1,3]的解 f=-2.*y./t+4*t //定义函数 [T,Y]=ode45(f,[1 3],2); //调用ode45命令 //[1,3]为区间，2为初值 Y_analysis = T.^2 + 1./(T.^2); // size(T); size(Y); size(Y_analysis); ? plot(T,Y,r,T,Y_analysis,g) ======================================== 传函模型： num1=[1] den1=[1000 50] sys=tf(num1,den1) ======================================== 分子或分母由多个多项式乘积时的传函模型： 例如：G(s) = 4(s+2)(s2+6s+6)/[s(s+1)3(s3+3s2+2s+5)] 用conv命令处理： num = 4*conv([1 2],[1 6 6]) den = conv([1 0],conv([1 1],conv([1 1],conv([1 1],[1 3 2 5])))) Sys=tf(num,den) 或者另一种方法： 例如 1/(((s+1)^2+s+2)^3+1) ? s=tf(s); //定义算子 G=1/(((s+1)^2+s+2)^3+1) 可以得出结果 ======================================== 零极点模型： h=zpk({[-3]},{[-1 -2 -5]},6) Zero/pole/gain: 6 (s+3) ----------------- (s+1) (s+2) (s+5) ======================================== 状态空间模型： ss 用法略 三种模型之间的转换： [A,B,C,D]= tf2ss([1],[1,100]) [z,p,k]=tf2zpk([3 23 46 22],[3 4 2 7 8 5]) 其它的看课件 ======================================== 两个传函之间的串联、并联与反馈： y1 y2 为两个不同的传递函数： y=series(y1,y2) 两模型串联 y=paralle(y1,y2) 两模型并联 yf=feedback(y,1) 加上单位负反馈后的传函 -1时为单位正反馈