| 查看: 2201 | 回复: 7 | ||
[求助]
求教matlab线性优化求最小值的问题
|
|
这两天遇到这样一个问题， 目前遇到一个函数 F(a,b,c,x)， 如果定义a,b,c，可以画出F在x一个区间上的曲线，但无法求得F的具体解析式。 我的目的是把F曲线与手头上的实验曲线达到最好吻合，同时求出相应的参数a,b,c。 目前我的思路是这样的：试验曲线有10000组数据 （x,y)，我用最小二乘的思想，求出每一个x上对应的F与y之差平方和，即sum((y-F)^2), 然后给出a, b , c的范围，用fmincon来求解，但是写了程序，发现问题，程序无法求出相应解。错误为 Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search). > In fmincon at 274 In main at 19 Optimization terminated: magnitude of directional derivative in search direction less than 2*options.TolFun and maximum constraint violation is less than options.TolCon. No active inequalities. 请问大神我该怎么办，是函数用错了吗？ 还是用别的方法？ 主程序如下： clear all;clc; global nr nt apt pidt vb freq nfrq ui=sqrt(-1); nr=50; nt=51; apt=100; pidt=10d0; vb=4212; %% A，b定义约束 45 < x1<55; 0.22 < x2<0.28; 0.0005 < x3<0.001; A=[1,0,0; 0,1,0;0,0,1;-1,0,0;0,-1,0;0,0,-1]; b=[55,0.28,0.001,-45,-0.22,-0.0005]'; x0= [49.13; 0.2592; 0.0007]; % Starting guess at the solution [x,fval] = fmincon(@myfun,x0,A,b) %%_______________________________________ function f = myfun(x); global nr nt apt pidt vb freq nfrq adf=load('y_para.dat'); fs=adf(:,1); con=adf(:,2); freq=fs; nfrq=length(freq); epson= 0.3302 ; ita=0.6052 ; c=1.0888 ; ym=admres(x(1),x(2),epson,ita,c,x(3)); f=sum((real(ym)-con).^2); |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有237人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab非线性拟合求助
已经有7人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- matlab无约束非线性方程求最小值的问题
已经有7人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- matlab解非线性方程组解不开
已经有3人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 用matlab拟合三元线性回归方程及分析
已经有9人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 关于matlab线性规划 背包问题的 解法 求救呀~~方程知道 但不知道如何解答····
已经有4人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- MATLAB非线性优化拟合怎么改才正确
已经有3人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab线性矩阵不等式求解做优化问题
已经有10人回复
- matlab求解非线性方程组
已经有16人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab非线性参数拟合问题
已经有7人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
2楼2013-05-14 15:54:36
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
3楼2013-05-14 16:01:07
4楼2013-05-14 21:35:05
5楼2013-05-14 21:36:51
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
6楼2013-05-14 21:41:15
7楼2013-05-16 15:06:56
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
8楼2013-05-16 16:10:41
