查看: 441|回复: 6|关注: 0

[已解决] 下面代码怎么改写输入,求大神

[复制链接]

新手

12 麦片

财富积分


050


1

主题

5

帖子

0

最佳答案
%% 清空环境变量
clc
clear
%% 网络参数配置
load traffic_flux input output input_test output_test
M=size(input,2); %输入节点个数
N=size(output,2); %输出节点个数
n=6; %隐形节点个数
lr1=0.01; %学习概率
lr2=0.001; %学习概率
maxgen=1000; %迭代次数
%权值初始化
Wjk=randn(n,M);Wjk_1=Wjk;Wjk_2=Wjk_1;
Wij=randn(N,n);Wij_1=Wij;Wij_2=Wij_1;
a=randn(1,n);a_1=a;a_2=a_1;
b=randn(1,n);b_1=b;b_2=b_1;
%节点初始化
y=zeros(1,N);
net=zeros(1,n);
net_ab=zeros(1,n);
%权值学习增量初始化
d_Wjk=zeros(n,M);
d_Wij=zeros(N,n);
d_a=zeros(1,n);
d_b=zeros(1,n);
%% 输入输出数据归一化
[inputn,inputps]=mapminmax(input');
[outputn,outputps]=mapminmax(output');
inputn=inputn';
outputn=outputn';
%% 网络训练
for i=1:maxgen
   
    %误差累计
    error(i)=0;
   
    % 循环训练
    for kk=1:size(input,1)
        x=inputn(kk,:);
        yqw=outputn(kk,:);
   
        for j=1:n
            for k=1:M
                net(j)=net(j)+Wjk(j,k)*x(k);
                net_ab(j)=(net(j)-b(j))/a(j);
            end
            temp=mymorlet(net_ab(j));
            for k=1:N
                y=y+Wij(k,j)*temp;   %小波函数
            end
        end
        
        %计算误差和
        error(i)=error(i)+sum(abs(yqw-y));
        
        %权值调整
        for j=1:n
            %计算d_Wij
            temp=mymorlet(net_ab(j));
            for k=1:N
                d_Wij(k,j)=d_Wij(k,j)-(yqw(k)-y(k))*temp;
            end
            %计算d_Wjk
            temp=d_mymorlet(net_ab(j));
            for k=1:M
                for l=1:N
                    d_Wjk(j,k)=d_Wjk(j,k)+(yqw(l)-y(l))*Wij(l,j) ;
                end
                d_Wjk(j,k)=-d_Wjk(j,k)*temp*x(k)/a(j);
            end
            %计算d_b
            for k=1:N
                d_b(j)=d_b(j)+(yqw(k)-y(k))*Wij(k,j);
            end
            d_b(j)=d_b(j)*temp/a(j);
            %计算d_a
            for k=1:N
                d_a(j)=d_a(j)+(yqw(k)-y(k))*Wij(k,j);
            end
            d_a(j)=d_a(j)*temp*((net(j)-b(j))/b(j))/a(j);
        end
        
        %权值参数更新      
        Wij=Wij-lr1*d_Wij;
        Wjk=Wjk-lr1*d_Wjk;
        b=b-lr2*d_b;
        a=a-lr2*d_a;
   
        d_Wjk=zeros(n,M);
        d_Wij=zeros(N,n);
        d_a=zeros(1,n);
        d_b=zeros(1,n);
        y=zeros(1,N);
        net=zeros(1,n);
        net_ab=zeros(1,n);
        
        Wjk_1=Wjk;Wjk_2=Wjk_1;
        Wij_1=Wij;Wij_2=Wij_1;
        a_1=a;a_2=a_1;
        b_1=b;b_2=b_1;
    end
end
%% 网络预测
%预测输入归一化
x=mapminmax('apply',input_test',inputps);
x=x';
%网络预测
for i=1:92
    x_test=x(i,:);
    for j=1:1:n
        for k=1:1:M
            net(j)=net(j)+Wjk(j,k)*x_test(k);
            net_ab(j)=(net(j)-b(j))/a(j);
        end
        temp=mymorlet(net_ab(j));
        for k=1:N
            y(k)=y(k)+Wij(k,j)*temp ;
        end
    end
    yuce(i)=y(k);
    y=zeros(1,N);
    net=zeros(1,n);
    net_ab=zeros(1,n);
end
%预测输出反归一化
ynn=mapminmax('reverse',yuce,outputps);
%% 结果分析
figure(1)
plot(ynn,'r*:')
hold on
plot(output_test,'bo--')
title('预测交通流量','fontsize',12)
legend('预测交通流量','实际交通流量')
xlabel('时间点')
ylabel('交通流量')
回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1635

帖子

343

最佳答案
  • 关注者: 81
发表于 2020-2-13 15:07:53 | 显示全部楼层
load traffic_flux input output input_test output_test这行代码即为输入训练样本和测试样本
回复此楼 已获打赏: 0 积分

举报

新手

12 麦片

财富积分


050


1

主题

5

帖子

0

最佳答案
 楼主| 发表于 2020-2-13 15:12:18 | 显示全部楼层
20141303 发表于 2020-2-13 15:07
load traffic_flux input output input_test output_test这行代码即为输入训练样本和测试样本 ...

谢谢回复,那怎么改变输入呢,就是我想自己重新输入数据
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1635

帖子

343

最佳答案
  • 关注者: 81
发表于 2020-2-13 22:22:53 | 显示全部楼层
取决于你想怎么输入,导入excel ,mat还是直接在脚本文件中手打样本
回复此楼 已获打赏: 0 积分

举报

新手

12 麦片

财富积分


050


1

主题

5

帖子

0

最佳答案
 楼主| 发表于 2020-2-14 13:44:54 | 显示全部楼层
20141303 发表于 2020-2-13 22:22
取决于你想怎么输入,导入excel ,mat还是直接在脚本文件中手打样本

好的,那导入excel,mat的话具体怎么操作呢,是要编写别的代码吗,比如我想输入三组数据,是一个20*3的矩阵。
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1635

帖子

343

最佳答案
  • 关注者: 81
发表于 2020-2-14 19:31:42 | 显示全部楼层 |此回复为最佳答案
导入excel使用xlsread函数,导入mat使用load函数,具体查询help
回复此楼 已获打赏: 0 积分

举报

新手

12 麦片

财富积分


050


1

主题

5

帖子

0

最佳答案
 楼主| 发表于 2020-2-19 21:34:35 | 显示全部楼层
20141303 发表于 2020-2-14 19:31
导入excel使用xlsread函数,导入mat使用load函数,具体查询help

好的,感谢您的帮助
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表