| 查看: 2220 | 回復(fù): 37 | ||||
ljw4010榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主
|
[求助]
急:求用Fortran或者其他語言編寫一個程序,處理下面的數(shù)據(jù)。。。。。。。。!
|
|||
|
我現(xiàn)在在處理程序畫圖時很費(fèi)勁,想對獲得數(shù)據(jù)文件處理下,因?yàn)槲覜]學(xué)過Fortran,自學(xué)了一下,但還是沒解決,其他語言寫的也可以,求助大神! 1.我的數(shù)據(jù)文件格式如下: T[K] N DOS(Ef) S s/t R_H kappa0 c chi 10.0000 -0.00088309 0.36829829E+01 -0.22021730E-03 0.26301947E+18 -0.79475299E-03 0.44466758E+11 0.18490799E-02 0.54983779E-10 20.0000 -0.00088309 0.26875931E+01 -0.90876783E-04 0.20008966E+18 -0.19857427E-04 0.68363230E+11 0.32057178E-02 0.40123462E-10 30.0000 -0.00088309 0.22924650E+01 -0.84664015E-04 0.16822569E+18 -0.83733589E-05 0.12349874E+12 0.56812697E-02 0.34224538E-10 40.0000 -0.00088309 0.21550792E+01 -0.98572721E-04 0.15610861E+18 -0.65840150E-05 0.18320424E+12 0.75723680E-02 0.32173486E-10 50.0000 -0.00088309 0.20689471E+01 -0.11499769E-03 0.15044739E+18 -0.61060992E-05 0.24851448E+12 0.91199376E-02 0.30887607E-10 60.0000 -0.00088309 0.19878188E+01 -0.13042206E-03 0.14698777E+18 -0.59478990E-05 0.32102167E+12 0.10595254E-01 0.29676431E-10 70.0000 -0.00088309 0.19072220E+01 -0.14418642E-03 0.14447949E+18 -0.58933740E-05 0.39984025E+12 0.12066150E-01 0.28473190E-10 80.0000 -0.00088309 0.18283895E+01 -0.15633123E-03 0.14246774E+18 -0.58795297E-05 0.48428944E+12 0.13542317E-01 0.27296288E-10 90.0000 -0.00088309 0.17522085E+01 -0.16711278E-03 0.14074025E+18 -0.58843546E-05 0.57466693E+12 0.15037461E-01 0.26158971E-10 10.0000 -0.00883088 0.84410600E+01 -0.11689296E-03 0.24789870E+19 -0.41920906E-06 0.38654460E+12 0.46876655E-02 0.12601779E-09 20.0000 -0.00883088 0.49614612E+01 -0.81871808E-04 0.14192156E+19 -0.65426762E-06 0.65286003E+12 0.87741598E-02 0.74070364E-10 30.0000 -0.00883088 0.46756154E+01 -0.55743996E-04 0.12778854E+19 -0.65291729E-06 0.97812153E+12 0.13236410E-01 0.69802931E-10 40.0000 -0.00883088 0.47012209E+01 -0.50147434E-04 0.12586296E+19 -0.64568125E-06 0.13259755E+13 0.17575579E-01 0.70185200E-10 50.0000 -0.00883088 0.47269535E+01 -0.51492758E-04 0.12561404E+19 -0.64100715E-06 0.16823697E+13 0.22066680E-01 0.70569365E-10 60.0000 -0.00883088 0.47408328E+01 -0.54860269E-04 0.12560297E+19 -0.63779752E-06 0.20366930E+13 0.26804368E-01 0.70776571E-10 70.0000 -0.00883088 0.47488656E+01 -0.58348815E-04 0.12551960E+19 -0.63559342E-06 0.23914729E+13 0.31943350E-01 0.70896493E-10 80.0000 -0.00883088 0.47576572E+01 -0.61496957E-04 0.12534166E+19 -0.63394892E-06 0.27522054E+13 0.37606365E-01 0.71027745E-10 90.0000 -0.00883088 0.47716152E+01 -0.64381860E-04 0.12511220E+19 -0.63256571E-06 0.31236204E+13 0.43859197E-01 0.71236126E-10 10.0000 -0.08830878 0.42141602E+02 0.87953597E-04 0.75126055E+19 -0.98769667E-07 0.23033392E+13 0.41537770E-01 0.62913800E-09 20.0000 -0.08830878 0.43692263E+02 -0.41588947E-05 0.82094250E+19 -0.82516579E-07 0.41555023E+13 0.80748205E-01 0.65228805E-09 30.0000 -0.08830878 0.44663624E+02 -0.21504742E-04 0.82540992E+19 -0.79165706E-07 0.55733828E+13 0.12447749E+00 0.66678964E-09 40.0000 -0.08830878 0.45312345E+02 -0.27007473E-04 0.80589511E+19 -0.79127587E-07 0.75597431E+13 0.17202093E+00 0.67647448E-09 50.0000 -0.08830878 0.45883812E+02 -0.29833271E-04 0.79234737E+19 -0.79037621E-07 0.99651412E+13 0.22031734E+00 0.68500600E-09 60.0000 -0.08830878 0.46378068E+02 -0.31377720E-04 0.78728970E+19 -0.78450172E-07 0.12484013E+14 0.26638841E+00 0.69238480E-09 70.0000 -0.08830878 0.46754191E+02 -0.32326518E-04 0.78761793E+19 -0.77551438E-07 0.14973971E+14 0.30934617E+00 0.69800000E-09 80.0000 -0.08830878 0.46993506E+02 -0.33260061E-04 0.79047459E+19 -0.76560603E-07 0.17417591E+14 0.34962559E+00 0.70157277E-09 90.0000 -0.08830878 0.47101837E+02 -0.34544692E-04 0.79414153E+19 -0.75615674E-07 0.19845552E+14 0.38805026E+00 0.70319006E-09 2.我需要提取第1列和第4列的數(shù)據(jù),并將第四列數(shù)據(jù)乘以10的6次方 3.以中間的空行為分割線,將每個分割后區(qū)域數(shù)據(jù) 依次輸出到 1.txt 、 2.txt、3.txt(數(shù)據(jù)多的話,可能更多塊,更多) 4.分割后的數(shù)據(jù)格式大致這樣: 1.txt為例: T[K] S 10.0000 -0.22021730E-03 20.0000 -0.90876783E-04 30.0000 -0.84664015E-04 40.0000 -0.98572721E-04 50.0000 -0.11499769E-03 60.0000 -0.13042206E-03 70.0000 -0.14418642E-03 80.0000 -0.15633123E-03 90.0000 -0.16711278E-03 5.數(shù)據(jù)附件如下: |

