【正文】 al_codebook(t,m)。 n=n+1。 end fprintf(fid,39。\n39。)。endt=(1:length(z))/fs。%時(shí)間的向量subplot(212),plot(t,z)。xlabel(39。時(shí)間39。),ylabel(39。幅度39。),title(39。碼書39。)。%余下fclose(fid)。cpu_time=toc。disp(39。cpu_time=39。)。disp(cpu_time)。disp(39。pj_distort 平均失真=39。)。disp(pj_distort)。disp(39。xd_shizhen 相對(duì)失真=39。)。disp(xd_shizhen)。disp(39。N1 每個(gè)碼書的樣本個(gè)數(shù)=39。)。disp(N1)。disp(39。j2循環(huán)次數(shù)=39。)。disp(j2)。disp(39。final_codebook=39。)。disp(final_codebook)。覆蓋算法的仿真程序clear allclc%[y,fs]=wavread(39。E:/wav/39。)。y=y(:,1)。save yleng=length(y)。fid=fopen(39。39。,39。wt39。)。for i=1:leng fprintf(fid,39。%d\n39。,y(i))。end。fclose(fid)。codebook_dimen=5。%碼書維數(shù)codebook_size=16。signal_num=200。%參加訓(xùn)練樣本的個(gè)數(shù)fid=fopen(39。39。,39。rt39。)。a=fscanf(fid,39。%f39。)。%把輸入文件中的數(shù)據(jù)賦給afclose(fid)。x=a(1001:1000+signal_num*codebook_dimen)。train_signal=zeros(signal_num,codebook_dimen)。r=1。for i=1:signal_num for j=1:codebook_dimen train_signal(i,j)=x(r)。 r=r+1。 endend。R=。y=train_signal。y_num=signal_num。k=1。disp(39。y=39。)。disp(y)。baoqiang_num=1000。t=0。tic。while (baoqiang_numcodebook_size)y=train_signal。 y_num=signal_num。k=1。R=R+。disp(39。R=39。)。disp(R)。i=0。mazi=zeros(signal_num,codebook_dimen)。num=zeros(1,signal_num)。while (y_numcodebook_dimen)mazi(k,:)=y(1,:)。 k=k+1。 p=0。 q=0。 for j=2:y_num A=0。 for m=1:codebook_dimen A=A+(y(j,m)y(1,m))^2。 end。 if AR p=p+1。 else q=q+1。 end。end。baoqiang=zeros(p,codebook_dimen)。shengyu=zeros(q,codebook_dimen)。p=1。q=1。for j=2:y_num A=0。 for m=1:codebook_dimen A=A+(y(j,m)y(1,m))^2。 end。 if AR baoqiang(p,:)=y(j,:)。 p=p+1。 else shengyu (q,:)=y(j,:)。 q=q+1。 end。end。disp(39。baoqiang=39。)。disp(baoqiang)。disp(39。shengyu=39。)。disp(shengyu)。y=shengyu。y_num=length(y)。disp(39。y_num=39。)。disp(y_num)。disp(39。小循環(huán)39。)i=i+1。disp(39。i=39。)。disp(i)。num(1,i)=p。end。disp(39。mazi=39。)。disp(mazi)。t=t+1。disp(39。t=39。)。disp(t)。baoqiang_num=i。disp(39。baoqiang_num=39。)。disp(baoqiang_num)。disp(39。 大循環(huán)* * * * * * * * 39。)end。disp(39。劃分胞腔次數(shù)=39。)。disp(t)。disp(39。最終碼書中碼字個(gè)數(shù)=39。)。disp(baoqiang_num)。disp(39。最終失真半徑=39。)。disp(R)。for c=1:igeshu(c,1)=num(1,c)。 c=c+1。end。disp(39。每個(gè)胞腔中的向量數(shù)=39。)。disp(geshu)。disp(39。cpu_time=39。)。disp(toc)。n=geshu。a=0。b=0。long=length(n)。for i=1:long a=a+n(i,1)。 i=i+1。end。a=a/long。for i=1:long b=b+(n(i,1)a)^2。 i=i+1。end。b=b/long。disp(39。方差為39。)。disp(b)。成對(duì)最鄰近（PNN）算法的仿真程序clear allclc%[y,fs]=wavread(39。E:/wav/39。)。y=y(:,1)。save yleng=length(y)。fid=fopen(39。39。,39。wt39。)。for i=1:leng fprintf(fid,39。%d\n39。,y(i))。end。fclose(fid)。codebook_dimen=5。%碼書維數(shù)signal_num=200。%參加訓(xùn)練樣本的個(gè)數(shù)fid=fopen(39。39。,39。rt39。)。a=fscanf(fid,39。%f39。)。%把輸入文件中的數(shù)據(jù)賦給afclose(fid)。x=a(1001:1000+signal_num*codebook_dimen)。train_signal=zeros(signal_num,codebook_dimen)。r=1。for i=1:signal_num for j=1:codebook_dimen train_signal(i,j)=x(r)。 r=r+1。 endend。y=train_signal。%怎樣把一個(gè)矩陣賦給另一個(gè)矩陣???y_num=signal_num。disp(39。y=39。)。disp(y)。tica=y。disp(39。a=39。)。disp(a)。codebook_size=16。geshu=length(a)。mazi_num=ones(geshu,1)。weishu=5。x=zeros(1,weishu)。while(geshucodebook_size)%如果胞腔大于六b=zeros(geshu)。t=geshu1。for i=1:t for j=1:geshu if ij b(i,j)=sum((a(i,:)a(j,:)).^2)。 else b(i,j)=0。 end。 end。end。zuida=max(max(b))。%找出矩陣中最大值disp(39。b1=39。)disp(b)。%把矩陣中零元素全置成最大值for i=1:geshu for j=1:geshu if i=j b(i,j)=zuida。 end。 end。end。disp(39。b2=39。)disp(b)。zuixiao=min(min(b))。for i=1:geshu for j=1:geshu if b(i,j)==zuixiao m=i。 n=j。 end end。end。disp(39。m=39。)。disp(m)。disp(39。n=39。)。disp(n)。%可以找出矩陣最小值元素，并把行數(shù)列數(shù)賦給m,ndisp(39。碼字間失真最小值=39。)。disp(zuixiao)。mazi_num(m,1)=mazi_num(m,1)+mazi_num(n,1)。 x(1,:)=a(m,:)+a(n,:)。 x(1,:)=x(1,:)/2。 a(m,:)=x(1,:)。 disp(39。x=39。)。 disp(x)。%過渡用的一行矩陣 disp(39。a1=39。)。disp(a)。%兩個(gè)問題，計(jì)算新碼字賦給yi，算每個(gè)胞腔中的碼字個(gè)數(shù)if n==geshu a(n,:)=0。 mazi_num(n,1)=0。else a(n,:)=a(geshu,:)。 a(geshu,:)=0。 mazi_num(n,1)=mazi_num(geshu,1)。 mazi_num(geshu,1)=0。end。disp(39。a2=39。)。disp(a)。disp(39。mazi_num39。)。disp(mazi_num)。geshu=geshu1。disp(39。39。)。end。for i=1:codebook_size final_book(i,:)=a(i,:)。end。k=signal_numcodebook_size。disp(39。循環(huán)次數(shù)=39。)。disp(k)。disp(39。final_book=39。)。disp(final_book)。disp(39。cup_time=39。)。disp(toc)。disp(39。每個(gè)胞腔中的碼字個(gè)數(shù)39。)。for i=1:codebook_size d(i,1)=mazi_num(i,1)。end。disp(d)。n=d。a=0。b=0。long=length(n)。for i=1:long a=a+n(i,1)。 i=i+1。end。a=a/long。for i=1:long b=b+(n(i,1)a)^2。 i=i+1。end。b=b/long。disp(39。方差為39。)。disp(b)。謝 辭在畢業(yè)設(shè)計(jì)完成之際，我真心的感謝在設(shè)計(jì)之中給予我?guī)椭臈钣裉m老師和我的同班同學(xué)，在這三個(gè)月中，楊老師給予了鼓勵(lì)和耐心的指導(dǎo)，特別是在課題的理解上，楊老師嚴(yán)謹(jǐn)?shù)目茖W(xué)思維和專業(yè)的知識(shí)對(duì)我產(chǎn)生了很大的影響，并使我受益非淺，當(dāng)我在編譯程序遇到瓶頸時(shí)，楊老師幫助我檢查程序并不斷的鼓勵(lì)我讓我不斷探索。同班同學(xué)王振宇也給我很多的幫助，使我做到對(duì)MATLAB編程語言從入門到熟練，同時(shí)還有互聯(lián)網(wǎng)上MATLAB論壇上的很多熱心網(wǎng)友的幫助，在這里我衷心的感謝他們。畢業(yè)論文（設(shè)計(jì)）原創(chuàng)性聲明本人所呈交的畢業(yè)論文（設(shè)計(jì)）是我在導(dǎo)師的指導(dǎo)下進(jìn)行的研究工作及取得的研究成果。據(jù)我所知，除文中已經(jīng)注明引用的內(nèi)容外，本論文（設(shè)計(jì)）不包含其他個(gè)人已經(jīng)發(fā)表或撰寫過的研究成果。對(duì)本論文（設(shè)計(jì)）的研究做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中作了明確說明并表示謝意。 作者簽名： 日期： 畢業(yè)論文（設(shè)計(jì)）授權(quán)使用說明本論文（設(shè)計(jì)）作者完全了解**學(xué)院有關(guān)保留、使用畢業(yè)論文（設(shè)計(jì)）的規(guī)定，學(xué)校有權(quán)保留論文（設(shè)計(jì)）并向相關(guān)部門送交論文（設(shè)計(jì)）的電子版和紙質(zhì)版。有權(quán)將論文（設(shè)計(jì)）用于非贏利目的的少量復(fù)制并允許論文（設(shè)計(jì)）進(jìn)入學(xué)校圖書館被查閱。學(xué)?？梢怨颊撐模ㄔO(shè)計(jì)）的全部或部分內(nèi)容。保密的論文（設(shè)計(jì)）在解密后適用本規(guī)定。 作者簽名： 指導(dǎo)教師