相关日志
-
-
分享
金融小程序
-
accumulation 2017-1-21 20:06
-
1.德州扑克计算器 2.雪球炒股
-
个人分类: 金融学|0 个评论
-
-
分享
SAS程序：字符变量统一转数值变量
-
liehuodgq 2015-12-25 14:21
-
/****************************************************************************************** *用途 ： 将a表中的数值的字符变量全部转换成数值变量，并将变量名称还原 *适用 ： 表中必须含有id变量且为唯一号，id变量为字符类型 *Author : Flame *Date : 11dec2015 *******************************************************************************************/ data a; input id $ x $ y $ z w $ t $ o $ p $ q $; cards; 1 3 app 34 2 9 苹果 3 梨 4 6 people 56 8 9.6 4,8 cc 香蕉 12 6 t23 9 . . ztc4 4q 菠萝 15 89 pear 90 0 3 yourar er5 番茄 ;run; %macro vartonum(); proc sql noprint; select max(varnum) into :maxno from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='A'; QUIT; %put maxno; proc sql noprint; select name into :all_var_ separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='A' ; QUIT; %put all_var_; %let all_var=%sysfunc(compress(all_var_,id)); %put all_var; data b; set a;run; data E; set a;run; %do i=1 %to maxno; proc sql noprint; select name into :vnam from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='A' and varnum=i; QUIT; data a_0; set a; if anyfirst(vnam) ne 0 and anydigit(vnam) ne 1 then v_id=2; else if anyfirst(vnam) =0 and anydigit(vnam)=1 then v_id=1; else v_id=0; run; proc sql noprint; select avg(case when v_id=0 then . else v_id end) into :var_avg from a_0 ;quit; %put var_avg; data a; set a; if var_avg=1 then aa_i.=input(vnam,8.); run; data b; set b; if var_avg ne 1 then aa_i.=vnam.; run; %end; data d; set a; drop all_var_; run; proc sql noprint; select name into :var_t_d separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='D' and name ne 'id'; QUIT; %put var_t_d; proc sql noprint; select name into :var_t_a separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='A' and name ne 'id'; QUIT; %put var_t_a; proc sql noprint; select name into :var_t_b separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='A' and name ne 'id'; QUIT; %put var_t_b; data c; set b; length missvar $30.; array char{*} _character_; do i=1 to dim(char); if missing(char{i}) then call catx(' ',missvar,vname(char{i})); end; drop i; if _n_=1 then call symput('missvar_',missvar); run; %put missvar_; proc sql noprint; select name into :var_t_c separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='C' and name ne 'id'; QUIT; %put var_t_c; data b; set b; drop all_var missvar_; run; proc sql noprint; select name into :var_t_B_ separated by ' ' from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='B' ; QUIT; %put var_t_B_; data a; set a; drop var_t_B_; run; proc sql; create table union_t as select a.*,b.* from a left join b on a.aa_1=input(compress(b.id),8.) ;quit; data last_t; retain var_t_d; set union_t; drop all_var id; run; proc sql noprint; select compress(b.name||'='||a.name) into :renam separated by ' ' from ( select varnum,name from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='E') a left join ( select varnum,name from sashelp.vcolumn where LIBNAME='WORK' AND MEMNAME='D') b on a.varnum=b.varnum ;quit; %put renam; data a; set last_t; rename renam; run; proc sql; drop table a_0,b,c,d,e,last_t,union_t; quit; %mend; %vartonum;
-
个人分类: SAS程序|23 次阅读|1 个评论
-
-
分享
E9安装程序命令行参数(适用于其它程序的安装)
-
mzhua1989 2015-11-13 00:02
-
E9安装程序命令行参数 Internet Explorer 9 有很多命令行参数，可配合安装程序一起使用。 要使用参数，您需要按以下格式在命令提示符 (cmd.exe) 下 运行安装程序(以 Windows 7 为例): 在Windows 7系统中免更新安装IE9 32位安装包“IE9-Windows7-x86-chs.exe”，一般的IE9安装需要检查更新。 帮助选项: 命令行参数 说明 /help 显示此文档。 安装模式: 命令行参数 说明 /passive 运行安装时不需要用户干预。 /quiet 等同于 /passive，但不显示任何用户界面。 安装选项: 命令行参数 说明 /update-no 不查找 Internet Explorer 更新。 /ieak-full:path 保留供 Internet Explorer 管理包使用。 /ieak-branding:path 保留供 Internet Explorer 管理包使用。 /closeprograms 自动关闭程序以便可以在不重新启动系统的情况下安装。 重新启动选项: 命令行参数 说明 /norestart 安装完成后不重新启动。 /forcerestart 如果安装后要求重新启动，则自动重新启动。 其他选项: 命令行参数 说明 /log:path 创建日志文件，位置: path。 /x:path 提取安装包内容，位置: path。 使用教程: 1、先把“IE9-Windows7-x86-chs.exe”放置到D盘目录下。 2、点任务栏左下放的“Windows图标”，在输入框中输入“cmd”，再按“Enter”按钮即可开启cmd.exe。 3、在打开的cmd窗口中输入“D:”，再按“Enter”进入D盘符下。 4、接着输入“IE9-Windows7-x86-chs.exe /update-no”，即可开启IE9免更新安装。 您可以在IE9-Windows7-x86-chs.exe后面添加不同的参数来控制IE9的安装。
-
个人分类: SAS|72 次阅读|0 个评论
-
-
分享
裂变势能曲面—TALYS
-
accumulation 2015-7-4 21:53
-
裂变产额的系统学 TALYS 程序 TALYS 程序是近几年发展的一个大型核反应程序，它包括 直接反应模 型，预平衡反应模型，裂变模型和系统学模型 ，可以用来计算 入射粒子为中子、 质子、氘、氚、3He和α ，靶核质量数大于12 的所有核反应 。与ALICE-91 程序 相比，它的最大的改进是采用将Hauser-Feshbach和Hill-Wheeler相结合的方法来 计算裂变几率，从而使 低能中子诱发裂变的计算 成为可能。将TALYS和温度相 关的Brosa模型相结合原则上可以 计算入射粒子能量在0 到200MeV范围内的裂 变碎片产额和裂变产物产额。
-
个人分类: 裂变模型|0 个评论
-
-
分享
裂变势能曲面
-
accumulation 2015-7-4 16:45
-
本工作独立编写了高维势能曲面计算程序和搜索算法的程序。这两个程序都 提供了不同的理论模型和算法供选择，势能曲面计算程序中可以选择不同宏观模 型和独立粒子势；搜索算法程序可以指定采用模拟降水算法或模拟泛洪算法。程 序的有效性的验证是通过对各种输入和输出进行分析检验来进行的。计算了U 和Pu 元素一系列同位素的裂变位垒，并与其他人的工作了比较，提出了进一步 改进程序的方向。 用Brosa 模型（多模式无规颈断裂模型）计算了低能中子诱发U和Pu元素 裂变的碎片质量分布，动能分布和累计产额。计算裂变碎片的质量分布和动能分 布时，将裂变系统的断前形状的几个主要参量与裂变后现象的物理量和分布相联 系，从而确定裂变核断前形状。计算累计产额时，使用了大型核反应程序Talys 计算了各个反应道的截面。 关键词：裂变位垒，势能曲面，裂变产额，Brosa 模型 （1） 势能曲面计算程序中可以选择 不同宏观模 型和独立粒子势； （2） 搜索算法程序可以指定采用模拟降水算法或模拟泛洪算法； （3） 计算U 和Pu 元素一系列同位素的裂变位垒； （4） 用Brosa 模型（多模式无规颈断裂模型）计算了低能中子诱发U和Pu元素 裂变的碎片质量分布，动能分布和累计产额； （5） 将裂变系统的断前形状的几个主要参量与裂变后现象的物理量和分布相联 系，从而确定裂变核断前形状； （6） 计算累计产额时，使用了大型核反应程序Talys 计算了各个反应道的截面；
-
个人分类: 裂变模型|0 个评论
-
-
分享
裂变势能曲面
-
accumulation 2015-7-4 16:18
-
独立编写了宏观-微观模型的高维势能曲面计算程序。核的表面用广义 Lawrence 形状描述。采用Strutinsky 的宏观势能（液滴模型能）加微观壳修正的 方法计算形变核的势能。计算液滴模型能时采用了两组不同的公式，即 Myers-Swiatecki 公式和LSD 公式。独立粒子势也采用了变形核Woods-Saxon势 和折叠Yukawa 势两种。比较不同理论模型的计算结果，确定了LSD 公式加折 叠汤川势的模型组合能得出最合理的结果。这是第一次将LSD 公式和折叠汤川 势与Lawrence 形状结合来计算裂变截面。 计算U 和Pu 同位素的势能曲面时，使用了多达400 万格点的大型网格，保 证了计算结果包含了裂变路径上所有的关键细节的信息。 （1） 宏观-微观模型的高维势能曲面计算程序； （2） 核的表面用广义 Lawrence 形状描述； （3） 采用Strutinsky 的宏观势能（液滴模型能）加微观壳修正的 方法计算形变核的势能； （4）液滴模型能： Myers-Swiatecki 公式和LSD 公式； （5）独立粒子势： 变形核Woods-Saxon势 和折叠Yukawa 势两种； （6） 比较不同理论模型的计算结果，确定了LSD 公式加折 叠汤川势的模型组合能得出最合理的结果； （7） 第一次将LSD 公式和折叠汤川 势与Lawrence 形状结合来计算裂变截面； （8） 400 万格点的大型网格；
-
个人分类: 裂变模型|0 个评论
-
-
分享
程序块
-
accumulation 2015-7-3 21:56
-
program abe x=2.5 call aa(x,y) call bb(y,z) write(*,1) z 1 format(2x,'z=',E15.7) stop end subroutine aa(x,y ) y=x**2 return end subroutine bb(y,z) z=sin(y) return end
-
个人分类: 裂变模型|0 个评论
-
-
分享
跪求面板数据门槛模型软件及程序，悬赏50论坛币
-
陈荣777 2015-7-1 09:17
-
跪求面板数据门槛模型软件及程序，悬赏50论坛币
-
15 次阅读|0 个评论
-
-
分享
IMSL程序库
-
accumulation 2015-5-16 15:08
-
1.线性系统 2.特征值问题 3.插值和近似 4.数值积分 5.微分方程 6.非线性方程求解 7.优化程序设计
-
个人分类: 裂变模型|0 个评论
-
-
分享
确定初始划界为极小的程序
-
accumulation 2015-5-2 13:25
-
#include cmath #include "nr.h" using namespace std; namespace { inline void shft3(DP a, DP b, DP c, const DP d) { a=b; b=c; c=d; } } void NR::mnbrak(DP ax, DP bx, DP cx, DP fa, DP fb, DP fc, DP func(const DP)) { const DP GOLD=1.618034,GLIMIT=100.0,TINY=1.0e-20; DP ulim,u,r,q,fu; fa=func(ax); fb=func(bx); if (fb fa) { SWAP(ax,bx); SWAP(fb,fa); } cx=bx+GOLD*(bx-ax); fc=func(cx); while (fb fc) { r=(bx-ax)*(fb-fc); q=(bx-cx)*(fb-fa); u=bx-((bx-cx)*q-(bx-ax)*r)/ (2.0*SIGN(MAX(fabs(q-r),TINY),q-r)); ulim=bx+GLIMIT*(cx-bx); if ((bx-u)*(u-cx) 0.0) { fu=func(u); if (fu fc) { ax=bx; bx=u; fa=fb; fb=fu; return; } else if (fu fb) { cx=u; fc=fu; return; } u=cx+GOLD*(cx-bx); fu=func(u); } else if ((cx-u)*(u-ulim) 0.0) { fu=func(u); if (fu fc) { shft3(bx,cx,u,cx+GOLD*(cx-bx)); shft3(fb,fc,fu,func(u)); } } else if ((u-ulim)*(ulim-cx) = 0.0) { u=ulim; fu=func(u); } else { u=cx+GOLD*(cx-bx); fu=func(u); } shft3(ax,bx,cx,u); shft3(fa,fb,fc,fu); } } int main() { }
-
个人分类: 裂变模型|0 个评论
-
-
分享
哪位可以帮忙解释下这个程序，wage是好像不是R自带程序
-
494226841 2015-4-29 00:31
-
eg2. data("wage1") attach(wage1) bw - npudensbw(~lwage+ordered(numdep),tol=.1,ftol=.1,data=wage1) numdep.seq - sort(unique(numdep)) lwage.seq - seq(min(lwage),max(lwage),length=50) wage1.eval - expand.grid(numdep=ordered(numdep.seq),lwage=lwage.seq) fhat - fitted(npudens(bws=bw,newdata=wage1.eval)) f - matrix(fhat,length(unique(numdep)),50) scatterplot3d(wage1.eval ,wage1.eval ,fhat, ylab="Log wage (lwage)", xlab="Number of Dependents (numdep)", zlab="Joint Density", angle=15,box=FALSE,type="h",grid=TRUE,color="blue") detach(wage1)
-
15 次阅读|0 个评论
-
-
分享
金融资产的定价理论与数值计算—附C++程序
-
accumulation 2015-4-18 20:58
-
第1章 货币的时间价值及应用 §1.1 复利计息 1.1.1累积函数 1.1.2 利率 1.1.3单利计息和复利计息 1.1.4 贴现率 1.1.5 复利的终值 1.1.6 计息次数 1.1.7 连续复利 §1.2 多期复利的终值和现值 1.2.1 多期复利的终值 1.2.1 多期复利的现值 1.2.3 年金的终值和现值 §1.3 固定收益证券定价 1.3.1 固定收益证券的基本特征和种类 1.3.2 固定收益证券定价 1.3.3 零息票债券定价 1.3.4 债券的收益率 §1.4 普通股定价 1.4.1 股票定价的基本模型-贴息贴现模型 1.4.2 贴息贴现模型的特殊形式 1.4.3 其它模型 §1.4 本章小结 第2章 远期、期货与互换 §2.1 远期合约定价 2.1.1 无收益证券远期合约 2.1.2 支付已知现金收益证券的远期合约 2.1.3 支付已知红利率证券的远期合约 §2.2 期货合约定价 2.2.1 期货合约与远期合约价格之间的关系 2.2.2 金融期货 2.2.3 利率期货 §2.3 金融互换 2.3.1 利率互换 2.3.2 货币互换 §2.4 本章小结 第3章 投资组合理论 §3.1 投资组合的风险与收益 3.1.1 金融风险定义及种类 3.1.2 证券的风险与收益的度量 3.1.3 证券组合与风险分散 §3.2 标准的均值-方差模型 §3.3 具有无风险资产的均值-方差模型 §3.5 本章小结 第4章 资本市场理论 §4.1 资本资产定价模型 4.1.1 标准资本资产定价模型 4.1.2 价格型资本资产定价模型 §4.2 套利定价模型 4.2.1 因素模型 4.3.2 套利原则 4.3.3 套利组合 4.3.4 套利定价理论模型 4.3.5 APT与CAPM的区别与联系 §4.4 本章小结 第5章 期权定价理论 §5.1 期权概述 5.1.1 期权概念 5.1.2 相关概念 5.1.3 影响期权价格的因素 §5.2 Black-Scholes期权定价理论 5.2.1 Black-Scholes偏微分方程 5.2.2 边界条件 5.2.3 Black-Scholes期权定价公式 5.2.3 看涨与看跌期权之间的评价关系 §5.3 红利的影响 5.3.1 欧式期权 5.3.2 美式期权定价 §5.4 风险对冲 5.3.1 对冲 5.3.2 对冲 5.3.3 对冲 5.3.4 对冲 5.3.5 对冲 §5.4 隐含波动率 5.4.1 二分法 5.4.2 Newton迭代法 §5.5 本章小 结
-
个人分类: 金融工程|0 个评论
-
-
分享
执行力管理格言 119/300 条
-
kychan 2015-4-17 18:18
-
执行力管理格言 119/300 条 人们不会做你希望的事，人们只会做你检查的事！ 你所做的事就是：确保每一个问题的检查程序， 并在你的日历上标明！
-
个人分类: 执行力管理格言 300|37 次阅读|1 个评论
-
-
分享
裂变模型—IMSL程序库
-
accumulation 2015-4-17 00:33
-
DBCONF函数 DBCONF函数
-
个人分类: 裂变模型|0 个评论
-
-
分享
执行力管理格言 102/300 条
-
kychan 2015-4-13 20:08
-
执行力管理格言 102/300 条 在解决问题的程序上， 先从最紧迫、 最重要的事情开始。 一开始就全抓， 反而解决不好。 记住， 重点只有一个！
-
个人分类: 执行力管理格言 300|13 次阅读|1 个评论
-
-
分享
金融工程
-
accumulation 2015-4-7 02:30
-
1.数理金融学与金融数学相关文献； 2.金融物理学等金融工程相关的交叉学科； 3.量化金融与相关程序编程； 4.金融工程、金融计量学相关经典文献。
-
个人分类: 金融工程|0 个评论
GMT+8, 2026-2-13 17:09
标签: 程序经管大学堂：名校名师名课
京公网安备 11010802022788号
论坛法律顾问：王进律师
知识产权保护声明
免责及隐私声明
