中易网

求高手!!!哪里错了???matlab的

答案:2  悬赏:40  
解决时间 2021-04-10 10:54
储存方程组编程如下将fun.m保存到工作路径中:
function f=fun(x);
%定义非线性方程组如下
%变量x1 x2
%函数f1 f2
syms x1 x2
f1=4*x1-x2+1/10*exp(x1)-1;
f2=-x1+4*x2+1/8*(x1)^2;
f=[f1 f2];
建立函数dfun
用来求方程组的雅克比矩阵将dfun.m保存到工作路径中:
function df=dfun(x);
%用来求解方程组的雅克比矩阵储存在dfun中
f=fun(x);
df=[diff(f,'x1');diff(f,'x2')];
df=conj(df');
编程牛顿法求解非线性方程组将newton.m保存到工作路径中:
function x=newton(x0,eps,N);
con=0;
%其中x0为迭代初值eps为精度要求N为最大迭代步数con用来记录结果是否收敛
for i=1:N;
f=subs(fun(x0),{'x1' 'x2'},{x0(1) x0(2)});
df=subs(dfun(x0),{'x1' 'x2'},{x0(1) x0(2)});
x=x0-f/df;
for j=1:length(x0);
il(i,j)=x(j);
end
if norm(x-x0) con=1;
break;
end
x0=x;
end
%以下是将迭代过程写入txt文档文件名为iteration.txt
fid=fopen('iteration.txt','w');
fprintf(fid,'iteration');
for j=1:length(x0)
fprintf(fid,' x%d',j);
end
for j=1:i
fprintf(fid,'\n%6d ',j);
for k=1:length(x0)
fprintf(fid,' %10.6f',il(j,k));
end
end
if con==1
fprintf(fid,'\n计算结果收敛!');
end
if con==0
fprintf(fid,'\n迭代步数过多可能不收敛!');
end
fclose(fid);
运行程序
在matlab中输入以下内容
newton([0.1 0.1 -0.1],0.00001,20)

上面的输入是错误的,是在matlab中输入以下内容
newton([0.1 0.1],0.00001,20)
最佳答案
错误提示是什么??
全部回答
ybest是不是拟合处理后的数据向量
我要举报
如以上问答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
十万元的轿车是不是看上去很垃圾?质量很差?
东莞长安哪里房子比较便宜
mos管是三极管吗
男女都属虎,同岁,没错开100天,2017年结婚
青春岛社区东区地址有知道的么?有点事想过去
缇奇米克穿越海贼王cp多弗朗明哥的小说叫什么
搬运工还需要面试,这个是面试什么的?会被淘
QQ中心发来的邮件收不到怎么回事
本月取得营业收入29000元,会计分录怎么做?
我的耳机插上出现个前面板麦克风已被选为录音
杭州市下城区地铁建设办公室怎么去啊,有知道
求 劲舞团 好看的舞步一套
如果受到刑讯逼供,受害人应该怎么做
请问谁有《终极23000单词放口袋》的mp3吗?很
魔兽世界法师T6和S4相比那个能好点
推荐资讯
西安碑林有王羲之的兰亭序吗
外交部的工作人员如果想要结婚的话,结婚对象
穿戴有钢丝文胸很痛是怎么回事阿?
安装VB未成功后出了个cd驱动器(z:)
火箭发射时,火箭尾部线下喷出高温气体,火箭
我下载了CF3.9到4.0的怎么还是进不去已经更新
我是真心的!~!~朱洋洋
两眉毛之间长粉刺是怎么回事
隔岸路2号院这个地址在什么地方,我要处理点
请问重庆农行卡在厦门取300扣多少?
DNF行家开价
高尔夫二手球怎么样翻新
手机登qq时,显示手机磁盘不足,清理后重新登
刺客的套装怎么选啊?