&& 查看话题
用MATLAB解方程作图，请大神帮忙
f1='a=exp(25.311-4.606*b)+c*exp(30.723-4.606*b)+c^2*exp(36.365-4.606*b)+c^3*exp(42.146-4.606*b)+c^4*exp(47.097-4.606*b)+exp(2.989-2.303*b)+exp(-9.614)+exp(-38.896+2.303*b)+exp(-68.337+4.606*b)+d*exp(15.507-4.606*b)+exp(-38.506+2.303*b)+exp(-68.783+4.606*b)';
f2='m+n=c+exp(21.417-2.303*b)*c+exp(30.723-4.606*b)*c+2*c^2*exp(36.365-4.606*b)+3*c^3*exp(42.146-4.606*b)+4*c^4*exp(47.097-4.606*b)';
f3='m=d*exp(15.507-4.606*b)+d+d*exp(23.782-2.303*b)+d*exp(134.483-4.606*b)';
f4='2*a+10^(-b)+c*exp(21.417-2.303*b)=2*m+10^(b-14)';
=solve(f1,f2,f3,f4)
上面是我写的解方程时的程序，f1、f2、f3、f4中是四个方程，其中abcdmn是未知数，我想通过解以上方程组得到m和n分别从0到10变化时，a关于m和n的三维图，如果得不到三维图，那就令n值固定（比如n=1）画出a关于m变化时的二维图也可以。
下面的附件是上面方程的来源
小白一枚，请能解决的大神把具体的过程写一下，不胜感激！！！
QQ截图06.png
MATLAB代码基本上就如下了，不过原方程不好解，初值难以给出。
function sol_three_D
format long
global m n
for n=0:1:10
& & for m=0:1:10
& && &x=fsolve(@nonlin,x0);
& && &a(m+1,n+1)=x(1);
=meshgrid(0:1:10,0:1:10);
surf(x,y,z),colorbar,%shading interp
function f= nonlin(x)
global m n
f1=exp(25.311-4.606*b)+c*exp(30.723-4.606*b)+c^2*exp(36.365-4.606*b)+...
& & c^3*exp(42.146-4.606*b)+c^4*exp(47.097-4.606*b)+exp(2.989-2.303*b)+...
& & exp(-9.614)+exp(-38.896+2.303*b)+exp(-68.337+4.606*b)+d*exp(15.507-4.606*b)+...
& & exp(-38.506+2.303*b)+exp(-68.783+4.606*b)-a;
f2=c+exp(21.417-2.303*b)*c+exp(30.723-4.606*b)*c+...
& & 2*c^2*exp(36.365-4.606*b)+3*c^3*exp(42.146-4.606*b)+...
& & 4*c^4*exp(47.097-4.606*b)-(m+n);
f3=d*exp(15.507-4.606*b)+d+d*exp(23.782-2.303*b)+d*exp(134.483-4.606*b)-m;
f4=2*a+10^(-b)+c*exp(21.417-2.303*b)-(2*m+10^(b-14));
f=; : Originally posted by 月只蓝 at
MATLAB代码基本上就如下了，不过原方程不好解，初值难以给出。
function sol_three_D
format long
global m n
for n=0:1:10
& & for m=0:1:10
& && &x=fsolve(@nonlin,x0);
& && &a( ... 再请问您一下，你里面的x0=是指的abcd的初值吗？这个初值的选取是由我具体的问题而定的吧？
我本来想通过分别确定m和n的值，这样就能解上面的方程组了，然后再用得到的点作图，我设m=1的话，让n从1到10变化，然后解这个方程，但是自己写的总是不对，能麻烦您再给指点一下吗？不然一个点一个点的做实在是得累死
小白一个，（你给的程序我都不会用，还是自己百度的:rol:）最后再麻烦一下大神！ : Originally posted by 一直在纠结 at
再请问您一下，你里面的x0=是指的abcd的初值吗？这个初值的选取是由我具体的问题而定的吧？
我本来想通过分别确定m和n的值，这样就能解上面的方程组了，然后再用得到的点作图，我设m=1的话，让n从1到10变化，然后 ... 确实是初值。大概算了一下，abcd的数量级差得很大，初值很难给出，一般的方程组这段程序就能把所有解和三维关系图给出。 : Originally posted by 月只蓝 at
确实是初值。大概算了一下，abcd的数量级差得很大，初值很难给出，一般的方程组这段程序就能把所有解和三维关系图给出。
... 想通过分别确定m和n的值，这样就能解上面的方程组了，然后再用得到的点作图，我设m=1的话，让n从1到10变化，然后解这个方程，但是自己写的总是不对，能麻烦您再给指点一下吗？不然一个点一个点的做实在是得累死，请大神再帮下忙吧 : Originally posted by 一直在纠结 at
想通过分别确定m和n的值，这样就能解上面的方程组了，然后再用得到的点作图，我设m=1的话，让n从1到10变化，然后解这个方程，但是自己写的总是不对，能麻烦您再给指点一下吗？不然一个点一个点的做实在是得累死，请 ... 上面的就是按你说的，让n从1到10变化求解的程序，只是原方程的初值难以给出，从而难以得到解。
另一个解决办法是用1stopt软件来解，该软件操作很简单，参见：
http://emuch.net/bbs/viewthread.php?tid=4586018&fpage=1 : Originally posted by 月只蓝 at
上面的就是按你说的，让n从1到10变化求解的程序，只是原方程的初值难以给出，从而难以得到解。
另一个解决办法是用1stopt软件来解，该软件操作很简单，参见：
http://emuch.net/bbs/viewthread.php?tid=4586018& ... 恩 好的 多谢啦！ : Originally posted by 月只蓝 at
MATLAB代码基本上就如下了，不过原方程不好解，初值难以给出。
function sol_three_D
format long
global m n
for n=0:1:10
& & for m=0:1:10
& && &x=fsolve(@nonlin,x0);
& && &a( ... 您好，帮忙看一下我的关于画图的MATLAB贴子好吗，谢谢啦！小弟比较菜，需要急用matlab里面怎么解方程？我是新手，比如说以下三个方程 x^2-4=0;3x+8=0;4x+2y=0，用什么命令来解？_百度作业帮
matlab里面怎么解方程？我是新手，比如说以下三个方程 x^2-4=0;3x+8=0;4x+2y=0，用什么命令来解？
matlab里面怎么解方程？我是新手，比如说以下三个方程 x^2-4=0;3x+8=0;4x+2y=0，用什么命令来解？
用solve， 你在主窗口里面打： help solve， 看看帮助就会了。比如：
[x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3 = 0')
如果把你这个看成方程组的话，两个未知数三个方程，解是求不出来的。一个一个的解了下：>> solve('x^2-4=0') ans =
>> solve('3*x+8=0') ans = -8/3
>> solve('4*x+2*y=0') ans = -1/2*y如何用matlab解方程
x=0.03-35.82ye^z x=0.028-0.58ye^2z x=0.026-0.41ye^3z
08-11-13 &匿名提问
MATLAB使用练习l        掌握MATLAB的基本使用方法实验1 数学建模初步l        通过实例了解数学建模的一般步骤;l        在以后的数学实验中用数学建模方法解决经过简化的实际问题;l        自觉培养用数学方法解决实际问题的意识和能力。实验2 差分方程与数值微分l         一阶常系数差分方程；l        高阶常系数差分方程；l         一阶常系数差分方程组；l      非线性差分方程；l        数值微分及其MATLAB实现；l        用差分方程或数值微分解决简单的实际问题。实验3 插值与数值积分   l        插值问题提法和求解思路;      l        Lagrange插值的原理和优缺点;      l 分段线性和三次样条插值的原理和优缺点;     l      用MATLAB实现分段线性和三次样条插值;    l      梯形、辛普森积分公式的原理及MATLAB实现l        数值积分公式的误差——收敛阶的概念l      高斯积分公式l        广义积分与多重积分l        用插值和数值积分解决简单的实际问题。实验4 常微分方程数值解l        欧拉方法的原理及龙格-库塔方法的思路l        局部截断误差和精度的概念l        龙格-库塔方法的MATLAB实现，包括求解微分方程组和高阶微分方程l        用微分方程解决简单的实际问题。实验5 线性方程组的解法l        主元素消去法和LU分解的原理l        方程组病态、向量和矩阵范数、条件数的概念l        迭代法的原理以及收敛的概念和条件l        用MATLAB解方程组，稀疏矩阵的处理l       拟合问题提法及最小二乘法的原理和结果;l        用线性方程组及最小二乘法解简单的实际问题。实验6 非线性方程近似解l        迭代法原理及收敛、收敛阶的概念l        用牛顿法解非线性方程和方程组l        非线性迭代法与混沌现象l        用非线性方程解简单的实际问题。实验7 无约束优化l        无约束优化模型及最优解的必要条件l        最速下降法、牛顿法、拟牛顿法的原理l        非线性最小二乘的解法l        MATLAB优化工具箱的用法，包括控制参数的功能，算法选择等l        用无约束优化（包括非线性最小二乘拟合）解决简单的实际问题。实验8 约束优化l        线性规划模型、解的性质和求解思路l        用MATLAB解线性规划，拉格朗日乘子的用途l        非线性规划模型、最优解的必要条件l        用MATLAB解非线性规划，包括控制参数的功能 l        用线性规划和非线性规划解决简单的实际问题实验9 整数规划l        整数规划模型、解的性质、松弛问题l        解整数规划的分枝定界法和动态规划法l        用LINDO解线性整数规划l        用LINGO解非线性整数规划l        用整数规划解决简单的实际问题实验10数据的统计描述和分析l        样本、频数、直方图、统计量（样本均值、标准差）的概念l        正态,  t, c2, F分布的由来，用MATLAB计算这些分布(密度、逆分布）l        正态总体下样本统计量的分布l        蒙特卡罗方法的原理及MATLAB实现l        用概率统计的基本思想解决简单的实际问题实验11 统计推断l        参数估计、置信区间、置信水平l        假设检验的思路，一总体和两总体均值的检验，双边和单边检验l        用MATLAB实现参数估计和假设检验l        实际问题的参数估计和假设检验实验12 回归分析l        回归分析要解决的问题（与数据拟合的关系）l        多元线性回归模型和系数的检验, 及MATLAB实现l        一元多项式、多元二项式回归, 逐步回归，非线性回归，及MATLAB实现l        用回归分析解决简单的实际问题实验13 数学建模综合l        综合应用各种方法，对实际问题建模、求解和分析
请登录后再发表评论!&& 查看话题
用matlab如何解方程？
请高手帮忙，请写出用matlab解下列方程的详细步骤，本人菜鸟一个，所以越详细越好，谢谢：
用 simulink画图吧 : Originally posted by sd116202 at
用 simulink画图吧 您好，谢谢您的回答，能给做一下吗？ 看看matlab里面的相关帮助文档吧，例子很详细的，应该养成看帮助文档的习惯
另外，有本matlab在微分方程求解方面的应用的书，http://emuch.net/bbs/viewthread.php?tid=571952
你这个问题应该是最简单的那种问题了:D http://ishare..cn/f/.html : Originally posted by SEEDEN at
http://ishare..cn/f/.html 谢谢提供资料，可惜该贴不能评分，下载好好看看 : Originally posted by mmccxx at
谢谢提供资料，可惜该贴不能评分，下载好好看看... 没事:D
顺便再说句，matlab的PDE工具箱使用时，要自己先把方程化成其要求的输入格式
可能在转化时有点小麻烦，比如算子符号——“倒三角”或者“正三角”，那个
这个找些相关的书看看（工程数学中的场论分析，一些数理方法的附录中也有），其实也简单，对照那些算子的相关公式变下形就好
跟matlab结合的一个功能强大的PDE求解软件是——COMSOL（最新版本要4个多G，有很多模型库），这个功能很强，使用方法跟matlab里的PDE工具箱有点像，只是更复杂罢了。如果有以后做相关的模拟仿真，可以关注一下这款软件:hand::D