matlab實(shí)習(xí)報(bào)告

1、實(shí)驗(yàn)報(bào)告題目一 : a sinusoid of frequency 0=0.1and duration of 300 samples, that is ,0n300, is input to a (causal) filter with transfer function ,where a=0.97。adjust the scale factor b such that the filters gain at 0 is unity。 determine and plot the input x(n) and the output y(n) of the filter over the inte

2、rval 0n450, by iterating the difference equation of the filter. at the same time, plot the filters magnitude response.1 實(shí)驗(yàn)?zāi)康慕o定一個(gè)系統(tǒng)的傳輸函數(shù)，輸入一個(gè)正弦序列，經(jīng)過系統(tǒng)的傳輸函數(shù)后求其差分方程，同時(shí)考慮零點(diǎn)，極點(diǎn)對(duì)幅頻特性曲線的影響。2 實(shí)驗(yàn)原理一般時(shí)域離散系統(tǒng)或網(wǎng)絡(luò)可以用差分方程，單位脈沖響應(yīng)以及系統(tǒng)函數(shù)進(jìn)行描述。3 實(shí)驗(yàn)步驟簡(jiǎn)述求系數(shù)b的值，利用其模值求出b的值。求出其差分方程y(n)，利用濾波函數(shù)將其傳輸函數(shù)進(jìn)行逆變換。求出傳輸函數(shù)的幅頻特性曲線，畫出輸入正弦

