matlab求数值解
用matlab 或mathematica求出积分微分的图或数值解
当然，最终结果并不是一个通常形式的解析解，但是只要把缺掉的参数赋上（具体说来就是a, b, c, d, f, V[0]，以下代码里的值是我随便赋的），就足以获得它的数值解了：eqn = a V'[u] - (b + c) V[u] + Integrate[(1 - d) c f Exp[-f x] V[u - x]\/(1 - d Exp@(...
用matlab求解二阶微分方程数值解,程序出现错误,求大神指点
1、把 G=1\/3*((5*Pp-2*P1)\/(P1-2*Pp)-P1*c1^2\/Pp*cp^2);改成 G=1\/3*((5*Pp-2*P1)\/(P1+2*Pp)-P1*c1^2\/(Pp*cp^2));表达式的错误导致出现高频振荡，积分步长非常小（大约在10^-8量级），所以很长时间算不出。2、从结果看，仿真时间取50秒就足够了（已进入稳态）。3...
如何用matlab求解微分方程并画图
如何用matlab求解微分方程并画图，可以先用dsolve（）或ode（）求出其微分方程（组）的解析解或数值解，然后用plot（）绘制其图形。例如：解微分方程 y'=y-2t\/y，y(0)=1，0<t<4 1、用dsolve（）求解，代码及结果如下 2、用ode45（）求解，代码及结果如下 3、当然喽，使用dsolve（）或ode（...
matlab怎么求解三角方程?
1、对于简单的三角方程，可以用solve()来求解，得到解析解。如 其求解代码和结果为 2、对于复杂的三角方程，可以用Broyden（拟牛顿法）等方法来求解，得到数值解。如 其求解过程如下 ①建立自定义函数文件，funm。m function y = funm(x)x=x(1,1)，y=x(2,1)y(1,1) = x(1,1) - 0.7...
matlab求方程数值解,最后得出结果总是有问题,请问有人帮我看一下吗...
由于方程有多个的解，使用subs（）无法用s数组值替代 x1、 x2、 x3、 y1、 y2、 y3、 z1、 z2、 z3的值。所以最有效的方法是，直接将s数组值赋值给 x1、 x2、 x3、 y1、 y2、 y3、 z1、 z2、 z3。按这个要求更改，可以得到a、b、c、d的解。syms a b c d s=[2,4,2^0....
matlab解三元一次方程组,数值解,fsolve()函数?
分析了题主的三元一次方程组，可以用vpasolve（）函数求解，而用fsolve（）来求解，需要找到k0，k1，k2附近值，否则可能得到的结果并不是你想要的结果。该三元一次方程组求解过程如下：m=1660;t=[14.34 17.43 21.49];delta_c=[1.04 1.05 1.05];V0=[34.72 30.56 26.39];V1=[29.17...
二元二阶微分方程组求数值解的matlab程序
方程或者初始条件有问题，请认真核实。原因：如果按照现在的初始条件，则第一个微分方程各项为0，方程不起作用。而且就算法而言，二阶（最高阶）导数的系数为0，会导致数值积分算法无法启动。参考代码大致如下：m=1;g=9.8;k=100;l0=1.1;x0=[0.1 0 0 0];% 定义x1=l, x2=l', x3=a, ...
用MATLABode45求微分方程 dy\/dx + ytanx = cosx 在[0,pi\/4]的数值解...
用MATLAB的ode45函数可以求得微分方程 dy\/dx + ytanx = cosx 在[0,pi\/4]的数值解。t0=0;tf=pi\/4;tspan=t0:pi\/40:tf;x0=[0,0]; %初值 [t,y]= ode45(@myodefun,tspan,x0); %myodefun为自定义微分方程函数 disp(' x y(x) dy\/dx')A=[t y];disp(A)figure(1)plo...
求一道用matlab编程解非线性方程组
对于非线性方程组F(X)=0，用fsolve函数求其数值解。fsolve函数的调用格式为： X=fsolve('fun',X0,option) 其中X为返回的解，fun是用于定义需求解的非线性方程组的函数文件名，X0是求根过程的初值，option为最优化工具箱的选项设定。最优化工具箱提供了20多个选项，用户可以使用optimset命令将它们显示...
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0, y(1)=2,y(2)=0,x...
这是你们老师出的题目吧，这个例子在好多书上都出现过。具体如下：function dydt = vdp1(t,y)dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];end 保存函数文件 [t,y] = ode45(@vdp1,[0 20],[2; 0]);This example uses @ to pass vdp1 as a function handle to ...
饶河县兰释回答： syms y1 y2; a=dsolve('Dy1=y2','(1+t^2)*Dy2=2*t*y2','y1(0)=1','y2(0)=3'); 运行结果: >> a.y1 ans = 3*t+t^3+1 >> a.y2 ans = 3+3*t^2 >> function dy=myfun(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=2*t*y(2)/(1+t^2); [t,y]=ode45(@myfun,[0 10],[1 3]) plot(t,y(:,1)); 网络故障,不小心发了两遍
饶庭18469222828问： matlab求数值解？
饶河县兰释回答：(1+x²)y''=2xy' y''/y'=2x/(1+x²) 两边积分 lny'=ln(1+x²)+c y'=c1(1+x²) y导|x=0=3 代入得 y'=c1*1=3 c1=3 y'=3(1+x²) 两边积分,y=3x+x^3+c2 y|x=0=1 代入得c2=1 因此特解,y=3x+x^3+1 希望能帮到你,祝学习进步
饶庭18469222828问： 用matlab解方程组的数值解,急 - ？
饶河县兰释回答： 解对初值很敏感, 需要提供初值试了试,当初值取很大时解收敛于V=0.072,R=35.622,C=1.2618
饶庭18469222828问： matlab中怎么求一元十次方程的数值解?列如:ax^10+bx^8+cx^6+dx^4+ex^2=0;其中:a,b,c,d,e均是已知的参数 - ？
饶河县兰释回答：[答案] 这个可以用solve这个函数来解 syms x; solve(a*x^10+b*x^8+c*x^6+d*x^4+e*x^2) 注意 要将上面的式子中的具体参数带进去,运行就可以了
饶庭18469222828问： Matlab 求数值解 - ？
饶河县兰释回答： ezplot('x-mod(x,100)-49/0.15',[0 1000])这是函数f(x)=x-mod(x,100)-49/0.15在0~1000范围里的图像可见,如果有数有解,解会出现在x=400附近这是函数值最接近0附件的图像ezplot('...
饶庭18469222828问： Matlab 求解方程组数值解dx/dt = yz;dy/dt = - xz;dz/dt = - 0.51xy;初始条件:x(0)=0 ,y(0) = 1 ,z(0) = 1 ,求0 - ？
饶河县兰释回答：[答案] clear;clc[T,Dxyz]=ode45('sdfd',[0 12],[0;1;1]);plot(T,Dxyz)--------------------------function dxyz=sdfd(t,x)dxyz=[0;0;0];dxyz(1) = x(2)*x(3);dxyz(2) = -x(1)*x(3);dxyz(3) = -0.51*x(1)*x(2);
饶庭18469222828问： 用MATLAB语言求微积分方程的数值解.(xd^2y)/dx^2 - 5dy/dx+y=0y(0)=0y'(0)=0 - ？
饶河县兰释回答：[答案] function dy=myfun03(x,y)dy=zeros(3,1) %初始化变量dydy(1)=y(2); %dy(1)表示y的一阶导数,其等于y的第二列值dy(2)=5/x*y(3)-y(1); %dy(2)表示y的二阶导数%ex0808 用ode23 ode45 ode113解多阶微分方程clear,clc[x23,y23]=ode23('myfun03',[1,10],[1 ...
饶庭18469222828问： matlab求方程数值解 - ？
饶河县兰释回答： 你的意思是不是这样: clear;a=1; for n=1:3x=eval(abs(solve('sqrt(x*x-a*a/4)+sqrt(1-a*a/4)=1','x'))) endx = 0.51760.5176
饶庭18469222828问： 用matlab 解这个方程的数值解,急 - ？
饶河县兰释回答： x=fsolve(fun,x0)求解fun(x)=0的解,x0是初值,fun是函数,x就是解 比如解方你的方程可定义一个函数 function y=f(x) y=(ln(0.2-x)-ln(0.04-x))/(ln(0.2-x)-ln(0.013-x))-0.55 再到主界面 x=fsolve(@f,0) 0是初始值 也可直接使用 fsolve(@(x)([(log(0.2-x)-log(0.04-x))/(log(0.2-x)-log(0.013-x))-0.55]),0)
饶庭18469222828问： 如何用matlab求解微分方程组的数值解 - ？
饶河县兰释回答： 在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; ...
