版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號(hào): 1530793
注冊(cè): 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

[求助] 用C讀取xyz文件的坐標(biāo)信息

有一個(gè)由pdb轉(zhuǎn)成的xyz文件。
如圖，我想用C將其中的坐標(biāo)信息提取出來，存在一個(gè)新的文件或是數(shù)組當(dāng)中。

7.688102 -2.282595 2.994640
7.303011 -1.469669 3.647740
…………
7.158374 3.922640 -0.572473

只能用C。其實(shí)用awk或其它的很容易實(shí)現(xiàn)。但這部分最后可能會(huì)在一個(gè)大的C當(dāng)中，所以只能用C。

我的思路是每次讀入一行，進(jìn)行判斷。但如何判斷我還沒有想好。因?yàn)槊啃械男畔㈤L(zhǎng)度不同。

現(xiàn)在有其它的工作，比較著急。非常感謝。

要求，程序要有程序說明和相應(yīng)的注釋。

捕獲.GIF

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : diala.xyz

2013-06-26 13:59:40, 1.41 K

» 猜你喜歡

生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有3人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有6人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)
學(xué)碩274求調(diào)劑已經(jīng)有8人回復(fù)
壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有7人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
22408 359分調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

需要將dem轉(zhuǎn)化成xyz坐標(biāo)的文本形式，自己經(jīng)驗(yàn)有限，求流程。已經(jīng)有5人回復(fù)
xyz文件打開方法已經(jīng)有4人回復(fù)
讀取Google截圖坐標(biāo)生成TXT坐標(biāo)文件的程序已經(jīng)有14人回復(fù)
發(fā)布一個(gè)自己寫的色坐標(biāo)計(jì)算軟件CIE1931xy.V1.6 已經(jīng)有15人回復(fù)
【求助】怎么導(dǎo)出xyz坐標(biāo) 已經(jīng)有4人回復(fù)
【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標(biāo)的腳本已經(jīng)有11人回復(fù)

唉。還是學(xué)吧。

1樓 2013-06-26 13:59:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應(yīng)助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時(shí)
蟲號(hào): 1128371
注冊(cè): 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2013-06-27 00:33:56

只要坐標(biāo)就只讀每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"

贊一下

回復(fù)此樓

要夢(mèng)遊，不要催眠。

2樓2013-06-26 15:17:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號(hào): 1530793
注冊(cè): 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by 綠遍山原 at 2013-06-26 15:17:39
只要坐標(biāo)就只讀每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"

綠遍山原:
您好。
1.能把整個(gè)程序?qū)懗鰜韱幔?br /> 2.我查了一下，查到的fgetl都是matlab中的函數(shù)，這是一個(gè)C函數(shù)嗎？
非常感謝。

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

3樓2013-06-26 17:41:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiaochong8693

木蟲之王 (文壇精英)

應(yīng)助: 58 (初中生)
金幣: 74351.8
散金: 3788
紅花: 233
帖子: 47068
在線: 2659.7小時(shí)
蟲號(hào): 752823
注冊(cè): 2009-04-20
專業(yè): 動(dòng)物學(xué)

祝樓主諸事順利，萬事如意

贊一下

回復(fù)此樓

4樓2013-06-26 18:11:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應(yīng)助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時(shí)
蟲號(hào): 1128371
注冊(cè): 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

3樓: Originally posted by meatball1982 at 2013-06-26 17:41:10
綠遍山原:
您好。
1.能把整個(gè)程序?qū)懗鰜韱幔?br /> 2.我查了一下，查到的fgetl都是matlab中的函數(shù)，這是一個(gè)C函數(shù)嗎？
非常感謝。...

是c的函數(shù)

回復(fù)此樓

要夢(mèng)遊，不要催眠。

5樓2013-06-26 22:30:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號(hào): 1119284
注冊(cè): 2010-10-11
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+2, 感謝參與 2013-06-27 00:34:06
meatball1982: 金幣+50, ★★★很有幫助, 程序完整，有說明。嘗試過，好用。 2013-06-27 15:56:29

CODE:

這個(gè)或許可以滿足你的要求

#include <stdio.h>

#include <stdlib.h>

#include <malloc.h>

#include <string.h>

#define BUFSIZE 512

int main(int argc,char *argv[])

