【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 7973  |  回复: 3
【悬赏金币】回答本帖问题,作者caldys将赠送您 15 个金币

caldys

新虫 (初入文坛)

[求助] 神经网络 NARX已有1人参与

利用NARX神经网络进行模拟仿真,所编程序如下:
clear
clc

%% 加载数据
load phdata
inputSeries = phInputs;
targetSeries = phTargets;

%% 建立非线性自回归模型
inputDelays = 1:2;
feedbackDelays = 1:2;
hiddenLayerSize = 10;
net = narxnet(inputDelays,feedbackDelays,hiddenLayerSize);

%% 网络数据预处理函数定义
net.inputs{1}.processFcns = {'removeconstantrows','mapminmax'};
net.inputs{2}.processFcns = {'removeconstantrows','mapminmax'};

%% 时间序列数据准备工作
[inputs,inputStates,layerStates,targets] = preparets(net,inputSeries,{},targetSeries);

%% 训练数据、验证数据、测试数据划分
net.divideFcn = 'dividerand';  
net.divideMode = 'value';  
net.divideParam.trainRatio = 70/100;
net.divideParam.valRatio = 15/100;
net.divideParam.testRatio = 15/100;

%% 网络训练函数设定
net.trainFcn = 'trainlm';  % Levenberg-Marquardt

%% 误差函数设定
net.performFcn = 'mse';  % Mean squared error

%% 绘图函数设定
net.plotFcns = {'plotperform','plottrainstate','plotresponse', ...
  'ploterrcorr', 'plotinerrcorr'};

%% 网络训练
[net,tr] = train(net,inputs,targets,inputStates,layerStates);

%% 网络测试
outputs = net(inputs,inputStates,layerStates);
errors = gsubtract(targets,outputs);
performance = perform(net,targets,outputs)

%% 计算训练集、验证集、测试集误差
trainTargets = gmultiply(targets,tr.trainMask);
valTargets = gmultiply(targets,tr.valMask);
testTargets = gmultiply(targets,tr.testMask);
trainPerformance = perform(net,trainTargets,outputs)
valPerformance = perform(net,valTargets,outputs)
testPerformance = perform(net,testTargets,outputs)

%% 网络训练效果可视化
figure, plotperform(tr)
figure, plottrainstate(tr)
figure, plotregression(targets,outputs)
figure, plotresponse(targets,outputs)
figure, ploterrcorr(errors)
figure, plotinerrcorr(inputs,errors)


%% close loop模式的实现
% 更改NARX神经网络模式
narx_net_closed = closeloop(net);
view(net)
view(narx_net_closed)

% 计算1500-2000个点的拟合效果
phInputs_c=phInputs(1500:2000);
PhTargets_c=phTargets(1500:2000);

[p1,Pi1,Ai1,t1] = preparets(narx_net_closed,phInputs_c,{},PhTargets_c);
% 网络仿真
yp1 = narx_net_closed(p1,Pi1,Ai1);
plot([cell2mat(yp1)' cell2mat(t1)'])


运行后出现如下错误:如何解决
Undefined function or method 'narxnet' for input arguments of type 'double'.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苗又又

新虫 (初入文坛)

老铁啊,你的数据类型都不转换的嘛。。那个double人家处理不了,你要加个num2cell改成元胞数组存储。。
2楼2018-03-24 19:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苗又又

新虫 (初入文坛)

哦不好意思我说错了,应该是tonndata函数把输入输出的数据改成元胞数组储存,否则没法干活儿
3楼2018-03-24 20:07:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jonchelle

新虫 (初入文坛)

【答案】应助回帖

[Xs,Xi,Ai,Ts] = preparets(net,InputCell,{},TargetCell)

函数内Input和Target必须是cell格式,可以使用tonndata()转换为标准NN参数形式
4楼2019-06-20 20:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 caldys 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[育儿交流] 小孩7岁,上一年级,不自信,学习较差 10+7 ljaxida 2024-04-16 36/1800 2024-04-20 07:19 by ljaxida
[论文投稿] 催稿后秒拒 +5 lizhengke06 2024-04-19 6/300 2024-04-19 23:44 by peterrjp
[考博] 申请24博士 材料/化工/环境 +4 满目_星辰 2024-04-17 4/200 2024-04-19 20:10 by 前行的道路
[论文投稿] 一审一个审稿人,小修,会怎么样呀? +9 林师妹 2024-04-18 9/450 2024-04-19 20:00 by nono2009
[基金申请] 基金和生小孩 +33 Ausy 2024-04-15 35/1750 2024-04-19 16:21 by feng6531
[教师之家] 美国教授开除中国研究生:我就不该录取你! +7 苏东坡二世 2024-04-15 7/350 2024-04-19 16:20 by feng6531
[考博] 25申博记录贴 +3 我属驴核动力驴 2024-04-18 4/200 2024-04-19 10:53 by 安塔瓦拉多
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +4 lin_2018 2024-04-17 4/200 2024-04-19 08:46 by TsolmonZ
[论文投稿] 投稿求助 5+3 我是洲洲啊 2024-04-17 5/250 2024-04-18 17:13 by topedit
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 4/200 2024-04-18 09:11 by 花花平头哥
[考研] 工学322,英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
[考研] 332求调剂 +3 木叶下1999 2024-04-16 5/250 2024-04-17 00:56 by angeliar
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[考研] 浙江海洋大学 船舶与海运学院 交通运输专硕 (交通信息工程及控制)接收调节 +4 joee 2024-04-15 8/400 2024-04-16 20:47 by TommyZiAng
[考研] 347求调剂 +3 寒辰ovo 2024-04-15 7/350 2024-04-16 19:05 by 寒辰ovo
[考研] 求调剂 +4 鹿萌月yy 2024-04-15 5/250 2024-04-16 16:16 by chenweiwade
[考研] 296求调剂 +3 Cclocomotive 2024-04-16 4/200 2024-04-16 10:04 by 19862091
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 328求调剂 +3 send rgsc 2024-04-14 7/350 2024-04-15 18:17 by zw_muchong
[考研] 300求调剂 +5 @Xqy668800 2024-04-14 6/300 2024-04-14 23:03 by lincunhui
信息提示
请填处理意见