非线性方程组建立与求解问题。
显示 更早的评论
小弟在编一个在function里求解非线性方程组的问题，遇到一点不懂的地方求大神指导。
问题简化如下：
function result=fun(n)
equ1='p1=1*(x1+1)^2';
equ2='p2=2*(x1+1)^2+1*(x2+2)^2';
equ3='p3=3*(x1+1)^2+2*(x2+2)^2+1*(x3+3)^2';
equ4='p4=4*(x1+1)^2+3*(x2+2)^2+2*(x3+3)^2+1*(x4+4)^2';
......
equn=..........
result=solve(equ1,equ2,equ3,.....equn)
end
想通过此程序求解p1~pn,x1~xn,由于这里是简化后的问题，方程不封闭这儿不讨论是否可以解出。
小弟想请教，在n很大的情况下，就会有很多方程，像我这样一个一个的输入会很麻烦，solve()内的东西也会很冗长，想请教各位大神，可不可以通过循环或其它方法将该程序简化。
采纳的回答
wafij
2022-11-24
clear
clc
n = 5;%方程个数
for in = 1:n
eval(['p' num2str(in) '=0;']);
for ix = 1:in
eval(['syms x' num2str(ix) ';']);
eval(['p' num2str(in) '= p' num2str(in) '+' num2str(in+1-ix) '*(x' num2str(ix) '+' num2str(ix) ')^2' ';']);
eval(['equ{' num2str(in) '}=char(p' num2str(in) ');'])
end
equ{in} = ['p' num2str(in) '=' equ{in}];
end
solve(equ{:})
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 大型文件和大数据 的更多信息