{

  FILE *ifp,*ofp;

  char buf[BUFSIZE];

  int number;       // record entry number

  int i;

  double *x,*y,*z;

  if((ifp = fopen("diala.xyz","r+")) == NULL)   // input file handle 

  {

    printf("Cannot open this file!\n");

        exit(1);

  }

  if((ofp = fopen(argv[1],"w+")) == NULL)   // output filename obtained from first commandline parameter

  {

    printf("Cannot write to this file!\n");

        exit(1);

  }

  fgets(buf,BUFSIZE,ifp);   // read record entry number

  number = atoi(buf);

  // memory allocation for three arrays

  x = (double *)malloc(sizeof(double) * number);

  y = (double *)malloc(sizeof(double) * number);

  z = (double *)malloc(sizeof(double) * number);

  i = 0;

  while(fgets(buf,BUFSIZE,ifp) != NULL)  // read info line by line

  {

    sscanf(buf,"%*d%*s%lf%lf%lf",&x[i],&y[i],&z[i]);

        fprintf(ofp,"%lf  %lf  %lf\n",x[i],y[i],z[i]);  // write into output file

        i++;

  }

  fclose(ifp);

  fclose(ofp);

  // release memory 

  free(x);

  free(y);

  free(z);

  return 0;

}

贊一下

回復(fù)此樓

6樓2013-06-26 22:47:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號(hào): 1530793
注冊(cè): 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

6樓: Originally posted by jackie1179 at 2013-06-26 22:47:21
這個(gè)或許可以滿足你的要求
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>

#define BUFSIZE 512

int main(int argc,char *argv[])
{
F ...

稍加修改，正是我想要的。非常感謝！

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

7樓2013-06-27 15:55:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

迷失的心

銅蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 146.9
帖子: 37
在線: 35.5小時(shí)
蟲號(hào): 2526371
注冊(cè): 2013-06-29
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

想知道buf的大小為什么是512呢？別的數(shù)不行嗎？

贊一下

回復(fù)此樓

努力才能做到最好

8樓2013-06-30 10:58:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號(hào): 1119284
注冊(cè): 2010-10-11
專業(yè): 生物信息學(xué)

引用回帖:

8樓: Originally posted by 迷失的心 at 2013-06-30 10:58:20
想知道buf的大小為什么是512呢？別的數(shù)不行嗎？...

buf的大小可以是其他數(shù)，此值大小不固定，但應(yīng)盡量設(shè)的大一點(diǎn)，否則有可能出錯(cuò)。

贊一下

回復(fù)此樓

9樓2013-07-01 14:45:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

迷失的心

銅蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 146.9
帖子: 37
在線: 35.5小時(shí)
蟲號(hào): 2526371
注冊(cè): 2013-06-29
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

9樓: Originally posted by jackie1179 at 2013-07-01 14:45:14
buf的大小可以是其他數(shù)，此值大小不固定，但應(yīng)盡量設(shè)的大一點(diǎn)，否則有可能出錯(cuò)。...

多謝，只是512這個(gè)數(shù)太常見了....

贊一下

回復(fù)此樓

努力才能做到最好

10樓2013-07-03 22:36:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 meatball1982 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 283求調(diào)劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +5	芝士雪baoo 2026-03-28	6/300	2026-03-28 15:25 by 1018329917
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +4	零八# 2026-03-28	4/200	2026-03-28 15:01 by 17865157980
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 086502化學(xué)工程342求調(diào)劑 +6	阿姨復(fù)古不過 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +9	vv迷 2026-03-22	9/450	2026-03-27 15:59 by 不吃魚的貓
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 329求調(diào)劑 +5	1() 2026-03-22	5/250	2026-03-26 20:40 by fmesaito
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17	RichLi_ 2026-03-25	17/850	2026-03-26 19:44 by plmuchong
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG
[考研] 332求調(diào)劑 +6	032500 2026-03-25	6/300	2026-03-25 22:45 by 418490947
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7	niko- 2026-03-22	7/350	2026-03-25 20:14 by qingfeng258
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 361求調(diào)劑 +3	Glack 2026-03-22	3/150	2026-03-23 22:03 by fuyu_
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin

24小時(shí)熱門版塊排行榜

[求助] 用C讀取xyz文件的坐標(biāo)信息

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: