matlab组织的培训讲义,matlab培训讲义.doc
matlab培訓講義.doc
x-0.6*sin(x)-0.3*cos(y) );f{2}=inline( y-0.6*cos(x)+0.3*sin(y) );holdonh=arrayfun(@(i)ezplot(f{i},[-1,1,-1,1]),[1,2]);set(h(1), color , r )-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81xyy-0.6cos(x)+0.3sin(y)=0X=0.62249Y=0.38153Level=0.0057698functiony=myfun(x)y(1)=x(1)-0.6*sin(x(1))-0.3*cos(x(2));y(2)=x(2)-0.6*cos(x(1))+0.3*sin(x(2));end隱函數的求解給定一個隱函數的自變量,通過數值方法求解得到因變量。輸入一個數組:>>f=@(x)@(y)(exp(y)+x.^y).^(1/y)-x.^2.*yf1=@(x)arrayfun(@(a)fsolve(f(a),rand),x)f=@(x)@(y)(exp(y)+x.^y).^(1/y)-x.^2.*yf1=@(x)arrayfun(@(a)fsolve(f(a),rand),x)f1就是一個匿名函數,輸入可以是一個數組。帶參數的非線性方程的求解a=[1368532112];b=-fliplr(a);xy=zeros(2,length(a));index=1:length(b);m=arrayfun(@(t)fsolve(@(x)myfun1(x,a(t),b(t)),xy(:,t)),index, UniOutput ,false);微分方程常微分方程求解:解析解的求解,用dsolve函數。基本用法:dsolve( eq1 , eq2 ,., cond1 , cond2 ,., v )具體可以查看下幫助文件。當微分方程不能得到解析解,matlab有許多函數求解微分方程得到數值解。函數:ode45基本用法:[T,Y]=solver(odefun,tspan,y0,options)下面看一個例子:functiondy=myfun_3(x,y)dy=zeros(2,1);dy(1)=cos(y(1)-sin(x+y(2))+dy(2))-sin(2*x/y(1)+y(2));dy(2)=-2*x*y(2)+y(1)+sin(x-dy(1))*y(1);clc;clear;closealloptions=odeset( RelTol ,1e-4, AbsTol ,[1e-41e-4]);[T,Y]=ode45(@myfun_3,[02.5],[10.25],options);plot(T,Y(:,1), r*- ,T,Y(:,2), b+- )00.511.522.50.20.40.60.811.21.41.61.8二階微分方程的求解(幫助文件exmple2)對于多階方程的求解主要是自己使用相關的算法(歐拉法、二階或四階龍格庫塔法求解得到數值解)Ode45ode23等帶參數的微分方程clear;clc;closealllorenzeq=@(t,x,beta,rou,sig)[-beta*x(1)+x(2)*x(3);rou*x(2)+rou*x(3);-x(1)*x(2)+sig*x(2)-x(3)];beta=2;rou=5;sig=20;t_final=100;x0=[001e-10];[t,x]=ode45(@(t,x)lorenzeq(t,x,beta,rou,sig),[0,t_final],x0);plot(t,x);figure;plot3(x(:,1),x(:,2),x(:,3));axis([1042-2020-2025])其他類型?單個高階微分方程?高階常微分方程?隱式微分方程?微分代數方程?延遲微分方程?邊值問題?偏微分方程?Simulink求解?見《高等應用數學問題的matlab求解》擬合線性回歸非線性擬合線性回歸Regress查看幫助文件非線性擬合CftoolSftool參考matlab技術論壇電子期刊第三期查看一個簡單的例子幫助文件functionyhat=hougen(beta,x)b1=beta(1);b2=beta(2);b3=beta(3);b4=beta(4);b5=beta(5);x1=x(:,1);x2=x(:,2);x3=x(:,3);yhat=(b1*x2-x3/b5)./(1+b2*x1+b3*x2+b4*x3);loadreactionbeta=nlinfit(reactants,rate,@hougen,beta);積分?quad?quadl?quad2d
創作挑戰賽新人創作獎勵來咯,堅持創作打卡瓜分現金大獎總結
以上是生活随笔為你收集整理的matlab组织的培训讲义,matlab培训讲义.doc的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 电脑版word文档怎么下载(word文档
- 下一篇: matlab 树状链表,多级树集合分裂(