




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、 部分程序運行結果:1.求A類DNA的中心坐標的運行結果,如下所示:aggcacggaaaaacgggaataacggaggaggacttggcacggcattacacggaggacgaggtaaaggaggcttgtctacggccggaagtgaagggggatatgaccgcttgg111nt=0.135135ng=0.396396nc=0.171171na=0.297297cggaggacaaacgggatggcggtattggaggtggcggactgttcggggaattattcggtttaaacgggacaaggaaggcggctggaacaaccggacggtggcagcaaa
2、gga111nt=0.153153ng=0.414414nc=0.162162na=0.27027gggacggatacggattctggccacggacggaaaggaggacacggcggacatacacggcggcaacggacggaacggaggaaggagggcggcaatcggtacggaggcggcgga111nt=0.0630631ng=0.45045nc=0.216216na=0.27027atggataacggaaacaaaccagacaaacttcggtagaaatacagaagcttagatgcatatgttttttaaataaaatttgtattattatggtatc
3、ataaaaaaaggttgcga111nt=0.288288ng=0.18018nc=0.108108na=0.423423cggctggcggacaacggactggcggattccaaaaacggaggaggcggacggaggctacaccaccgtttcggcggaaaggcggagggctggcaggaggctcattacggggag111nt=0.108108ng=0.423423nc=0.234234na=0.234234atggaaaattttcggaaaggcggcaggcaggaggcaaaggcggaaaggaaggaaacggcggatatttcggaagtggata
4、ttaggagggcggaataaaggaacggcggcaca111nt=0.126126ng=0.396396nc=0.126126na=0.351351atgggattattgaatggcggaggaagatccggaataaaatatggcggaaagaacttgttttcggaaatggaaaaaggactaggaatcggcggcaggaaggatatggaggcg111nt=0.189189ng=0.36036nc=0.0990991na=0.351351atggccgatcggcttaggctggaaggaacaaataggcggaattaaggaaggcgttctcgcttt
5、tcgacaaggaggcggaccataggaggcggattaggaacggttatgagg111nt=0.189189ng=0.369369nc=0.162162na=0.279279atggcggaaaaaggaaatgtttggcatcggcgggctccggcaactggaggttcggccatggaggcgaaaatcgtgggcggcggcagcgctggccggagtttgaggagcgcg111nt=0.153153ng=0.432432nc=0.207207na=0.207207tggccgcggaggggcccgtcgggcgcggatttctacaagggcttcct
6、gttaaggaggtggcatccaggcgtcgcacgctcggcgcggcaggaggcacgcgggaaaaaacg110nt=0.136364ng=0.409091nc=0.272727na=0.1818180.135135A類坐標中心0.1541770.3832510.1759212.求B類DNA的中心坐標的運行結果,如下所示:gttagatttaacgttttttatggaatttatggaattataaatttaaaaatttatattttttaggtaagtaatccaacgtttttattactttttaaaattaaatatttatt110na=0.5 nt=0.1 n
7、g=0.0454545gtttaattactttatcatttaatttaggttttaattttaaatttaatttaggtaagatgaatttggttttttttaaggtagttatttaattatcgttaaggaaagttaaa110na=0.5 nt=0.145455 ng=0.0272727gtattacaggcagaccttatttaggttattattattatttggattttttttttttttttttttaagttaaccgaattattttctttaaagacgttacttaatgtcaatgc110na=0.518182 nt=0.127273 ng=0.1gt
8、tagtcttttttagattaaattattagattatgcagtttttttacataagaaaatttttttttcggagttcatattctaatctgtctttattaaatcttagagatatta110na=0.5 nt=0.118182 ng=0.0818182gtattatatttttttatttttattattttagaatataatttgaggtatgtgtttaaaaaaaatttttttttttttttttttttttttttttttaaaatttataaatttaa110na=0.645455 nt=0.0636364 ng=0gttatttttaaattta
9、attttaattttaaaatacaaaatttttactttctaaaattggtctctggatcgataatgtaaacttattgaatctatagaattacattattgat110na=0.463636 nt=0.0909091 ng=0.0818182gtatgtctatttcacggaagaatgcaccactatatgatttgaaattatctatggctaaaaaccctcagtaaaatcaatccctaaacccttaaaaaacggcggcctatccc110na=0.263636 nt=0.136364 ng=0.245455gttaattatttattcctt
10、acgggcaattaattatttattacggttttatttacaattttttttttttgtcctatagagaaattacttacaaaacgttattttacatactt110na=0.5 nt=0.0909091 ng=0.118182gttacattatttattattatccgttatcgataattttttacctcttttttcgctgagtttttattcttactttttttcttctttatataggatctcatttaatatcttaa110na=0.563636 nt=0.0727273 ng=0.145455gtatttaactctctttacttttttt
11、ttcactctctacattttcatcttctaaaactgtttgatttaaacttttgtttctttaaggattttttttacttatcctctgttat110na=0.563636 nt=0.0636364 ng=0.172727XB類坐標中心x=0.501818y=0.100909z=0.101818Press any key to continue3.利用歐氏距離判斷第21個DNA所屬類別:tttagctcagtccagctagctagtttacaatttcgacaccagtttcgcaccatcttaaatttcgatccgtaccgtaatttagcttagatttg
12、gatttaaaggatttagattga113t=0.362832g=0.168142c=0.194690.09016120.0324626屬于B類4.利用相關度判斷第21個字符串的類別的運行結果:tttagctcagtccagctagctagtttacaatttcgacaccagtttcgcaccatcttaaatttcgatccgtaccgtaatttagcttagatttggatttaaaggatttagattga與A類的相關度:aggcacggaaaaacgggaataacggaggaggacttggcacggcattacacggaggacgaggtaaaggaggcttgtctac
13、ggccggaagtgaagggggatatgaccgcttgg兩字符串的相關度為:34cggaggacaaacgggatggcggtattggaggtggcggactgttcggggaattattcggtttaaacgggacaaggaaggcggctggaacaaccggacggtggcagcaaagga兩字符串的相關度為:34gggacggatacggattctggccacggacggaaaggaggacacggcggacatacacggcggcaacggacggaacggaggaaggagggcggcaatcggtacggaggcggcgga兩字符串的相關度為:34atggataacg
14、gaaacaaaccagacaaacttcggtagaaatacagaagcttagatgcatatgttttttaaataaaatttgtattattatggtatcataaaaaaaggttgcga兩字符串的相關度為:38cggctggcggacaacggactggcggattccaaaaacggaggaggcggacggaggctacaccaccgtttcggcggaaaggcggagggctggcaggaggctcattacggggag兩字符串的相關度為:38atggaaaattttcggaaaggcggcaggcaggaggcaaaggcggaaaggaaggaaacggcggata
15、tttcggaagtggatattaggagggcggaataaaggaacggcggcaca兩字符串的相關度為:38atgggattattgaatggcggaggaagatccggaataaaatatggcggaaagaacttgttttcggaaatggaaaaaggactaggaatcggcggcaggaaggatatggaggcg兩字符串的相關度為:38atggccgatcggcttaggctggaaggaacaaataggcggaattaaggaaggcgttctcgcttttcgacaaggaggcggaccataggaggcggattaggaacggttatgagg兩字符串的相關
16、度為:38atggcggaaaaaggaaatgtttggcatcggcgggctccggcaactggaggttcggccatggaggcgaaaatcgtgggcggcggcagcgctggccggagtttgaggagcgcg兩字符串的相關度為:38tggccgcggaggggcccgtcgggcgcggatttctacaagggcttcctgttaaggaggtggcatccaggcgtcgcacgctcggcgcggcaggaggcacgcgggaaaaaacg兩字符串的相關度為:38x平均相關度:36.8與B類的相關度:gttagatttaacgttttttatggaatttatg
17、gaattataaatttaaaaatttatattttttaggtaagtaatccaacgtttttattactttttaaaattaaatatttatt兩字符串的相關度為:40gtttaattactttatcatttaatttaggttttaattttaaatttaatttaggtaagatgaatttggttttttttaaggtagttatttaattatcgttaaggaaagttaaa兩字符串的相關度為:40gtattacaggcagaccttatttaggttattattattatttggattttttttttttttttttttaagttaaccgaattattttcttta
18、aagacgttacttaatgtcaatgc兩字符串的相關度為:40gttagtcttttttagattaaattattagattatgcagtttttttacataagaaaatttttttttcggagttcatattctaatctgtctttattaaatcttagagatatta兩字符串的相關度為:40gtattatatttttttatttttattattttagaatataatttgaggtatgtgtttaaaaaaaatttttttttttttttttttttttttttttttaaaatttataaatttaa兩字符串的相關度為:44gttatttttaaatttaatttt
19、aattttaaaatacaaaatttttactttctaaaattggtctctggatcgataatgtaaacttattgaatctatagaattacattattgat兩字符串的相關度為:44gtatgtctatttcacggaagaatgcaccactatatgatttgaaattatctatggctaaaaaccctcagtaaaatcaatccctaaacccttaaaaaacggcggcctatccc兩字符串的相關度為:44gttaattatttattccttacgggcaattaattatttattacggttttatttacaattttttttttttgtcctataga
20、gaaattacttacaaaacgttattttacatactt兩字符串的相關度為:44gttacattatttattattatccgttatcgataattttttacctcttttttcgctgagtttttattcttactttttttcttctttatataggatctcatttaatatcttaa兩字符串的相關度為:44gtatttaactctctttactttttttttcactctctacattttcatcttctaaaactgtttgatttaaacttttgtttctttaaggattttttttacttatcctctgttat兩字符串的相關度為:44x平均相關度:42.4
21、屬于B類Press any key to continue其余字符串的判定同理。附程序:程序1.求兩類DNA 的坐標中心:#include "iostream.h"/求兩類DNA的中心坐標的函數#include "string.h"#define N 10void main() char ch1150; float aN,bN,cN; double f=0,d=0,e=0; int k=0; while(cin>>ch1&&k<N) int na=0,ng=0,nt=0,nc=0; int l=strlen(ch1);
22、cout<<l<<endl; for(int i=0;i<l;i+) if(ch1i='a') na+; else if(ch1i='t') nt+; else if(ch1i='g') ng+; else if(ch1i='c') nc+; cout<<"na="<<(float)nt/l; cout<<" nt="<<(float)ng/l; cout<<" ng="<&l
23、t;(float)nc/l<<endl; ak=(float)nt/l; bk=(float)ng/l; ck=(float)nc/l; k+; for(int j=0;j<N;j+) f+=aj; d+=bj; e+=cj; f=f/N; d=d/N; e=e/N; cout<<"x="<<f<<endl;/nt cout<<"y="<<d<<endl;/ng cout<<"z="<<e<<endl;/nc程
24、序2 .通過歐氏距離判斷DNA類別:#include "iostream.h"#include "string.h"void main()double a3,b3,c3;a0=0.154177;/nta1=0.383251;/nga2=0.175921;/ncb0=0.501818;/ntb1=0.100909;/ngb2=0.101818;/ncfloat d1,d2;char ch1150,ch2150; while(cin>>ch1) int na=0,ng=0,nt=0,nc=0; int l=strlen(ch1); cout<
25、;<l<<endl; for(int i=0;i<l;i+) if(ch1i='a') na+; else if(ch1i='t') nt+; else if(ch1i='g') ng+; else if(ch1i='c') nc+; c0=(float)nt/l; c1=(float)ng/l; c2=(float)nc/l; cout<<"t="<<c0<<endl; cout<<"g="<<c1<
26、<endl; cout<<"c="<<c2<<endl; d1=(c0-a0)*(c0-a0)+(c1-a1)*(c1-a1)+(c2-a2)*(c2-a2); cout<<d1<<endl; d2=(c0-b0)*(c0-b0)+(c1-b1)*(c1-b1)+(c2-b2)*(c2-b2); cout<<d2<<endl; if(d1>d2)/距離B類中心更近些 cout<<"屬于B類"<<endl; else if(d1<d
27、2) cout<<"屬于A類"<<endl; else if(d1=d2) cout<<"無法判別類別"<<endl; 程序3. 通過計算相關度來確定DNA 所屬類別。#include "iostream.h"/相關度計算#include "string.h"#define N1 120#define N2 120#define N3 120#define N 10int min(int a,int b)if(a>b)return b;elsereturn a;void tell(char ch1,float *x)/最后要多輸入一個不相關的字符串結束,
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
評論
0/150
提交評論