matlab怎么重新打开新的代码,方程求解程序代码求助-程序代码修改或新的代码...
很简单的方程求解程序,调用mulDNewton函数求解,之前在Matlab 2011b版本上运行成功,现在在Matlab 2018a版本上总是出错,程序代码和出错的提示如下,mulDNewton函数代码也如下。
请教大神该程序应该如何修改,请给出正确的代码,必有重谢!@所有人
mulDNewton迭代函数代码如下:
function [r,m]=mulDNewton(F,x0,eps)
%非线性方程组:F
%初始解:x0
%解的精度:eps
%求得的一组解:r
%迭代步数:m
if nargin==2
eps=1.0e-4;
end
x0 = transpose(x0);
dF = jacobian(F);
m=1;
tol=1;
while tol>eps
tol=1;
w=1;
Fx = subs(F,findsym(F),x0);
dFx = subs(dF,findsym(dF),x0);
F1=norm(Fx);
while tol>=0 %下面的循环是选取下山因子w的过程
r=x0-w*inv(dFx)*Fx; %核心的迭代公式
Fr = subs(F,findsym(F),r);
tol=norm(Fr)-F1;
w=w/2;
end
tol=norm(r-x0);
m=m+1;
x0=r;
if(m>10000) %迭代步数控制
disp('迭代步数太多,可能不收敛!');
return;
end
end
方程程序代码如下:
clc
clear
syms x y z
f1=x^2+y+z^2-12
f2=x^3+y^2+z-8
f3=x+y+z-6
x0=[1 1 1]% 初值
f=[f1;f2;f3]%方程组
[r,b]=mulDNewton(f,x0,1.0e-6)%调用迭代函数
%求解方程,并将结果赋予N01,N02,N03
N01=r(1);
N02=r(2);
N03=r(3);
% 求解的正确结果应该是X=1,y=2, z=3
错误提示代码如下:
警告: FINDSYM will be removed in a future release. Use SYMVAR instead.
> In sym/findsym (line 32)
In mulDNewton (line 19)
In equation (line 11)
错误使用 sym>convertChar (line 1448)
Character vectors and strings in the first argument can only specify a variable or number. To evaluate character vectors and strings representing symbolic
expressions, use 'str2sym'.
出错 sym>tomupad (line 1214)
S = convertChar(x);
出错 sym (line 211)
S.s = tomupad(x);
出错 sym/subs>normalize (line 168)
X = {sym(X)};
出错 sym/subs>mupadsubs (line 157)
[X2,Y2,symX,symY] = normalize(X,Y); %#ok
出错 sym/subs (line 145)
G = mupadsubs(F,X,Y);
出错 mulDNewton (line 19)
Fx = subs(F,findsym(F),x0);
出错 equation (line 11)
[r,b]=mulDNewton(f,x0,1.0e-6)%调用迭代函数
matlab怎么重新打开新的代码,方程求解程序代码求助-程序代码修改或新的代码...相关推荐
- matlab油膜压力程序,编的关于雷诺方程求解油膜压力的程序无法运行
程序是基于活塞环-缸套润滑模型的雷诺方程求解,采用有限差分法迭代求解,但程序无法在command window界面运行,求大家帮忙看一下 ita1=0.8*10^-6; ita2=0.4*10-6; ...
- Matlab 隐函数方程求解最小二乘法拟合一阶线性拟合二阶拟合传感器实验
九层妖塔 起于垒土 Matlab 最小二乘法拟合一阶线性拟合&传感器实验 一.代码 二.数据处理结果 三.Notes 一.代码 %电容传感器位移实验数据 最小二乘法一阶线性拟合 x = ...
- 使用MATLAB低版本打开高版本simulink文件
目录 (一)拟解决问题 (二)解决办法 (一)拟解决问题 由于MATLAB软件每年更新两次版本,所以当两个MATLAB软件版本不一致的时候,低版本的MATLAB在打开高版本创建的simulink文件时 ...
- 莱昂纳多(Leonardo)方程求解计算国土面积 (matlab程序答案)
上机题目 莱昂纳多(Leonardo)方程求解 莱昂纳多(Leonardo)在1225年研究了方程 x3+2x2+10x-20=0,并且得到它的一个根 x*=1.368808107,无人知道他是怎么做 ...
- 打开游戏要运行 19.8 亿次 if 语句?黑客嘲讽 RockStar 游戏代码太烂了
贾浩楠 发自 凹非寺 量子位 报道 | 公众号 QbitAI 博客原文 https://nee.lv/2021/02/28/How-I-cut-GTA-Online-loading-times-by- ...
- 用自然语言从GitHub搜代码,跳过论坛提问环节,来自Facebook新研究
晓查 发自 凹非寺 量子位 报道 | 公众号 QbitAI "如何关闭或隐藏Android软键盘?" 如果你是个Android入门开发者,去Stack Overflow论坛去寻找 ...
- matlab lyap,Matlab的Lyapunov、Sylvester和Riccati方程的Matlab求解
一.连续Lyapunov方程连续Lyapunov方程可以表示为 Lyapunov方程来源与微分方程稳定性理论,其中要求C为对称正定的n×n方阵,从而可以证明解X亦为n×n对称矩阵,这类方程直接求解比较 ...
- PCR主成分回归预测MATLAB代码 代码注释清楚。 main为主程序,读取EXCEL数据
PCR主成分回归预测MATLAB代码 代码注释清楚. main为主程序,读取EXCEL数据,也可以换自己数据集. 很方便,容易上手. ID:9624654486820873总有刁民膜拜朕
- lyapunov函数 matlab,科学网-[转载]Matlab的Lyapunov、Sylvester和Riccati方程的Matlab求解-吴雄君的博文...
一.连续Lyapunov方程连续Lyapunov方程可以表示为 Lyapunov方程来源与微分方程稳定性理论,其中要求C为对称正定的n×n方阵,从而可以证明解X亦为n×n对称矩阵,这类方程直接求解比较 ...
最新文章
- sql 集合查询 数据更新操作语句
- java8新特性_Java8新特性之Date API|乐字节
- formdata.append加多个值_呆哥数学每日一题 —— 多个最值函数混合向量
- 为修复一个代码执行安全漏洞,TensorFlow决定不再支持YAML
- 自动点击器如何设置最快_铁粉技巧 | iPhone如何设置自动开关机,iPhone更改字体...
- WinForm始终置顶并获取焦点
- HT中特殊字符φ的显示问题
- 下载 6g 概念及愿景白皮书_6G,到底有多6?6G概念及愿景白皮书正式发布!
- 太赞了!阿里巴巴AI每天服务全球10亿人
- php中json_encode中文编码问题(转)
- 李晨 | 无人机市场浅析
- echar地图使用小总结
- LUP分解求解线性方程组及求逆矩阵 java
- 大话西游2服务器修改,大话西游2:9.17维护解读:五倍次数修改全服上线,灵兽村要变样啦...
- 油猴插件的介绍和安装详解脚本的介绍和添加举例
- android提示“SIM not provisioned”
- 企业建网站力求大气上档次!
- 2022 年的前端行业,咋样啦?
- chrome浏览器手机模式下不显示鼠标
- android临时文件夹,【报Bug】安卓端,保存视频到相册成功,临时文件夹中有,但是相册中没有...
热门文章
- mysql sql实现原理_Mysql的数据库原理
- 信阳农林技术学院经纬度_信阳无人机创业团队:让梦想在田野起飞
- python计算入门_Python入门教程02-01(计算机核心基础)
- MATLAB GUI引用了不存在的字段
- windows7 安装MySQL 解决报错 服务没有响应控制功能
- Verilog HDL设计实现m序列+选择器
- 图像的旋转和尺度缩放在对数极坐标系下的研究
- 小甲鱼 OllyDbg 教程系列 (十六) : 简单病毒的逆向分析
- linux 禁止 密码 登陆,CentOS设置证书登录并禁止密码登录
- mysql profile 调试sql_SQL Server profile使用技巧