| 查看: 1370 | 回复: 4 | ||
wndg木虫 (正式写手)
|
[求助]
牛顿迭代法解非线性方程组中的矩阵维数 已有1人参与
|
|
f1.m %牛顿迭代法的 方程函数 function f=f1(x0) x=x0(1); y=x0(2); f1=x^2-2*x-y+0.5; f2=x^2+4*y^2-4; %最后方程函数 以行向量输出 f=[f1 f2]; df1.m function f=df1(x0); x=x0(1); y=x0(2); f=[2*x-2 -1 2*x 8*y]; new_ton.m %牛顿迭代法 计算非线性方程组 %输入 x0为迭代初值 %tol为误差容限 如果缺省 默认为10的-10次方 %data用来存放计算的中间数据便于计算收敛情况分析 function [x,n,data]=new_ton(x0,tol) if nargin==1 tol=1e-10; end x1=x0-f1(x0)/df1(x0); n=1; %迭代过程 while (norm(x1-x0)>tol)&(n<500) x0=x1; x1=x0-f1(x0)/df1(x0); n=n+1; %data用来存放中间数据 data(:,n)=x1; end x=x1; new_main.m %牛顿迭代法的主函数 x0=[1 1]; [x,n,data]=new_ton(x0); disp('计算结果为') x disp('迭代次数为') n %抽取data1中第一个变量数据 画出曲线 subplot(2,1,1) plot(data(1,  ),title('x在迭代中的变化')%抽取data中的第二个变量数据 画出其变化曲线 subplot(2,1,2) plot(data(2, ),title('y在迭代中的变化')%以下为数据存储部分 num=(1:n)'; a=[num data']; save data1.txt a -ascii 以上是书上的代码，我的问题是这个式子 f1(x0)/df1(x0) 左边是行向量，右边是2阶矩阵，为什么能除呢？？ 不应该是2维列向量才能除吗？ 对于AX=B,B不应该是列向量吗？？ 代码里X也是行向量，搞不清楚这种处理什么意思，求大神  |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有14人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
极限拉扯，如何抉择
已经有15人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 注册化工工程师须知
已经有23人回复
匿名
用户注销 (职业作家)
Matlab专家
- 应助: 223 (大学生)
- 金币: 2891.3
- 散金: 5302
- 红花: 59
- 沙发: 10
- 帖子: 3308
- 在线: 945小时
- 虫号: 0
- 注册: 2012-05-01
- 性别: GG
- 专业: 力学
2楼2016-10-13 17:56:01
昊客天下
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 184.8
- 散金: 5
- 帖子: 27
- 在线: 8.4小时
- 虫号: 4046101
- 注册: 2015-08-31
- 专业: 凝聚态物性I:结构、力学和
3楼2016-10-13 18:01:15
wndg
木虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 4679.8
- 散金: 1766
- 红花: 4
- 帖子: 620
- 在线: 154.8小时
- 虫号: 2609331
- 注册: 2013-08-21
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2016-10-13 18:05:17
匿名
用户注销 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1040.9
- 散金: 1688
- 红花: 21
- 沙发: 1
- 帖子: 1040
- 在线: 270.6小时
- 虫号: 0
- 注册: 2014-06-06
- 性别: MM
- 专业: 管理科学与工程
 本帖仅楼主可见 |
5楼2016-10-17 18:53:19
10