Matlab知识点总结.doc

1、符号积分改换傅里叶改换及其反改换1傅里叶变换f=f(x) F=F(w) syms x w u vf=sin(x)*exp(-xA2);F1=fourier(f)F1 = transform:fourier(sin(x)/exp(xA2), x, -w)>> f=x;F2=fourier(f)F2 = pi*dirac(w, 1)*2*i>> h=x*exp(-abs(x);F3=fourier(h)F3 = -(w*4*i)/(wA2 + 1)A2>> h=x*exp(-abs(x);F3=fourier(h,u)F3 =-(u*4*i)/(uA2 + 1)

2、A22傅里叶反改换syms w v x tg=exp(-abs(x);IF2=ifourier(g) IF2 = 1/(pi*(tA2 + 1)拉普拉斯改换及其反改换1拉普拉斯改换syms x s t vf1=sqrt(t);L1=laplace(f1)L1 =piA(1/2)/(2*sA(3/2)2拉普拉斯反改换syms a s t u v xf=exp(x/sA2);IL1=ilaplace(f)IL1 =ilaplace(exp(x/sA2), s, t)Z改换及其反改换方程的解析解线性方程组的解析解包含求解线性方程组和非线性方程组的函数solve (),也有求解常微分方程组的函数 ds

3、olve)L1='x+y+z=10'L2='3*x+2*y+z=14'L3='2*x+3*y-z=1' %L1 L2 L3分别是三个字符串g=solve(L1,L2,L3) g =x: 1x1 symy: 1x1 symz: 1x1 sym来明g是一个结构数组，其中每个元素为一>> g.x嗡号类型的量，用如下办法检查方程解的详细值 ans =1一般求解办法：L1='x+y+z=10'L2='3*x+2*y+z=14'L3='2*x+3*y-z=1'x y z=solve(L1,L2,L3

4、) x =1 y =2 z =7线性方程组的解析解>> f=sym('a*xA2+b*x+c=0');xf=solve(f) xf =-(b + 8八2 - 4*a*c)A(1/2)/(2*a)-(b - (bA2 - 4*a*c)A(1/2)/(2*a)>>syms y z u v weq1=u*yA2+v*z+w;eq2=y+z+w;y z=solve(eq1,eq2,y,z)y =(v + 2*u*w + (vA2 + 4*u*w*v - 4*u*w)A(1/2)/(2*u) - w(v + 2*u*w - (vA2 + 4*u*w*v - 4*u

5、*w)A(1/2)/(2*u) - w z =-(v + 2*u*w + (vA2 + 4*u*w*v - 4*u*w)A(1/2)/(2*u)-(v + 2*u*w - (vA2 + 4*u*w*v - 4*u*w)A(1/2)/(2*u)常微分方程组的解析解在微分方程组的表达式equ中，大写字母D表明对自变量 (设为x)的微分算子:D=d/dx,D2=d2/d2M分算子D后边的字母则表明为因变量，即带求解的不知道函 数。>> y=dsolve('Dy+a*x=0','x')y =C2 - (a*xA2)/2>> y=dsolve(&#

6、39;D2y+2*x=2*y','x')y = x + C4*exp(2A(1/2)*x) + C5/exp(2A(1/2)*x)>> y=dsolve('D2y+2*x=2*y','y(2)=5','Dy(1)=2','x')y =x + (exp(2A(1/2)*x)*(6*exp(2A(1/2) + 2A(1/2)/(2*exp(2A(1/2)*(exp(2*2A(1/2)+ 1) - (2A(1/2)*exp(2*2A(1/2)*(exp(2A(1/2)-3*2A(1/2)/(2*exp

7、(2A(1/2)*x)*(exp(2*2A(1/2) + 1)MATLAB!序设计全局变量global A B C变量名区别大小写脚本文件是m文件中最简略的一种输入顿号输出参数，用指令句子能够操控MATLA脂令作业空间的一切数据编程核算向量元素的平均值:x=input(输入向量：x=' ); %average_1dm向量元素的平均值m,n=size(x);if (m=1|n=1)average=sum(x)/length(x) else error仍须输入向量。')End假如m文件的第一个可执行句子以function 开端，该文件便是函数文件，每一 个函数文件都界说一个函数。函

8、数有自己独立的作业空间，它与 MATLAB的工 作空间区别开。Function y=average_2(X)f% average_2(X以计算向量元素的平均值%输入参数x为输入向量,MATLAB勺程序操控结构For循环一般方法For循环操控变量=表达式1:表达式2:表达式3For x=-2.0:0.25:-0.75While1环挑选结构：MATLAB中的挑选结构句子有if 句子、A=；for k=1:5for j=1:5if k=jA(k,k)=5;elseif abs(k-j)=1A(k,j)=1;elseA(k,j)=0;endendendA成果:5 1 0 0 01 5 1 0 00 1

9、 5 1 00 0 1 5 10 0 0 1 5程序流的操控1 .brea子停止本层for或while循环，跳转到本层循环句子输出参数y为核算的平均值switch 句子和try 句子end的下一条句子。2 .returr6j 子 停止被调用函数的运转，返回到调用函数。3 .pauseJ 子Pause暂停程序运转，按恣意键持续 Pause(n):程序暂停运转n秒后持续。Pause on/ of答应/制止这以后的程序暂停4 .continue)子在for循环或while 循环中遇到该句子，将越过这以后的循环体句子，进行下一次 循环。数据的输入与输出键盘输入句子inputx=input( '

10、 prom/1 );x=input( prompt 's') %不至于将输入的数字看成是数值型数据。屏幕输出句子disp(x)M数据文件的存储/加载(save/load)式化文本文件的存储/读取(fprintf/fxcanf二进制数据文件的存储/读取(fwrite/fread)：据文件行存储/读取(fgetl/fgets)MATLABC件操作MATLA敢据可视化运用函数figure来树立图形窗口，最简略的调用方法figureClos封闭当时图形窗口根本二维图形制作(1) plot x轴和y轴均为线性刻度(2) Loglog x轴和y轴均为对数刻度(3) Semilogxx为对数

11、刻度，y为线性刻度(4) Semilogyx为线性刻度，y为对数刻度(5) Plotyy :制作双纵坐标图形，制作两条具有不同纵坐标的曲线，调用格局Plotyy(x1, y1, x2, y0Plot是最根本的二维绘图函数：Plot(Y) , plot(X, Y)> > x=0:0.05:5;y=sin(x.八2);plot(x,y);> > x=0:0.05:5;y1=0.2*x-0.8;y2=sin(x.A2);plot(x,y1,x,y2);> > x=0:pi/180:2*pi;舷不同纵坐标的曲线，颜色不同加以区分y1=exp(-0.3*x).*cos

12、(2*x);y2=10*exp(-1.5*x);plotyy(x,y1,x,y2);条形图：> > x=-2.9:022.9;bar(x,exp(-x.*x);极坐标：t=0:0.1:2*pi;polar(t,abs(cos(2*t);针状图：>> x=0:0.1:4;y=(x.A0.8).*exp(-x);stem(x,y)阶梯图：>> x=0:0.25:10;stairs(x,sin(2*x)+sin(x)饼图：>> x=43,78,88,43,21;pie(x)颜色和线型>> x=0:0.2:8;y1=0.2+sin(-2*x)

13、;y2=sin(x.A0.5);plot(x,y1,'g-+',x,y2,'r-d');%y1 绿色、实线、加号标记；%y遗用赤色，虚线，菱形表明坐标轴及标示函数axis指定坐标轴的刻度规模其调用格局为Axis (xmin,xmax,ymin,yn> axhold on/of持原有图形做写原有图形axis on/of现/撤销坐标轴xlabel(' option ) x轴加标注，options示任意选项ylabel('option ) y轴加标注Title(" option )图形加标题Legend(option )图形加标注Gri

14、d on/of现撇销网格线Box on/of坐标加/不加边框线>>x=0:0.05:5;y1=exp(0.4.Ax)-1.5;y2=sin(x*4);plot(x,y1,x,y2, 'r-.');line(0,5,0,0);胭直线，代替横坐标xlabel( 'input' );ylabel( 'output');title( 'two function');legend( 'y1=exp(0.4Ax)-1.5' , 'y2=sin(x*4)')grid on %f网格线在一个图形窗口顶用

15、函数subplot能够一起画出多个子图形x=linspace(0,2*pi,100);subplot(2,2,1);plot(x,sin(x);xlabel( 'x' );ylabel( 'y' );title( 'sin(x)');subplot(2,2,2);plot(x,cos(x);xlabel( 'x' );ylabel( 'y' );title( 'cos(x)');subplot(2,2,3);plot(x,exp(x);xlabel( 'x' );ylabel( 'y' );title( 'exp(x)');subplot(2,2,4);plot(x,exp(-x);xlabel( 'x' );ylabel(