原来是小飞 发表于 2013-10-28 09:18

MATLAB 代码的详细解释问题

如下的两个matlab程序代码，本人能看懂一小部分，搞不清楚它的意思，希望有人帮忙解释一下每一句代码的详细解释，不胜感激。clear all
data=load('b0.txt');   
%A=data(:,);
A=data(100:300,);
m1=200;
m2=m1+1;
for i=1:m2
for j=1:24
    X((i-1)*24+j,1)=A(i,j);
   end
end
roll=X;%取X数据
len=m1*24;
oldroll=roll;
AverRoll=roll-mean(roll);
roll=AverRoll;
mtbp=6;
taow=mtbp;
deltaT=0.5;
tao=6;
taowdivtao=round(taow/tao);
d=2;%taowdivtao+1;
prestep=24;
N = length(roll)-prestep;
MaxStep = 24;
Step = 1 : 1 : MaxStep;
= FunctionChaosPredict(roll(1:N),N,mtbp,deltaT,tao,d,MaxStep);
subplot(2,1,1);plot(roll(1:len+MaxStep)+mean(oldroll));
hold on;
subplot(2,1,1);
plot(Step + N,PredictedData+mean(oldroll),'r:');
hold on;
subplot(2,1,2);
plot(Step + N,PredictedData+mean(oldroll),'r:');
hold on;
subplot(2,1,2);
plot(Step + N,oldroll(Step + N),'b-');%AOLMM多步预测函数
function = FunctionChaosPredict(Data,N,mtbp,deltaT,tao,d,MaxStep)
%Data是一维信号时间序列，N是信号数据长度，mtbp,deltaT,tao,d分别是重构相空间的平均时间序列、采样周期、时延及嵌入维

roll=Data;%取横摇数据
M = N - (d - 1)*tao;
for i = 1 : M
    for j = 1 : d
      MatrixX(i,j) = roll(i + (j - 1)*tao);
    end
end

%计算相空间中第M点与各点的距离
for j = 1 : (M - 1)
         Dis(j) = norm(MatrixX(M,:) - MatrixX(j,:),2);
end
%排序计算相空间中第M点的(m+1)个参考邻近点
for i = 1 : (d + 1)
    NearDis(i) = Dis(i);
    NearPos(i) = i;
end
for i = (d + 2) : (M - mtbp)
    for j = 1 : (d + 1)
      if (abs(i-j)>mtbp) %& abs(i-j)<10*mtbp
      if(Dis(i) < NearDis(j))
            NearDis(j) = Dis(i);
            NearPos(j) = i;
            break;
      end
      end
    end
end
SortedDis = sort(NearDis);
MinDis = SortedDis(1);
%计算第M点的(m+1)个参考邻近点的权P
SumP = 0;
for i = 1 : (d + 1)
    P(i) = exp(-NearDis(i)/MinDis);
    SumP = SumP + P(i);
end
P = P/SumP;
%用最小二乘法计算a[],b[]
for step=1:1:MaxStep
    aCoe1 = 0;
    aCoe2 = d;
    bCoe1 = 0;
    bCoe2 = 0;
    e = 0;
    f = 0;
    for i = 1 : (d + 1)
      aCoe1 = aCoe1 + P(i)*sum(MatrixX(NearPos(i),:));
      bCoe1 = bCoe1 + P(i)*(MatrixX(NearPos(i),:)*MatrixX(NearPos(i),:)');
      e = e + P(i)*(MatrixX(NearPos(i) + step,:)*MatrixX(NearPos(i),:)');
      f = f + P(i)*sum(MatrixX(NearPos(i) + step,:));
    end
    bCoe2 = aCoe1;
    CoeMatrix = ;
    ResultMatrix = ;
    abResult = pinv(CoeMatrix)*ResultMatrix;
    a = abResult(1);
    b = abResult(2);
for j = 1 : d
%   MatrixX(M + step,j) = a + b*MatrixX(M,j); %以历史上相近点的演化规律作为中心点的演化规律以中心点为基准进行预报
   
    MatrixX(M + step,j) = 0;
    for i = 1 : (d + 1)
      MatrixX(M + step,j) = MatrixX(M + step,j) + P(i)*(a + b*MatrixX(NearPos(i),j)); %以历史上相近点的演化加权和直接作为中心点的演化点进行预报
    end
end
%误差修正
    if M-tao+step+(d-1)*tao < N+1
      for j=1:d-1
            err(j)=MatrixX(M + step,j)-roll(M+step+(j-1)*tao);
      end
      ppp=1:d-1;ttt=err;neterr=newrbe(ppp,ttt);xxx=2:d;errp=sim(neterr,xxx);
      PredictedData(step) = MatrixX(M + step,d) - errp(d-1);
      roll(N+step)=PredictedData(step);
    else    PredictedData(step) = MatrixX(M + step,d);
    end%      roll(N+k)=PredictedData(k);
FChaosPredict(step) = PredictedData(step);
% FChaosPredict(step) = MatrixX(M + step,d);
end

ChaChing 发表于 2013-11-7 14:48

每一句代码的详细解释
这要求有点难...

muzhicai 发表于 2013-11-8 12:05

这个全部解释有点困难吧，我学习别人程序就是一边运行，一边学习

原来是小飞 发表于 2013-11-8 17:10

ChaChing 发表于 2013-11-7 14:48 static/image/common/back.gif
这要求有点难...

嗯，一边运行，一边学习吧

原来是小飞 发表于 2013-11-8 17:10

muzhicai 发表于 2013-11-8 12:05 static/image/common/back.gif
这个全部解释有点困难吧，我学习别人程序就是一边运行，一边学习

您的方法不错，支持，学习中。

牛小贱 发表于 2013-12-7 10:34

原来是小飞 发表于 2013-11-8 17:10 static/image/common/back.gif
您的方法不错，支持，学习中。

这确实是个好办法，其实学习Matlab编程，关键是调试里面的Bug，只是最重要的。遇到bug，自己耐心调试，不会的可以百度之类的……我之前学习编程，就是慢慢的调试运行！！

drago0119 发表于 2013-12-8 21:04

学习学习好方法

尤今 发表于 2013-12-8 22:06

嗯，在学习之中。。
页: [1]
查看完整版本: MATLAB 代码的详细解释问题