ziguanzuichufangzheng.m

本程序是用AR模型系数与LM神经网络完成微弱信号识别的程序

clear all
fp1=fopen('e:\zyh\z01.dat','r');z01=fscanf(fp1,'%f');fclose(fp1);z01=z01(1:90);z01=z01-sum(z01)/90;
fp1=fopen('e:\zyh\z02.dat','r');z02=fscanf(fp1,'%f');fclose(fp1);z02=z02(1:90);z02=z02-sum(z02)/90;
fp1=fopen('e:\zyh\z03.dat','r');z03=fscanf(fp1,'%f');fclose(fp1);z03=z03(1:90);z03=z03-sum(z03)/90;
fp1=fopen('e:\zyh\z04l.dat','r');z04=fscanf(fp1,'%f');fclose(fp1);z04=z04(1:90);z04=z04-sum(z04)/90;
fp1=fopen('e:\zyh\z05l.dat','r');z05=fscanf(fp1,'%f');fclose(fp1);z05=z05(1:90);z05=z05-sum(z05)/90;
fp1=fopen('e:\zyh\z06l.dat','r');z06=fscanf(fp1,'%f');fclose(fp1);z06=z06(1:90);z06=z06-sum(z06)/90;
fp1=fopen('e:\zyh\z07.dat','r');z07=fscanf(fp1,'%f');fclose(fp1);z07=z07(1:90);z07=z07-sum(z07)/90;
fp1=fopen('e:\zyh\z08.dat','r');z08=fscanf(fp1,'%f');fclose(fp1);z08=z08(1:90);z08=z08-sum(z08)/90;
fp1=fopen('e:\zyh\z09.dat','r');z09=fscanf(fp1,'%f');fclose(fp1);z09=z09(1:90);z09=z09-sum(z09)/90;
fp1=fopen('e:\zyh\z10.dat','r');z10=fscanf(fp1,'%f');fclose(fp1);z10=z10(1:90);z10=z10-sum(z10)/90;
fp1=fopen('e:\zyh\z11.dat','r');z11=fscanf(fp1,'%f');fclose(fp1);z11=z11(1:90);z11=z11-sum(z11)/90;
fp1=fopen('e:\zyh\z12l.dat','r');z12=fscanf(fp1,'%f');fclose(fp1);z12=z12(1:90);z12=z12-sum(z12)/90;
fp1=fopen('e:\zyh\z13l.dat','r');z13=fscanf(fp1,'%f');fclose(fp1);z13=z13(1:90);z13=z13-sum(z13)/90;
fp1=fopen('e:\zyh\z14l.dat','r');z14=fscanf(fp1,'%f');fclose(fp1);z14=z14(1:90);z14=z14-sum(z14)/90;
fp1=fopen('e:\zyh\z15.dat','r');z15=fscanf(fp1,'%f');fclose(fp1);z15=z15(1:90);z15=z15-sum(z15)/90;
fp1=fopen('e:\zyh\b01.dat','r');b01=fscanf(fp1,'%f');fclose(fp1);b01=b01(1:90);b01=b01-sum(b01)/90;
fp1=fopen('e:\zyh\b02.dat','r');b02=fscanf(fp1,'%f');fclose(fp1);b02=b02(1:90);b02=b02-sum(b02)/90;
fp1=fopen('e:\zyh\b03.dat','r');b03=fscanf(fp1,'%f');fclose(fp1);b03=b03(1:90);b03=b03-sum(b03)/90;
fp1=fopen('e:\zyh\b04.dat','r');b04=fscanf(fp1,'%f');fclose(fp1);b04=b04(1:90);b04=b04-sum(b04)/90;
fp1=fopen('e:\zyh\b05.dat','r');b05=fscanf(fp1,'%f');fclose(fp1);b05=b05(1:90);b05=b05-sum(b05)/90;
fp1=fopen('e:\zyh\b06.dat','r');b06=fscanf(fp1,'%f');fclose(fp1);b06=b06(1:90);b06=b06-sum(b06)/90;
fp1=fopen('e:\zyh\b07.dat','r');b07=fscanf(fp1,'%f');fclose(fp1);b07=b07(1:90);b07=b07-sum(b07)/90;
fp1=fopen('e:\zyh\b08.dat','r');b08=fscanf(fp1,'%f');fclose(fp1);b08=b08(1:90);b08=b08-sum(b08)/90;
fp1=fopen('e:\zyh\b09.dat','r');b09=fscanf(fp1,'%f');fclose(fp1);b09=b09(1:90);b09=b09-sum(b09)/90;
fp1=fopen('e:\zyh\b10.dat','r');b10=fscanf(fp1,'%f');fclose(fp1);b10=b10(1:90);b10=b10-sum(b10)/90;
fp1=fopen('e:\zyh\b11.dat','r');b11=fscanf(fp1,'%f');fclose(fp1);b11=b11(1:90);b11=b11-sum(b11)/90;
fp1=fopen('e:\zyh\b12.dat','r');b12=fscanf(fp1,'%f');fclose(fp1);b12=b12(1:90);b12=b12-sum(b12)/90;
fp1=fopen('e:\zyh\b13.dat','r');b13=fscanf(fp1,'%f');fclose(fp1);b13=b13(1:90);b13=b13-sum(b13)/90;
fp1=fopen('e:\zyh\b14.dat','r');b14=fscanf(fp1,'%f');fclose(fp1);b14=b14(1:90);b14=b14-sum(b14)/90;
fp1=fopen('e:\zyh\b15.dat','r');b15=fscanf(fp1,'%f');fclose(fp1);b15=b15(1:90);b15=b15-sum(b15)/90;
N=90;
%正常人参数
n=3;
    Y=z01(n+1:N);
    m=N-n;
   for i=1:m
        for j=1:n
            X(i,j)=z01(n+i-j);
    end
