应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【分享】自已写的洛仑兹展宽的程序。可以从能级得到ＤＯＳ，及做ＩＲ等
北京石油化工学院2026年研究生招生接收调剂公告
191/1返回列表
查看: 2314  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价12次，作者jjf_sxnu增加金币 11 个
当前主题已经存档。

jjf_sxnu

铜虫 (正式写手)

[资源] 【分享】自已写的洛仑兹展宽的程序。可以从能级得到ＤＯＳ，及做ＩＲ等

/*###############################################################################################################
#
# This program is written for Lorentz boarding by Jianfeng Jia( ShanXi Normail University) . It can
# be used for Freq. NMR. DOS and so on.
# Using it by:        Loz inputfile outputfile delta steplength
# Parameter:        inputfile give the isolated input data want to boarding line by line for
#                    　　 For example:
#                             12.1
#                             13.5
#                             4.
#                             ........
#                      outputfile is the result give in X and Y sets, It can be imported to
#                         Orign to give a plot.
#                           delta is a boarding parameter for Lorentz.
#                      steplength is the step length of X in output file.  
#################################################################################################################*/
#include
#include
#include
#include
#include
#include
#include
#define PI 3.14159265
using namespace std;
main(int argc, char * argv[])
{          if(argc!=5)
        { cout<<"You must and only can provide two file. and then a delta and a steplength"<         exit (1);
        }
        ifstream inf(argv[1]);
        if(inf==NULL)
        {cout<<"Can not open file: "<         exit(1);
        }
        vector inputed_xvalues;
        string s_x;       
        double halfwidth=atof(argv[3]);
        while(!inf.eof())
        {
        inf>>s_x;
        if(!s_x.empty()) inputed_xvalues.push_back(atof(s_x.c_str()));
        }
        sort(inputed_xvalues.begin(),inputed_xvalues.end());
        double max_X=*max_element(inputed_xvalues.begin(),inputed_xvalues.end())+10*halfwidth;
        double min_X=*min_element(inputed_xvalues.begin(),inputed_xvalues.end())-10*halfwidth;
        map output_XX_YY;
        double xx=min_X;
        double step=atof(argv[4]);
        while(xx         {       
                double yy=0;
                vector::iterator iter;
                for(iter=inputed_xvalues.begin();iter!=inputed_xvalues.end();iter++)
                yy=yy+(halfwidth/PI)/((xx-*iter)*(xx-*iter)+halfwidth*halfwidth);
                output_XX_YY.insert(map::value_type(xx,yy));
                xx+=step;
        }
        ofstream of(argv[2]);
        map::iterator m_iter;
        for(m_iter=output_XX_YY.begin();m_iter!=output_XX_YY.end(); m_iter++)
        {       
                of<first<<"\t"<second<         }
}

[ Last edited by zeoliters on 2009-6-13 at 15:37 ]
回复此楼

» 收录本帖的淘帖专辑推荐

双咪唑

» 猜你喜欢
1楼 2009-03-11 21:46:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

erylingjet

铁杆木虫 (著名写手)
不过,怎么用呢?嘿嘿
赞一下(1人) 回复此楼
5楼2009-03-12 18:13:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gwdavid

木虫 (著名写手)

★★★★★ 五星级,优秀推荐

支持原创！
回复此楼
2楼2009-03-11 22:17:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wgpig

木虫 (小有名气)

★★★★★ 五星级,优秀推荐

支持原创！
回复此楼
3楼2009-03-12 08:45:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

★★★★★ 五星级,优秀推荐



不支持,自己都会后悔的
赞一下 回复此楼
4楼2009-03-12 18:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjf_sxnu

铜虫 (正式写手)
引用回帖:
Originally posted by erylingjet at 2009-3-12 18:13:
不过,怎么用呢?嘿嘿

C＋＋编译后产生可执行文件，如Loz

就用  Loz   inputfile  outputfile   delta   step

delta 是展宽系数，step 是生成数据在x轴上的变化值（一般可取0.001）

输入输出都是文本文件，输入文件中只提供一列数据， 如画DOS，就提供每个能级的数据。一个能级一行如：

－20.6
－17.8
－5.6
－0.004
0.03

等，输出文件直接导入Origin画图。
赞一下 回复此楼
6楼2009-03-12 19:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

benjamming

木虫 (正式写手)

★★★ 三星级,支持鼓励

那做红外怎么用呢？
赞一下 回复此楼
7楼2009-03-13 15:39:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)

★★★★★ 五星级,优秀推荐

最近刚好在研究c++，源代码收藏，呵呵。感谢
赞一下 回复此楼
8楼2009-03-13 17:02:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csfn

荣誉版主 (知名作家)

★★★★★ 五星级,优秀推荐

非常强烈的支持
回复此楼
10楼2009-03-27 12:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)

★★★ 三星级,支持鼓励

