签到
苹果/安卓/wp
苹果/安卓/wp
客户端
0.0
0.00
推广加币
数据VIP
升级SVIP
注册
|
登录
项目交易
CDA数据分析师
CDA网校
CDA社区
CDA认证考试
CDA俱乐部
CDA Live
在线教育
JG学术培训
经管云课堂
CDA网校
CDA数据分析研究院
统计软件培训
金融科技
就学培训网
经管题库
培训证书查询
成为签约讲师
经管文库
专家入驻
学术博客
就学平台
美国在职研究生
论坛BBS
服务一览
VIP服务
数据VIP
数据库
兑换商城
广告服务
案例库
软件销售
校园代理
文献下载
会员课服务
我的
帖子
收藏
好友
我的空间
关注的人
关注的贴
找人
文库
任务
道具
勋章
网站地图
搜索
搜索
用户
人大经济论坛
›
标签
›
Start
标签: Start
经管大学堂:名校名师名课
相关帖子
版块
作者
回复/查看
最后发表
【求助】用R软件做copula的参数估计过程中出错,求解
R语言论坛
snowsummery
2013-1-13
13
10860
囊金如未足4
2018-8-28 20:14:51
Nature:Young scientists: Public engagement should start early
- [!reward_solved!]
求助成功区
qijiongli
2013-7-1
1
859
宝木0413
2013-7-1 16:48:13
using the format procedure to match userdefine data
SAS专版
zhou.wen
2013-6-14
0
1549
zhou.wen
2013-6-14 17:08:30
经典:Dale Carnegie - How To Stop Worrying And Start Living
商学院
zhang_nianbing
2006-12-17
1
2215
swjswj11
2013-1-21 23:51:31
[下载]venture start guide
商学院
MEI眼泪
2007-3-23
0
1327
MEI眼泪
2011-11-18 03:18:42
A company's cultural strategy: early or late start?
企业文化
yuanjinhui
2010-1-30
0
1465
yuanjinhui
2010-1-30 08:58:44
GS-US Portfolio Strategy-2010 Outlook-cyclical start,defensive finish-091208
行业分析报告
milfoil
2009-12-9
1
1110
yjwang05
2009-12-9 08:57:47
MS-China Coal Sector - Start of a New Upcycle-090619
金融学(理论版)
boeing
2009-6-25
0
1412
boeing
2009-6-25 16:08:07
[分享][下载]Statistical Analysis with R--a quick start
R语言论坛
smart_liu
2009-5-30
3
1935
xdl1125
2009-6-20 23:32:48
[建议]Five ways to start the world economic recovery
真实世界经济学(含财经时事)
yingtao5876
2008-12-22
0
2194
yingtao5876
2008-12-22 02:57:00
重金求QUICKBOOKS SOFTWATE,不要SIMPLE START EDITION
会计与财务管理
snowtea
2008-11-18
0
2628
snowtea
2008-11-18 13:37:00
BLOOMBERG_Quick_Start_Guide
金融学(理论版)
新生菜鸟
2008-11-1
0
2585
新生菜鸟
2008-11-1 10:24:00
[free]BLOOMBERG API for C Quick Start Guide
金融学(理论版)
cia801027
2008-2-28
0
4854
cia801027
2008-2-28 19:45:00
Malpe study guide of start
MATLAB等数学软件专版
ygh7814
2007-3-10
1
2146
ygh7814
2007-3-10 01:55:00
更多...
相关日志
分享
工作论文上传方法
冰霜雪莎
2015-7-9 04:31
https://hq.ssrn.com/login/pubSignInJoin.cfm? 先注册个号,填好信息。登陆以后左边选my papers然后点start submission按钮然后填表就好了 常见问题: http://www.ssrn.com/en/index.cfm/ssrn-faq/
0
个评论
分享
C++:变尺度法DFP求函数的最小值
accumulation
2015-4-24 16:53
#include "iostream.h" #include "math.h" #include "stdlib.h" double fun(double *x); double fcon(double *x); double funs(double *x,double M); void Grad(double (*pf)(double *x,double M),double m,int n,double *point,double *grad); double Search(double (*pf)(double *x,double M),double m,int n,double *start,double *direction); double DFP(double (*pf)(double *x,double M),double m,int n,double *min_point); double DFP(double (*pf)(double *x,double M),double m,int n,double *min_point) { int i,j,k; double e=1E-5,g_norm; double *g0=new double ; double *g1=new double ; double *dg=new double ; double *p=new double ; double t; double *x0=new double ; double *x1=new double ; double *dx=new double ; double **H=new double* ; for(i=0; in; i++) H = new double ; double **tempH=new double* ; for (i=0; in; i++) tempH = new double ; double *gH=new double ; double *Hg=new double ; double num1; double num2; for(i=0;in;i++) for(j=0;jn;j++) { if(i==j) H =1.0; else H =0.0; tempH =0.0; } for(i=0;in;i++) x0 =min_point ; Grad(pf,m,n,x0,g0); g_norm=0.0; for(i=0;in;i++) g_norm=g_norm+g0 *g0 ; g_norm=sqrt(g_norm); if(g_norme) { for(i=0;in;i++) min_point =x0 ; delete g1; delete p; delete x1; delete H ; delete tempH ; delete gH; delete =-g0 ; k=0; do { t=Search(pf,m,n,x0,p); for(i=0;in;i++) x1 =x0 +t*p ; Grad(pf,m,n,x1,g1); g_norm=0.0; for(i=0;in;i++) g_norm=g_norm+g1 *g1 ; g_norm=sqrt(g_norm); for(i=0;in;i++) { dx =x1 -x0 ; dg =g1 -g0 ; } for(i=0;in;i++) { gH =0.0; Hg =0.0; } for(i=0;in;i++) { for(j=0;jn;j++) { gH =gH +dg *H ; Hg =gH ; } } num1=0.0; num2=0.0; for(i=0;in;i++) { num1=num1+dx *dg ; num2=num2+gH *dg ; } for(i=0;in;i++) for(j=0;jn;j++) tempH =0.0; for(i=0;in;i++) { for(j=0;jn;j++) { tempH =tempH +H ; tempH =tempH +dx *dx /num1; tempH =tempH -Hg *gH /num2; } } for(i=0;in;i++) { for(j=0;jn;j++) { H =tempH ; } } for(i=0;in;i++) p =0.0; for(i=0;in;i++) { for(j=0;jn;j++) { p =p -H *g1 ; } } for(i=0;in;i++) { g0 =g1 ; x0 =x1 ; } k=k+1; }while(g_norme); for(i=0;in;i++) min_point =x1 ; delete g1; delete p; delete x1; delete H ; delete tempH ; delete gH; delete ; temp_point=new double ; Grad(pf,m,n,start,grad); diver_a=0; for(i=0;in;i++) diver_a=diver_a+grad *direction ; do { b=a+step; for(i=0;in;i++) temp_point =start +b*direction ; Grad(pf,m,n,temp_point,grad); diver_b=0; for(i=0;in;i++) diver_b=diver_b+grad *direction ; if( fabs(diver_b)1E-10 ) { delete temp_point; return b; } if( diver_b-1E-15 ) { a=b; diver_a=diver_b; step=2*step; } }while(diver_b=1E-15 ); for(i=0;in;i++) temp_point =start +a*direction ; value_a=pf(temp_point,m); for(i=0;in;i++) temp_point =start +b*direction ; value_b=pf(temp_point,m); do { s=3*(value_b-value_a)/(b-a); z=s-diver_a-diver_b; w=sqrt( fabs(z*z-diver_a*diver_b) ); t=a+(w-z-diver_a)*(b-a)/(diver_b-diver_a+2*w); value_b=pf(temp_point,m); for(i=0;in;i++) temp_point =start +t*direction ; value_t=pf(temp_point,m); Grad(pf,m,n,temp_point,grad); diver_t=0; for(i=1;i=n;i++) diver_t=diver_t+grad *direction ; if(diver_t1E-6) { b=t; value_b=value_t; diver_b=diver_t; } else if(diver_t-1E-6) { a=t; value_a=value_t; diver_a=diver_t; } else break; }while( (fabs(diver_t)=1E-6) (fabs(b-a)1E-6) ); delete temp_point; return t; } void Grad(double (*pf)(double *x,double M),double m,int n,double *point,double *grad) { double h=1E-3; int i; double *temp; temp = new double ; for(i=0;in;i++) { temp =point ; } for(i=0;in;i++) { temp +=0.5*h; grad =4*pf(temp,m)/(3*h); temp -=h; grad -=4*pf(temp,m)/(3*h); temp +=(3*h/2); grad -=(pf(temp,m)/(6*h)); temp -=(2*h); grad +=(pf(temp,m)/(6*h)); temp =point ; } delete ; double *x1=new double ; double e=1E-3,t; double M0,c; M0=1; c=2.5; for(i=0;in;i++) x0 =min_point ; k=0; do{ DFP(funs,M0,n,x0); for(i=0;in;i++) x1 =x0 ; t=fabs(fcon(x1)); if(te) { for(i=0;in;i++) min_point =x1 ; delete x1; return pf(min_point); } M0=c*M0; }while(k10); for(i=0;in;i++) min_point =x1 ; delete x1; return pf(min_point); } double fun(double *x) { //return -(0.5428-log(x )-((log(x )-5.9123)*(log(x )-5.9123))/0.0469-((x -127.2)*(x -127.2))/158.5624-0.1068*x -exp(-0.1068*(x -46.1355))); double a1=345,a2=181.52,b1=42.435,b2=10.70968; return log(2.51*b1*x )+0.5*((log(x )-a1)*(log(x )-a1))/(b1*b1)+log(2.51*b2)+0.5*((x -a2)*(x -a2)/(b2*b2)); } double fcon(double *x) { //return x -x -x ; return x -x ; } double funs(double *x,double M) { return fun(x)+M*fcon(x)*fcon(x); } void main() { int n=3; //double min_point ={201.15,201.15}; double mi ={5.836,181.52}; double min_value=Dev(fun,n,min_point); cout"本程序功能:使用变尺度法(DFP)求函数的最小值。"endl; cout"最小值点为:"endl; //cout"("min_point ","min_point ","min_point ")"endl; cout"("min_point ","min_point ")"endl; cout"最小值为:"endl; coutmin_valueendl; system("pause"); }
个人分类:
裂变模型
|
0
个评论
分享
A new start
w936528
2014-9-20 02:46
I still remember the first time I stepped on the land of America. I told myself it is a new start. I did it but I fail to believe it. Gradually, I realize it is not my fault but the nature of the world.
8 次阅读
|
0
个评论
更多...
京ICP备16021002-2号
京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明
GMT+8, 2024-4-25 10:13
积分 0, 距离下一级还需 积分