| 查看: 1326 | 回复: 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也是行向量，搞不清楚这种处理什么意思，求大神  |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 注册化工工程师须知
已经有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
