|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
有关最优化函数fminunc的一个问题,谢谢大家了!
现在我将问题描述如下:
我要优化的能量函数V是坐标x,y,z的函数,其中坐标x,y,z又是碳纳米管原子编号的函数,例如,对于一个有25(b=25)层,每一层有14(a=14)个原子的碳纳米管,则坐标x(a,b),y(a,b),z(a,b)都是14行25列的矩阵,我求出了V对x,y,z的偏导数dV_x,dV_y,dV_z
请问如何利用最优化函数fminunc编写程序,谢谢!
这样对吗》?
function [f,g1,g2,g3] = myfun(x,y,z)
f=Vt(x,y,z);
g1=dVt_x;
g2=dVt_y;
g3=dVt_z;
options = optimset('GradObj','on');
x0;
y0;
z0; %初始坐标都给定了
[x,y,z] = fminunc(@myfun,x0,y0,z0,options);
x;
y;
z; |
|