【发布时间】:2011-04-26 03:44:25
【问题描述】:
我正在尝试编写一个程序,它将一个文件复制到一个包含 DNA 碱基的字符串中(到目前为止还不错)。然后将碱基转换为蛋白质,首先找到第一个 ATG 序列,然后读取 3 个序列并转换它们,将它们写入另一个文件。
到目前为止,程序在进入第一个 for 循环之前就崩溃了。我不知道是什么导致了问题。
int proteina(char DNA_origem[], char proteina_destino[]){
char aminocidosING [64][14]={"Isoleucine","Isoleucine","Isoleucine","Leucine","Leucine","Leucine","Leucine","Leucine","Leucine","Valine","Valine","Valine","Valine","Phenylalanine","Phenylalanine","Methionine","Cysteine","Cysteine","Alanine","Alanine","Alanine","Alanine","Glycine","Glycine","Glycine","Glycine","Proline","Proline","Proline","Proline","Threonine","Threonine","Threonine","Threonine","Serine","Serine","Serine","Serine","Serine","Serine","Tyrosine","Tyrosine","Tryptophan","Glutamine","Glutamine","Asparagine","Asparagine","Histidine","Histidine","Glutamic acid","Glutamic acid","Aspartic acid","Aspartic acid","Lysine","Lysine","Arginine","Arginine","Arginine","Arginine","Arginine","Arginine","Stop codons","Stop codons","Stop codons"};
char aminocidosPT [64][18]={"Isoleucina","Isoleucina","Isoleucina","Leucina","Leucina","Leucina","Leucina","Leucina","Leucina","Valina","Valina","Valina","Valina","Fenilalanina","Fenilalanina","Metionina","Cisteína","Cisteína","Alanina","Alanina","Alanina","Alanina","Glicina","Glicina","Glicina","Glicina","Prolina","Prolina","Prolina","Prolina","Treonina","Treonina","Treonina","Treonina","Serina","Serina","Serina","Serina","Serina","Serina","Tirosina","Tirosina","Triptofano","Glutamina*","Glutamina","Asparagina","Asparagina","Histidina","Histidina","Ácido glutâmico","Ácido glutâmico","Ácido aspártico","Ácido aspártico","Lisina","Lisina","Arginina","Arginina","Arginina","Arginina","Arginina","Arginina","Códons Stop","Códons Stop","Códons Stop"};
char codoes[64][3]={"ATT","ATC","ATA","CTT","CTC","CTA","CTG","TTA","TTG","GTT","GTC","GTA","GTG","TTT","TTC","ATG","TGT","TGC","GCT","GCC","GCA","GCG","GGT","GGC","GGA","GGG","CCT","CCC","CCA","CCG","ACT","ACC","ACA","ACG","TCT","TCC","TCA","TCG","AGT","AGC","TAT","TAC","TGG","CAA","CAG","AAT","AAC","CAT","CAC","GAA","GAG","GAT","GAC","AAA","AAG","CGT","CGC","CGA","CGG","AGA","AGG","TAA","TAG","TGA"};
char proteinas[64][1] = {"I","I","I","L","L","L","L","L","L","V","V","V","V","F","F","M","C","C","A","A","A","A","G","G","G","G","P","P","P","P","T","T","T","T","S","S","S","S","S","S","Y","Y","W","Q","Q","N","N","H","H","E","E","D","D","K","K","R","R","R","R","R","R",".",".","."};
/* a esta altura suponho que tenhas definido na main as strings dos aminoácidos*/
char **string1;
FILE * ficheiro;
FILE * ficheiro_close;
int f_cmp;
int k, i, start=0; /* variavel de comprimento */
char proteina_origem;
ficheiro = fopen(DNA_origem,"r"); /* DNA origem e a variavel onde ta guardada o nome do ficheiro do utilizador */
ficheiro_close = fopen(proteina_destino,"w+");
fscanf(ficheiro,"%c",string1); /* isto lê os conteudos da stream para a string, copiando pra lá o ficheiro. */
for(i=1;i<=f_cmp;i++) {
if (strncmp(string1[i],codoes[15],3)==0) {
fputs(proteinas[15],ficheiro_close);
for(k=i+2;k<=f_cmp;k+3) {
if ((strncmp(string1[k],codoes[k],3))==0) {
fputs(proteinas[k],ficheiro_close);
if (k==61&&k==62&&k==63) {
return(0);
}
}
}
}
}
}
此外,如果我在 char 定义中不使用 **,编译器会发出警告。你能解释一下吗?这应该只是一个简单的项目,但我被困在最后一个功能中..
不要介意 var 名称和 cmets,它是葡萄牙语。
非常感谢您的宝贵时间!
【问题讨论】:
-
在使用之前确保
ficheiro是一个有效的文件指针(通过检查if(ficheiro) -
我之前有一个 printf 来测试,打印出变量,它就像一个魅力。
-
还有一点代码缺失,但还是不行,以防万一,在for循环之前
f_cmp=strlen(*string1); -
好的,那可能不是这样,但你还是应该检查你的文件指针:)
-
String1 在很多方面都困扰着我,因为首先它被用作 fscanf 中的一个字符。然后你尝试通过取消引用它来像字符串一样使用它,但它不是一个字符串......并且你跳过了第 0 个成员(数组从 C 中的 0 到 size-1),所以你可能会超过结尾。跨度>