龙格库塔法matlab实例
鲍残妹2614龙格库塔法在matlab中能不能实现一阶常微分方程的计算? -
苏君徐15572108636 ______ ode45和ode23都是用龙哥库塔法求解常微分方程的命令,直接用即可,具体用法可参看其帮助文档.[t,y]=ode45(@(t,p)-3*p./(100+t),[0 200],7); plot(t,y);xlabel('t');
鲍残妹2614龙格库塔法求解一阶常微分方程组 的Matlab代码 -
苏君徐15572108636 ______ clc f=@(y,t)([1.65e-3*(14.93-y(1))-1.70e13*exp(-9064.23/y(2))*y(1);4.93-0.01655*y(2)+1.117e13*exp(-9064.23/y(2))*y(1)]); [Ca T]=ode45(f,[0 1],[.8 300])
鲍残妹2614Matlab用四阶龙格库塔法求解不可微分方程组的初值
苏君徐15572108636 ______ %% function dx=aaa(t,x)%定义matlab函数来描述方程 dx=[120-2*x(1)+2*x(2);2*x(1)-5*x(2)] %% x0=[0;0] [t,y]=ode45('aaa',[0,1],x0); plot(t,y)
鲍残妹2614急求MATLAB编程源代码用四阶龙格库塔法解如下微分方程 y'=y - 2x/y(0<x<1),y(0)=1,步长为h=0.2 -
苏君徐15572108636 ______ % 以下另存为文件 myrk4.m function [x,y]=myrk4(ufunc,y0,h,a,b)%参数: 函数名称,初始值向量,步长,时间起点,时间终点 n=floor((b-a)/h);%求步数 x(1)=a;%时间起点 y(:,1)=y0;%赋初值%按龙格库塔方法进行求解 for ii=1:n x(ii+1)=x(ii)+h; ...
鲍残妹2614matlab用四阶龙格库塔法解微分方程组; -
苏君徐15572108636 ______ function dy=test(t,y) dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);]; [t,y]=ode45('test',[0.01,1],[1,1]) plot(t,y(:,1),t,y(:,2)); x=3000*sin(70*2*pi/360)-y(1).*sin(y(2)); z=3000*cos(70*2*pi/360)-...
鲍残妹2614求编程达人帮忙用matlab编程用龙格库塔方法解微分方程 -
苏君徐15572108636 ______ clearclcn=4;h0=(1-0)/n;h=h0;m=0;y(1)=1;y(n+1)=2;y0(3)=33;n=n/2;while abs(y(n+1)-y0(3))>((1e-8)/15) n=n*2;x=0:h:1; for k=1:n...
鲍残妹2614用四阶龙格库塔法求解 -
苏君徐15572108636 ______ 初值给一下. 在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt=ode_...
鲍残妹2614用matlab编程实现四阶龙格库塔解二元二阶微分方程组 -
苏君徐15572108636 ______ 求解二阶微分方程,初始条件还需要给出y1'(0)和y2'(0).这里暂时按照0处理. function zd530003514 a=0.1; b=0.1; Y0 = [b-1; 0; b; 0]; % 解方程 [t,Y]= ode45(@ode,[0 10],Y0); y1=Y(:,1); y2=Y(:,3); % 绘图 subplot 211 plot(t,y1); subplot 212 plot(t,y2...
鲍残妹2614可不可以给给遗传算法的实例,有matlab实现的最好!谢谢咯! -
苏君徐15572108636 ______ 点matlab里面的start菜单,从toolbox里面选取optimization toolbox,在早期版本里面不叫这个名字,从2010b以后的版本把非线性方程组求解的算法都放在一起了.里面第一个下拉菜单选择最后一个,里面写的ga.fitness function是你的函数.遗传算法求解的goal是0,如果你要求趋近某一个值,就要看你的fitness function怎么写了.下面number of variables是你要求解的维数.接着点击运行就行了.因为不知道你的具体目的是什么,也没有代码发上来,所以我只能说这种很没用的话.
鲍残妹2614matlab算法 -
苏君徐15572108636 ______ 拉格朗日function y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s;end SOR迭代法的Matlab程序 function [x]=SOR_iterative(A,b)% 用SOR...