| 查看: 2214 | 回复: 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); |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有297人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