end
p1=inv(X'*X)*(X'*Y);
    Y=z02(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z02(n+i-j);
    end
end
p2=inv(X'*X)*(X'*Y);
    Y=z03(n+1:N);
    m=N-n;
    
    for i=1:m
        for j=1:n
            X(i,j)=z03(n+i-j);
    end
end
p3=inv(X'*X)*(X'*Y);
    Y=z04(n+1:N);
    m=N-n;
    
    for i=1:m
        for j=1:n
            X(i,j)=z04(n+i-j);
    end
end
p4=inv(X'*X)*(X'*Y);
    Y=z05(n+1:N);
    m=N-n;   
    for i=1:m
        for j=1:n
            X(i,j)=z05(n+i-j);
    end
end
p5=inv(X'*X)*(X'*Y);
    Y=z06(n+1:N);
    m=N-n;
  
    for i=1:m
        for j=1:n
            X(i,j)=z06(n+i-j);
    end
end
p6=inv(X'*X)*(X'*Y);
    Y=z07(n+1:N);
    m=N-n;
    
    for i=1:m
        for j=1:n
            X(i,j)=z07(n+i-j);
    end
end
p7=inv(X'*X)*(X'*Y);
    Y=z08(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z08(n+i-j);
    end
end
p8=inv(X'*X)*(X'*Y);
    Y=z09(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z09(n+i-j);
    end
end
p9=inv(X'*X)*(X'*Y);
    Y=z10(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z10(n+i-j);
    end
end
p10=inv(X'*X)*(X'*Y);
    Y=z11(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z11(n+i-j);
    end
end
p11=inv(X'*X)*(X'*Y);
    Y=z12(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z12(n+i-j);
    end
end
p12=inv(X'*X)*(X'*Y);
    Y=z13(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z13(n+i-j);
    end
end
p13=inv(X'*X)*(X'*Y);
    Y=z14(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z14(n+i-j);
    end
end
p14=inv(X'*X)*(X'*Y);
    Y=z15(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=z15(n+i-j);
    end
end
p15=inv(X'*X)*(X'*Y);
%非正常人参数
    Y=b01(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b01(n+i-j);
    end
end
p16=inv(X'*X)*(X'*Y);
    Y=b02(n+1:N);
    m=N-n;
   for i=1:m
        for j=1:n
            X(i,j)=b02(n+i-j);
    end
end
p17=inv(X'*X)*(X'*Y);
    Y=b03(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b03(n+i-j);
    end
end
p18=inv(X'*X)*(X'*Y);
    Y=b04(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b04(n+i-j);
    end
end
p19=inv(X'*X)*(X'*Y);
    Y=b05(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b05(n+i-j);
    end
end
p20=inv(X'*X)*(X'*Y);
    Y=b06(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b06(n+i-j);
    end
end
p21=inv(X'*X)*(X'*Y);
    Y=b07(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b07(n+i-j);
    end
end
p22=inv(X'*X)*(X'*Y);
    Y=b08(n+1:N);
    m=N-n;
     for i=1:m
        for j=1:n
            X(i,j)=b08(n+i-j);
    end
end
p23=inv(X'*X)*(X'*Y);
    Y=b09(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b09(n+i-j);
    end
end
p24=inv(X'*X)*(X'*Y);
    Y=b10(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b10(n+i-j);
    end
end
p25=inv(X'*X)*(X'*Y);
    Y=b11(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b11(n+i-j);
    end
end
p26=inv(X'*X)*(X'*Y);
    Y=b12(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b12(n+i-j);
    end
end
p27=inv(X'*X)*(X'*Y);
    Y=b13(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b13(n+i-j);
    end
end
p28=inv(X'*X)*(X'*Y);
    Y=b14(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b14(n+i-j);
    end
end
p29=inv(X'*X)*(X'*Y);
    Y=b15(n+1:N);
    m=N-n;
    for i=1:m
        for j=1:n
            X(i,j)=b15(n+i-j);
    end
end
p30=inv(X'*X)*(X'*Y);
%网络仿真
p=[p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 p16 p17 p18 p19 p20 p21 p22 p23 p24 p25 ];
ps=[p11 p12 p13 p14 p15 p26 p27 p28 p29 p30];
t=[1;1 ;1 ;1 ;1 ;1 ;1 ;1 ;1 ;1 ;-1;-1 ;-1 ;-1 ;-1 ;-1 ;-1 ;-1 ;-1 ;-1 ]';
net=newff([0 4;-5 0;-2 4],[7 1],{'tansig','tansig' },'trainlm');
net.trainparam.goal=0.001;
net.trainparam.epochs=5000;
[net,tr]=train(net,p,t);
Y=sim(net,p);
YS=sim(net,ps)