matlab fsolve函数 误差,matlab关于fsolve函数的运用出现的问题
在网上查找了很多相同问题的答案,但是都没有合适的方法解决我的问题,求大神们指点
K1=442.3926125;K2=45483.17187;K3=0.000253321;K4=22987873.07;K5=6.301410878;
K6=2.537198529;K7=3.82405E+17;K8=1.79031E+23;K9=1.65563E-22;K10=5.751160296;
K11=115.997776;K12=2819719195;K13=16056077273;K14=2.107894631;K15=1.36646E+12;
K16=1.0507E+15;K17=8785.72203;K18=238699.3468;K19=585427.5027;K20=294600207.6;
F=@(x)[(0.5*x1+K1*x1*x2+2*K2*x1^2*x2+3*K3*x1^3*x2+3*K4*x1^3*x2^2+2*K7*x1^2*x6+3*K8*x1^3*x6+4*K9*x1^4*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+2*K19*x1^2*x2^2*x3+3*K20*x1^3*x2^2*x3)*x7-0.2;
(x2+K1*x1*x2+K2*x1^3*x2+2*K4*x1^3*x2^2+K5*x2*x5+K6*x2*x5^2+K10*x2*x3+K11*x2*x3^2+K14*x2*x4^2+K17*x1*x2*x3+2*K18*x1*x2^2*x3+2*K19*x1^2*x2^2*x3+2*K20*x1^3*x2^2*x3)*x7-0.172;
(0.5*x3+K10*x2*x3+2*K11*x2*x3^2+2*K15*x3^2*x6+3*K16*x3^3*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+K19*x1^2*x2^2*x3+K20*x1^3*x2^2*x3)*x7-0.05;
(0.5*x4+3*K12*x4^3*x6+4*K13*x4^4*x6+2*K14*x2*x4^2)*x7-0.41;
(0.5*x5+K5*x2*x5+2*K6*x2*x5^2)*x7-0.04;
(x6+K6*x2*x5^2+K7*x1^2*x6+K8*x1^3*x6+K9*x1^4*x6+K12*x4^3*x6+K13*x4^4*x6+K15*x3^2*x6+K16*x3^3*x6)*x7-0.03;
x1+x2+x3+x4+x5+x6+K1*x1*x2+K2*x1^2*x2+K3*x1^3*x2+K4*x1^3*x2^2+K5*x2*x5+K6*x2*x5^2+K7*x1^2*x6+K8*x1^3*x6+K9*x1^4*x6+K10*x2*x3+K11*x2*x3^2+K12*x4^3*x6+K13*x4^4*x6+K14*x2*x4^2+K15*x3^3*x6+K16*x3^3*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+K19*x1^2*x2^2*x3+K20*x1^3*x2^2*x3-1];
x0=[0.7;0.3;0.5;0.7;0.2;0.5;0.5];
x=fsolve(F,x0)
报错
Undefined function or variable 'x1'.
Error in
@(x)[(0.5*x1+K1*x1*x2+2*K2*x1^2*x2+3*K3*x1^3*x2+3*K4*x1^3*x2^2+2*K7*x1^2*x6+3*K8*x1^3*x6+4*K9*x1^4*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+2*K19*x1^2*x2^2*x3+3*K20*x1^3*x2^2*x3)*x7-0.2;(x2+K1*x1*x2+K2*x1^3*x2+2*K4*x1^3*x2^2+K5*x2*x5+K6*x2*x5^2+K10*x2*x3+K11*x2*x3^2+K14*x2*x4^2+K17*x1*x2*x3+2*K18*x1*x2^2*x3+2*K19*x1^2*x2^2*x3+2*K20*x1^3*x2^2*x3)*x7-0.172;(0.5*x3+K10*x2*x3+2*K11*x2*x3^2+2*K15*x3^2*x6+3*K16*x3^3*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+K19*x1^2*x2^2*x3+K20*x1^3*x2^2*x3)*x7-0.05;(0.5*x4+3*K12*x4^3*x6+4*K13*x4^4*x6+2*K14*x2*x4^2)*x7-0.41;(0.5*x5+K5*x2*x5+2*K6*x2*x5^2)*x7-0.04;(x6+K6*x2*x5^2+K7*x1^2*x6+K8*x1^3*x6+K9*x1^4*x6+K12*x4^3*x6+K13*x4^4*x6+K15*x3^2*x6+K16*x3^3*x6)*x7-0.03;x1+x2+x3+x4+x5+x6+K1*x1*x2+K2*x1^2*x2+K3*x1^3*x2+K4*x1^3*x2^2+K5*x2*x5+K6*x2*x5^2+K7*x1^2*x6+K8*x1^3*x6+K9*x1^4*x6+K10*x2*x3+K11*x2*x3^2+K12*x4^3*x6+K13*x4^4*x6+K14*x2*x4^2+K15*x3^3*x6+K16*x3^3*x6+K17*x1*x2*x3+K18*x1*x2^2*x3+K19*x1^2*x2^2*x3+K20*x1^3*x2^2*x3-1]
报错为
Error in fsolve (line 218)
fuser = feval(funfcn{3},x,varargin{:});
Error in gongcun (line 13)
x=fsolve(F,x0)
Caused by:
Failure in initial user-supplied objective function evaluation. FSOLVE cannot
continue.
请问是格式有误吗?原来用syms solve函数解不出来,然后用fsolve数值函数迭代求解,困扰啊
matlab fsolve函数 误差,matlab关于fsolve函数的运用出现的问题相关推荐
- matlab浮点运算的误差,MATLAB的eps,浮点数相对精度
MATLAB的eps函数,官方帮助文档是这么描述的:Floating-point relateive accuracy,也就是浮点数的相对精度. 大家都知道,在数学中,实数有无穷多个,数值可以无限大, ...
- 用matlab求残余误差,matlab在测量误差分析中的应用
matlab在测量误差分析中的应用 MATLAB在测量误差分析中的应用 在技术测量中,按照误差的特点与性质,误差可分为:系统误差,粗大误差和随机误差.在假定不含有系统误差的情况下,可借助MATLAB对 ...
- matlab拟合公式不准确,Matlab拟合函数误差:函数值和YDATA大小不相等
我正在尝试使用fit函数来估计4参数模型(P B A R)并在消息后遇到错误,我不知道它是什么意思. Error using fit>iFit (line 367) Function value ...
- matlab指令vpa(j10),matlab中vpa函数
Matlab中矩阵函数_IT/计算机_专业资料.Matlab中矩阵函数 矩阵转置... Matlab 中 solve 函数主要是用来求解线性方程组的解析解或者精确解.对于得 出的结果是符号变量,可以通 ...
- matlab 解函数方程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
- matlab中sum参数,matlab中sum函数的用法
答:两个sum连起来可以是答案更精准,以下是用法:a=sum(x);%列求和a=sum(x,2)... sum 函数是求和函数.它是excel函数中最为常用的函数之一,sum函数分别出现在数学函数.全 ...
- matlab解方程教程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
第7章 MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7.4 函数极值 7.1 线性方程组求解 7.1.1 直接解法 1.利用 ...
- matlab求函数极值教程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
- MATLAB 求函数极值的内置函数一览表(实则优化算法函数汇总)
MATLAB 求函数极值的内置函数一览表 收集了几乎所有的 MATLAB 内置的优化函数,可收藏,需要时查阅. 文章目录 MATLAB 求函数极值的内置函数一览表 简介 概览 求函数极值(或最大值最小 ...
最新文章
- Ubuntu 14.04 安装 Sublime Text 3
- 小小串联电阻,大大的作用。常见串联电阻情况分类整理汇总
- cheat sheet for cs221 at stanford: good luck
- PST转换软件 v6.3
- 虚拟无线接入网:行业的演进方向
- 关于vmware中的centos无法ping通宿主机或者是外网的解决思路
- 是否允许此网站打开你计算机上的程序
- [软件更新]迅雷v5.9.8.1084发布
- 华为USG防火墙配置命令
- addrule android用法,RelativeLayout.LayoutParams 使用addRule失效的问题解决办法
- angular2完整视频教程收藏(大漠穷秋)
- 谈谈我对 iPhone5 全景照相机 的理解
- 【巴什博弈】HDOJ2188悼念512汶川大地震遇难同胞——选拔志愿者
- iOS Tagged Pointer
- 大学四年如何规划之出国留学
- oracle中修改用户密码
- 谷粒商城--后端笔记(混乱)
- Python基础——流程控制之循环语句
- Vitalik:为何需要区块链投票?
- BetaFlight模块设计之三十五:RSSI信号强度链路稳定性分析
热门文章
- Android 模拟器下使用虚拟摄像头
- WEB安全白帽(一) 米斯特
- (转)完美解决 Android WebView 文本框获取焦点后自动放大有关问题
- TI OMAP Bootloader Project 关于启动过程
- 庖丁解牛linux内核视频下载,庖丁解牛 - dataCube的个人空间 - OSCHINA - 中文开源技术交流社区...
- 移花接木—让Windows启动无限精彩
- 如何彻底卸载百度卫士?
- USB 驱动程序介绍(一)
- WINPE U盘版制作-深山红叶版,完美者等
- 關於RemotelyAnywhere 1326登陸錯誤解決方法