学习中  谢谢！
回复此楼
11楼2009-03-29 05:38:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

★★★★★ 五星级,优秀推荐

鼓励原创
回复此楼
12楼2009-06-13 22:28:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuzhaonan

铁杆木虫 (著名写手)
学习中谢谢，请楼主讲一下红外用法
赞一下 回复此楼
13楼2009-06-16 08:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lerel12345

木虫 (职业作家)
支持原创！
回复此楼
15楼2009-07-21 17:56:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjf_sxnu

铜虫 (正式写手)
如下是画PDOS，很多朋友问红外的，其实与此类似，看一看，想一想就明白了，仅多乘了一个强度因子。
#include
#include
#include
#include
#include
#include
#include
#define PI 3.14159265
using namespace std;
main(int argc, char * argv[])
{          if(argc!=7)
        {
        exit (1);
        }
        ifstream inf(argv[1]);
        if(inf==NULL)
        //argv[1] is input filename;
        {cout<<"Can not open file: "<         exit(1);
        }
        vector > inputed_line;
        //argv[3] is detal for extention.
        double halfwidth=atof(argv[3]);
        while(!inf.eof())
        {        bool havedata=1;
                vector tmp;
                //argv[5] the number of colum of the inputfile.
                for(int i=0;i                 {        
                        string s_x;
                        inf>>s_x;
                        if(!s_x.empty()) tmp.push_back(atof(s_x.c_str()));
                        else {havedata=0; break;}
                }
                if(havedata) inputed_line.push_back(tmp);
        }
        for(int i=0;i                 {
                        for(int j=0;j.size();j++)
                                cout<[j]<<"\t";
                        cout<                 }
        double max_X=inputed_line[0][0];
        double min_X=inputed_line[0][0];
        for(int i=0;i         {
                if(inputed_line[0]>max_X) max_X=inputed_line[0];
                if(inputed_line[0][0];
        }
        //cout<         max_X+=10*halfwidth;
        min_X-=10*halfwidth;
        vector > out_line;
        double xx=min_X;
        //argv[4] is the step for x;
        double step=atof(argv[4]);
        while(xx         {       
               
                //vector::iterator iter;
                //for(iter=inputed_xvalues.begin();iter!=inputed_xvalues.end();iter++)
                vector tmp;
                tmp.push_back(xx);
                for(int j=1;j                 {
                        double yy=0.0;       
                        for(int i=0;i                         {        double tmp=0.0;
                                tmp=(halfwidth/PI);
        //                        cout<         //                        cout<[0]<         //                        cout<                                
                                double tmp2=(xx-inputed_line[0])*(xx-inputed_line[0])+halfwidth*halfwidth;
        //                        cout<                                 tmp=tmp/tmp2;
        //                        cout<                                 tmp=tmp*inputed_line[j];
        //                        cout<[j]<<"\t"<         //                        tmp=tmp*atoi(argv[6]);
        //                        cout<                                 yy+=tmp;
                        }
        //                cout<                         tmp.push_back(yy);
                }
                out_line.push_back(tmp);
                xx+=step;
        }
        ofstream of(argv[2]);
        for(int i=0;i         {
                for(int j=0;j.size();j++)
                        of<[j]<<"\t";
                of<         }

}
赞一下 回复此楼
16楼2009-08-05 22:06:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mingdong

荣誉版主 (知名作家)

★★★★★ 五星级,优秀推荐

非常谢谢了
回复此楼
17楼2009-10-15 10:14:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

严重支持原创贴
回复此楼
19楼2009-12-07 14:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
fox-fox9楼
2009-03-13 23:18   回复  
 
dalongmao14楼
2009-07-14 16:26   回复  
 
zhangfan219218楼
2009-12-06 21:40   回复  
谢谢分享
相关版块跳转 我要订阅楼主 jjf_sxnu 的主题更新
191/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 考研调剂 +4 美丽的youth_ 2026-04-04 5/250 2026-04-04 17:16 by imissbao
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200 +3 what张 2026-04-04 4/200 2026-04-04 09:50 by rzh123456
[考研] 327求调剂 +9 小卡不卡. 2026-03-29 9/450 2026-04-03 19:48 by 百灵童888
[考研] 生物学308分求调剂（一志愿华东师大） +7 相信必会光芒万� 2026-04-02 7/350 2026-04-03 16:48 by rzh123456
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 调剂求助 +8 想换手机不想解� 2026-04-02 11/550 2026-04-03 14:16 by 星空星月
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 求调剂推荐 +3 南山南@ 2026-04-01 3/150 2026-04-02 12:09 by xiaoranmu
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 一志愿北京科技，085601总分305求调剂 +9 半生瓜！ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 288资源与环境专硕求调剂，不限专业，有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
[考研] 318求调剂 +10 陈晨79 2026-03-30 10/500 2026-03-31 17:37 by 544594351
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
信息提示
关闭
请填处理意见
关闭