3、函數(shù)的圖形，畫出差分方程的圖形開始利用其模值求出b的值調(diào)用filter函數(shù)求差分方程調(diào)用繪圖程序繪制正弦圖形調(diào)用繪圖程序繪制差分方程結(jié)束4 實(shí)驗(yàn)程序框圖 5 實(shí)驗(yàn)源程序 當(dāng)a=0.97時(shí)close all;clear all;clctica1=0.97;w=0.1*pi;b=abs(1-a1*exp(-j*w);n=300;n=1:n;xn=sin(w*n);x1=xn zeros(1,150);a=1 -a1;yn=filter(b,a,x1);figure(name,實(shí)驗(yàn)一,color,r);subplot(211);plot(xn);title( x(n);xlabel(n);ylabe

4、l(幅度);t1=tocsubplot(212);plot(yn);title( y(n);xlabel(n);ylabel(幅度);figure(name,實(shí)驗(yàn)一,color,r);freqz(b,a);pauseclose all;6.實(shí)驗(yàn)結(jié)果當(dāng)a=0.35時(shí)close all;clear all;clctica1=0.35;w=0.1*pi;b=abs(1-a1*exp(-j*w);n=300;n=1:n;xn=sin(w*n);x1=xn zeros(1,150);a=1 -a1;yn=filter(b,a,x1);figure(name,實(shí)驗(yàn)一,color,r);subplot(21

5、1);plot(xn);title( x(n);xlabel(n);ylabel(幅度);t1=tocsubplot(212);plot(yn);title( y(n);xlabel(n);ylabel(幅度);figure(name,實(shí)驗(yàn)一,color,r);freqz(b,a);pauseclose all;實(shí)驗(yàn)結(jié)果7. 實(shí)驗(yàn)總結(jié): 當(dāng)a值變小后，它的相頻響應(yīng)的曲線會(huì)變的平緩題目 （二）一：實(shí)驗(yàn)題目： it is desired to extract signal x(n)=s(n)+v(n)=s+v(n),where v(n) is zero-mean white caussian no

6、ise of variance .(1)first-order iir smoother; to this end,the following iir lowpass filter is used:h(z)=b/(1-az-1),where the parameter a is restrict to the range 0alimit,break,endendfor r=1:4; if val(r)limit,break,endenddisp(touch-tone symbol=,setstr(tm(r,s-4)子程序?yàn)閒unction xf=gfft(x,n,k)if length(x)n

7、 xe=x zeros(1,n-length(x);else xe=x;endx1=xe 0;d1=2*cos(2*pi*k/n);w=exp(-i*2*pi*k/n);y=filter(1,1 -d1 1,x1);xf=y(n+1)-w*y(n);七 程序結(jié)果及圖表：輸入a后，分別輸入需要檢測(cè)的點(diǎn)，會(huì)有相應(yīng)的圖形。八 實(shí)驗(yàn)總結(jié)：本實(shí)驗(yàn)是采用dft而不采用fft進(jìn)行信號(hào)檢測(cè)，這是因?yàn)閷?shí)驗(yàn)中只要求對(duì)整數(shù)點(diǎn)個(gè)個(gè)別值進(jìn)行檢測(cè)，如果用fft則每個(gè)點(diǎn)都各進(jìn)行一次檢測(cè)。實(shí)驗(yàn)中對(duì)信號(hào)進(jìn)行205點(diǎn)的dft，這是因?yàn)闉榱耸姑總€(gè)點(diǎn)都落在整數(shù)點(diǎn)上。題目 （七）一：實(shí)驗(yàn)題目： 該題目的目的是說明一個(gè)pn擴(kuò)頻信號(hào)在抑制

8、正弦干擾中的有效性?，F(xiàn)考慮下圖所示的二進(jìn)制通信系統(tǒng)，對(duì)信號(hào)發(fā)生器的輸入乘上一個(gè)二進(jìn)制（1）pn序列。同一個(gè)二進(jìn)制pn序列用來與解調(diào)器輸入相乘，因此消除了這個(gè)pn序列在期望信號(hào)上的影響。信號(hào)將傳送信號(hào)受到一寬帶加性噪聲序列w(n)和一個(gè)正弦干擾序列i(n)=asinw0n,0w0n)=);w0=pi*0.3;b=4;t=100;n=1:n;p=shuju(n); pq=pseudo_random(b,n);pp=sinxl(a,w0,n);pp1=50*noises(n);d=input(是否加入pn序列?/n 1:是 2:否 0:退出 );if d=1for k=1:19; t=t+n; g=

9、xor(p,pq)*t; x=g+sinxl(a,w0,n); i=1:n; sh=x+pp1; q=sht/3; jtout=xor(q,pq);i=1:n; p0=p(i)=jtout(i); p1=length(find(p0); cc(k)=p1/ni=1:n; me=0; me=me+g(i).2; mp=me/n;j=1:n; se=0; se=se+pp(j).2+pp1(j).2;sp=se/n; sn(k)=mp/sp;endelse if d=2 for k=1:19; t=t+n; g=p*t; x=g+sinxl(a,w0,n); i=1:n; sh=x+pp1; p0

10、=p(i)=sh(i); p1=length(find(p0); cc(k)=p1/n; i=1:n; me=0; me=me+p(i).2; mp=me/n;j=1:n; se=0; se=se+pp(j).2+pp1(j).2; sp=se/n; sn(k)=mp/sp; end; else if d=0 break; end; end;end;%畫信噪比與差錯(cuò)率曲線i=0;figure(i+1);k=1:19;stem(sn(k),cc(k);xlabel(信噪比); ylabel(差錯(cuò)率); title(信噪比與差錯(cuò)率曲線); grid %子函數(shù) %二進(jìn)制數(shù)據(jù)發(fā)生function s=

11、shuju(n)s=rand(1,n);n=1:n;s(n)=s(n)0.5;%產(chǎn)生擴(kuò)頻序列pn function ps_ra=pseudo_random(b,n)a=zeros(1,b); a(b)=1;for n=1:1000 y(n)=a(b);temp=mod(a(b-1)+a(b),2); for j=b:-1:2; a(j)=a(j-1); end a(1)=temp;endps_ra=y(1:n);% w(n)function w=noises(n)w=randn(1,n);%正弦i(n)function i=sinxl(a,w0,n)n=0:n-1;i=a*sin(w0*n);

12、七 程序結(jié)果及圖表 八 實(shí)驗(yàn)總結(jié)相同采樣點(diǎn)數(shù)情況下，隨信噪比增大，差錯(cuò)率減小。在信噪比過大時(shí)，實(shí)際不容易實(shí)現(xiàn)，因此對(duì)信噪比應(yīng)控制在實(shí)際范圍內(nèi)。 不同采樣點(diǎn)數(shù)情況下，隨采樣點(diǎn)數(shù)增加，差錯(cuò)率也增加。信噪比過大時(shí)，差錯(cuò)率無義。在相同信噪比與采樣點(diǎn)數(shù)時(shí)，加上pn序列的差錯(cuò)率小于不加pn序列的差錯(cuò)率，由此可以說明pn擴(kuò)頻系統(tǒng)在正弦干擾信號(hào)存在下優(yōu)于常規(guī)的二進(jìn)制通信系統(tǒng)。題目八：設(shè)原n點(diǎn)序列的dft為驗(yàn)證以下序列的dft與x(k)的關(guān)系：1：, 2：x2(n)=x(2n)，3：x3(n)=x(2n+1)，4：x4(n)=x(n)nr2n(n)，5：x5(n)=x(n)+x(n+n/2)，6：x6(n)=x

13、(n)-x(n+n/2)7：二：實(shí)驗(yàn)?zāi)康?1 根據(jù)數(shù)字信號(hào)處理上學(xué)習(xí)的設(shè)計(jì)濾波器的方法，用matalab語言實(shí)現(xiàn)高通數(shù)字濾波器的設(shè)計(jì)，其中要用雙線性變換法和低通模擬巴特沃斯濾波器法。2 驗(yàn)證一個(gè)pn擴(kuò)頻信號(hào)在抑制正弦干擾中的有效性。3 一個(gè)序列不同條件下的的n/2點(diǎn)和2n點(diǎn)的dft與n點(diǎn)的dft x(k)的關(guān)系。三：實(shí)驗(yàn)原理（1） 設(shè)歸一化低通原形濾波器的系統(tǒng)函數(shù)為hal(p)，p為模擬內(nèi)的拉普拉斯變量。在模擬域內(nèi)從低通向高通變換，即 以1/p 代替p，有hah(p)= hal(1/p)，進(jìn)行反歸一化，即將p=s/c代入上式，得到反歸一化的高通濾波器的傳輸函數(shù)h(s) ，再進(jìn)行雙線性變換，得到

14、數(shù)字高通濾波器的系統(tǒng)函數(shù)為：h(z)=h(s)|s=(1-1/z)/(1+1/z)。（2） 對(duì)信號(hào)發(fā)生器的輸出乘上一個(gè)二進(jìn)制（1）pn序列。同一個(gè)二進(jìn)制pn序列用來與解調(diào)器輸入相乘，因此消除了這個(gè)pn序列在期望信號(hào)上的影響。用和不用pn序列執(zhí)行一個(gè)仿真系統(tǒng)，對(duì)不同的m值測(cè)量差錯(cuò)率，畫出每種情況的差錯(cuò)率曲線，作比較說明這個(gè)pn序列對(duì)于正弦干擾信號(hào)的效果。（3） 通過輸出原序列的x(k)與題目中不同序列的dft得到不同情況下時(shí)域與頻域的變化關(guān)系。如題目一中n為偶數(shù)時(shí)為x(n/2)n為奇數(shù)時(shí)為0，可得知實(shí)域插0，頻域拓展一個(gè)周期。四：實(shí)驗(yàn)步驟簡(jiǎn)述(一)1.確定所需高通數(shù)字濾波器的技術(shù)指標(biāo). 由題目中

15、可知有如下數(shù)據(jù): fs=20000; %采樣頻率 t=1/fs;%采樣周期fs=4000; %數(shù)字阻帶截止頻率 fp=5000;%數(shù)字通帶截止頻率ap=0.5;%通帶最大衰減 as=10;%阻帶最小衰減2.將高通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成低通模擬濾波器的技術(shù)指標(biāo),轉(zhuǎn)換公式為 =2/ ttan(w/2)3.將高通濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技術(shù)指標(biāo)4.設(shè)計(jì)模擬低通濾波器g(p),即確定它的階數(shù)及截止頻率.5.將模擬低通通過頻率變換,轉(zhuǎn)換成模擬高通濾波器,具體是將g(s)的變量換成1/s,得到模擬高通ha(s).6.采用雙線性變換法,將模擬高通ha(s)轉(zhuǎn)換成數(shù)字高通濾波器h(z).（二）1分別

16、設(shè)計(jì)一個(gè)二進(jìn)制序列， pn序列，正弦信號(hào)，噪聲信號(hào)； 2對(duì)二進(jìn)制數(shù)據(jù)發(fā)生器的輸出乘上一個(gè)pn序列，進(jìn)行差分調(diào)制；3調(diào)制后序列加一正弦信號(hào)和噪聲信號(hào)；4設(shè)置限濾除掉噪聲，再將該序列乘上一個(gè)二進(jìn)制（1）pn序列進(jìn)行解調(diào)；5設(shè)計(jì)檢測(cè)器和差錯(cuò)計(jì)數(shù)器，讓解調(diào)后的序列通過檢測(cè)器和差錯(cuò)計(jì)數(shù)器；6對(duì)二進(jìn)制數(shù)據(jù)發(fā)生器的輸出加一正弦信號(hào)和噪聲信號(hào)；7門限濾除掉噪聲；8設(shè)計(jì)檢測(cè)器和差錯(cuò)計(jì)數(shù)器，讓濾掉噪聲的序列通過檢測(cè)器和差錯(cuò)計(jì)數(shù)器；9畫出該有pn序列和沒有pn序列加入系統(tǒng)的信噪比與差錯(cuò)率曲線；（三）1首先應(yīng)該先確定n的數(shù)值，在這里我們?cè)O(shè)為8。求原函數(shù)x(n) 8點(diǎn)dft為x(k) ；輸出為圖象以便于與后面的序列做比

17、較從而得到時(shí)域與頻域之間的變化關(guān)系。2 找出x1(n)與x(n)的關(guān)系，確定其長(zhǎng)度為m并寫出對(duì)應(yīng)各點(diǎn)的值，3 讓其進(jìn)行m點(diǎn) 的dft變化得到x1(k)；4 畫出函數(shù)x(n)，x1(n)，x(k)和x1(k)的圖形，分析圖形之間的關(guān)系；5 得出結(jié)論；6 其他程序同上。五：程序框圖開始頻率歸一化將數(shù)字指標(biāo)轉(zhuǎn)換成模擬指標(biāo)設(shè)計(jì)滿足條件的模擬低通濾波將低通模擬轉(zhuǎn)換成高通模擬將高通模擬轉(zhuǎn)換成高通數(shù)字繪制該濾波器的幅度頻譜圖結(jié)束開始x(n)的n點(diǎn)dft x(k)x1(n)各點(diǎn)的值確定長(zhǎng)度m和nx1(n)的n點(diǎn)dft x1(k)畫出函數(shù)x(n)，x1(n)，x(k)和x1(k)的圖形結(jié)束開 始調(diào)制器每比特樣本

18、數(shù)選擇是否經(jīng)過pn序列二進(jìn)制序列s和pn序列，讓兩個(gè)序列相乘進(jìn)行得到g序列正弦序列q噪聲序列re,與g序列相加，得到x1序列一門限濾波器， x1序列經(jīng)過該門限濾除掉噪聲得到序列r2讓r2與pn序列相乘進(jìn)行解調(diào)得到序列r3結(jié) 束檢測(cè)器和差錯(cuò)計(jì)數(shù)器，讓解調(diào)后的序列通過檢測(cè)器和差錯(cuò)記數(shù)器二進(jìn)制序列與原始信號(hào)的幅度因子相乘得到序列h正弦序列q和噪聲信號(hào)re,，并讓h序列與這兩個(gè)序列相加得到序列x2門限濾波器， x2序列經(jīng)過該門限濾除掉噪聲得到序列r6六：源程序 4. fs=20000;t=1/fs;fs=4000;fp=5000;ap=0.5;as=10;ws=(fs/fs)*2*pi;wp=(fp/

19、fs)*2*pi; %將已知頻率歸一化wp=2*tan(wp/2)/t;ws=2*tan(ws/2)/t; %將頻率非線性變換n,wn=buttord(wp,ws,ap,as,s); % h模擬低通濾波器的階和截止頻率 fprintf(濾波器階數(shù)n=%.0fn,n); % butterworth模擬低通濾波器的階nz,p,k=buttap(n); %butterworth模擬低通濾波器原形num,den=zp2tf(z,p,k); %變零極點(diǎn)增益形式為傳遞函數(shù)形式numa,dena=lp2hp(num,den,wn); %模擬低通轉(zhuǎn)換到模擬高通numd,dend=bilinear(numa,d

20、ena,fs); %雙線性變換法disp(分子系數(shù)numd);fprintf(%.4e ,numd);fprintf(n);disp(分母系數(shù)dend);fprintf(%.4e ,dend);fprintf(n);figure(name,實(shí)驗(yàn)四,color,w); h,w=freqz(numd,dend,100); %計(jì)算數(shù)字高通濾波器的頻率響應(yīng)（n=100點(diǎn)均勻分布在上半單位圓0-pi，n點(diǎn)頻率記錄在w中，對(duì)應(yīng)的頻率響應(yīng)記錄在h中）plot(w*2*fs/(pi),abs(h);axis(0 6000 0 1.5);xlabel( 頻率 f/hz );ylabel(歸一化 幅度 );tit

21、le(高通濾波器);grid on;7. function =ti7()w0=0.2*pi;n=4;while(1)n=input(請(qǐng)選擇樣本數(shù)m：n1-m=50n2-m=100n3-m=500n4-m=1000n0-quitn); if n=1 m=50; %原始信號(hào)的長(zhǎng)度 a=60; %正弦信號(hào)的幅度 z=50; %原始信號(hào)的幅度因子 b=20; %噪聲幅度 p(n,a,m,z,b,w0);elseif n=2 m=100; a=110; z=100; b=40; p(n,a,m,z,b,w0); elseif n=3 m=500; a=550; z=500; b=80; p(n,a,m,

22、z,b,w0); elseif n=4 m=1000; a=1100; z=1000; b=130; p(n,a,m,z,b,w0); elseif n=0 break;end end function =p(n,a,m,z,b,w0)s=randout(m);r=pn(n,m);for k=1:19 %循環(huán)增加z（信號(hào)的幅度） z=z+b; %使調(diào)制后的序列幅度z大于噪聲幅度b g=xor(s,r)*z; %調(diào)制（邏輯異或）差分調(diào)制 q=sinout(a,w0,m); x=g+q; %調(diào)制后序列加正弦信號(hào) re=b*wout(m); %生成噪聲信號(hào) x1=re+x; %加入噪聲 r2=x1z

23、/3; %設(shè)置一個(gè)門限，濾除噪聲 r3=xor(r2,r); %解調(diào) i=1:m; r4=s(i)=r3(i); %檢測(cè) （當(dāng)不等時(shí)取1 且計(jì)入r4 當(dāng)相等時(shí)取0 不計(jì)入r4） r5=length(find(r4); %計(jì)算錯(cuò)誤的個(gè)數(shù)（計(jì)算r4中元素取1的個(gè)數(shù)） cc(k)=r5/m; %錯(cuò)誤的概率 i=1:m; me=0; me=me+g(i).2; %計(jì)算信號(hào)的總能量 mp=me/m; %計(jì)算信號(hào)的平均功率 j=1:m; se=0; se=se+q(j).2+re(j).2;%計(jì)算噪聲的總能量 sp=se/m; %計(jì)算噪聲的平均功率 sn(k)=mp/sp; %計(jì)算信噪比 %不加pn序列

24、h=s*z; %原始信號(hào)乘幅度因子 x2=h+re+q; %信號(hào)序列加正弦信號(hào)和噪聲信號(hào)（不與pn序列調(diào)制） r6=x2z/3; %設(shè)置一個(gè)門限，濾除噪聲 i=1:m; r7=s(i)=r6(i); %檢測(cè) （當(dāng)不等時(shí)取1 且計(jì)入r7 當(dāng)相等時(shí)取0 不計(jì)入r7） r8=length(find(r7); %計(jì)算錯(cuò)誤的個(gè)數(shù)（計(jì)算r7中元素取1的個(gè)數(shù)） cd(k)=r8/m; %錯(cuò)誤的概率 i=1:m; ma=0; ma=ma+h(i).2; %計(jì)算信號(hào)的總能量 md=ma/m; %計(jì)算信號(hào)的平均功率 j=1:m; sa=0; sa=sa+q(j).2+re(j).2;%計(jì)算噪聲的總能量 st=sa

25、/m; %計(jì)算噪聲的平均功率 sd(k)=md/st; %計(jì)算信噪比end figure; stem(sn,cc,.); %輸出有pn序列的圖hold on;stem(sd,cd,-r); %輸出沒有pn序列的圖xlabel(信噪比); ylabel(差錯(cuò)率);title(信噪比與差錯(cuò)率曲線); grid %打開網(wǎng)絡(luò)線function s=randout(m)n=0:m;s=rand(1,m); %1*m的0-1的均勻分布的隨機(jī)數(shù)矩陣j=1:m;s(j)=s(j)0.5; %使每個(gè)元素的值大于0.5 function p_s=pn(n,m) a=zeros(1,n); %1*n的零矩陣 a(n

26、)=1; %令第n列數(shù)為1 for i=1:1000 y(i)=a(n); temp=mod(a(n-1)+a(n),2);%除法求余（與除數(shù)同號(hào)） for j=n:-1:2 a(j)=a(j-1); %改變a(n)序列 end a(1)=temp; %改變a(n)序列 end p_s=y(1:m); %取a(n)序列的最后一位元素生成pn序列function in=sinout(a,w0,m)j=1:m;in=a*sin(w0*j);function w=wout(m)w=rand(1,m); %1*m的0-1的均勻分布的隨機(jī)數(shù)矩陣 8.dft8:n=8; a=1.2; n=0:7; xn=

27、a.n; xk=dft8(xn,n); subplot(3,1,1) stem(n,xn,.k);axis(0,8,0,4.5) subplot(3,1,2) stem(n,abs(xk),.k);axis(0,8,0,10) subplot(3,1,3) stem(n,angle(xk),.k);axis(0,8,-4.5,4.5)x1n: n=8; a=1.2; n=0:15; m=2*n; xn=a.n; yn=xn(:,1) 0 xn(:,2) 0 xn(:,3) 0 xn(:,4) 0 xn(:,5) 0 xn(:,6) 0 xn(:,7) 0 xn(:,8) 0; yk=dft8(

28、yn,m); subplot(3,1,1) stem(n,yn,.k);axis(0,16,0,5) subplot(3,1,2) stem(n,abs(yk),.k);axis(0,16,0,5) subplot(3,1,3) stem(n,angle(yk),.k);axis(0,16,-4.5,4.5)x2n: n=8; a=1.2; n=0:7; m=n/2; xn=a.n; m=0:3; x2n=xn(:,1) xn(:,3) xn(:,5) xn(:,7); x2k=dft8(x2n,m); subplot(3,1,1) stem(m,x2n,.k);axis(0,3,0,1.5)

29、 subplot(3,1,2) stem(m,abs(x2k),.k);axis(0,3,0,5) subplot(3,1,3) stem(m,angle(x2k),.k);axis(0,3,-1.5,1.5)x3n: n=8; a=1.2; n=0:7; m=n/2; xn=a.n; m=0:3; x3n=xn(:,2) xn(:,4) xn(:,6) xn(:,8); x3k=dft8(x3n,m); subplot(3,1,1) stem(m,x3n,.k);axis(0,3,0,5) subplot(3,1,2) stem(m,abs(x3k),.k);axis(0,3,0,5) su

30、bplot(3,1,3) stem(m,angle(x3k),.k);axis(0,3,-4.5,4.5)x4n: n=8; a=1.2; n=0:15; m=2*n; xn=a.n; x4n=xn(:,1) xn(:,2) xn(:,3) xn(:,4) xn(:,5) xn(:,6) xn(:,7) xn(:,8) xn(:,1) xn(:,2) xn(:,3) xn(:,4) xn(:,5) xn(:,6) xn(:,7) xn(:,8); x4k=dft8(x4n,m); subplot(3,1,1) stem(n,x4n,.k);axis(0,16,0,5) subplot(3,1,

31、2) stem(n,abs(x4k),.k);axis(0,16,0,5) subplot(3,1,3) stem(n,angle(x4k),.k);axis(0,16,-8,8)x5n: n=8; a=1.2; n=0:7; m=n/2; xn=a.n; m=0:3; x5n=xn(:,1)+xn(:,5) xn(:,2)+xn(:,6) xn(:,3)+xn(:,7) xn(:,4)+xn(:,8); x5k=dft8(x5n,m); subplot(3,1,1) stem(m,x5n,.k);axis(0,3,0,5) subplot(3,1,2) stem(m,abs(x5k),.k);axis(0,3,0,5) subplot(3,1,3) stem(m,angle(x5k),.k);axis(0,3,-4.5,4.5)x6n: n=8; a=1.2; n=0:7; m=n/2; xn=a.n; m=0:3; wn=exp(-j*2*pi/n); wnnk=wn.m; yn=xn(:,1)-xn(:,5) xn(:,2)-xn(:,6) xn(:,3)-xn(:,7) xn(:,4)-xn(:,