木蟲 (正式寫手)
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
Shell Script + AWK 1.??????????????????(case.trace_fixdoping) 2.??????????1?к??4?е???????????????????????10??6?η? > awk '{print $1,$4}' case.trace_fixdoping ?????1?к??4?е????? awk '{print $1,$4*1000000}' case.trace_fixdoping ?????????????10??6?η? 3.???м????????????????????????????? ????????? 1.txt ?? 2.txt??3.txt??????????????????飬???? > ???? awk '{if($1<800){print $1,$4 >> n".txt"} if($1==800){print $1,$4 >> n".txt"; n=n+1}}' case.trace_fixdoping ???????M?Y?????? T=10 ~ T=800 ????? awk ??y????????Y?????? 800 ????????????n?? ???? question 2 => awk '{if($1<800){print $1,$4*1000000 >> n".txt"} if($1==800){print $1,$4*1000000 >> n".txt"; n=n+1}}' case.trace_fixdoping ??????? awk '{if($1<800){printf ("%s %15.8f\n", $1,$4*1000000) >> n".txt"} if($1==800){print ("%s %15.8f\n", $1,$4*1000000) >> n".txt"; n=n+1}}' case.trace_fixdoping 15.8f ??? ???? 15 λ????С???c???@????? 8 λ ??g???W?????????? ???? 15.8e ????? |
版主 (知名作家)

銀蟲 (小有名氣)
|
可執(zhí)行文件(把你的數(shù)據(jù)文件復(fù)制到dist目錄下,運(yùn)行plot_case.exe即可): http://pan.baidu.com/s/1eQmsNJ0 源代碼: http://pan.baidu.com/s/1bnChxFp 我也是新手,寫的很粗糙,有問題給我留言 |
銀蟲 (小有名氣)
榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主


榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主

榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主

|
T[K] S 10 -220.217300000000 20 -90.8767830000000 30 -84.6640150000000 40 -98.5727210000000 50 -114.997690000000 60 -130.422060000000 70 -144.186420000000 80 -156.331230000000 90 -167.112780000000 T[K] S 10 -116.892960000000 20 -81.8718080000000 30 -55.7439960000000 40 -50.1474340000000 50 -51.4927580000000 60 -54.8602690000000 70 -58.3488150000000 80 -61.4969570000000 90 -64.3818600000000 T[K] S 10 87.9535970000000 20 -4.15889470000000 30 -21.5047420000000 40 -27.0074730000000 50 -29.8332710000000 60 -31.3777200000000 70 -32.3265180000000 80 -33.2600610000000 90 -34.5446920000000 |
榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主

榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主

至尊木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|