24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yujunhui

银虫 (著名写手)

应助: 22 (小学生)
金币: 312.9
散金: 447
红花: 12
帖子: 1218
在线: 111.2小时
虫号: 669946
注册: 2008-12-07
性别: GG
专业: 金属结构材料

[交流] 【求助】postreg函数运行出现错误,请指点？

postreg函数运行出现错误,请指点？
源程序

clear all
clc

s = [             0.0758                0.3030    0.3535 0.3788    0.4545                0.8081    1.0354 1.0606    1.1111 1.1364       1.1616 1.1869    1.2121          1.2879 1.3384    1.5152 1.5657 1.6919    1.7424       1.9697    2.0202             2.2222    2.3990 2.4242 2.4495       2.5000       2.8283    2.8788 2.9293       2.9545       3.0051             3.2576       3.2828       3.3081    3.3333       3.6364          3.6869                3.7121    3.7374    3.7626    3.7879    3.8131       3.8636    3.8889          3.9141    3.9394       3.9646    3.9899       4.0152    4.0404    4.0657          4.0909             4.1162       4.1414          4.2424 4.4192 4.4444 4.4697 4.4949 4.5202    4.5455    4.7980          4.8232       4.8485 4.8737 4.8990       4.9242    4.9495       4.9747    5.0000       5.0253 5.0505    5.0758 5.1010    5.1263    5.1515    5.1768    5.2020       5.2273 5.4545       5.4798    5.5051       5.6566       5.6818 5.7323    5.9091 5.9343    6.0354 6.0606 6.0859 6.1111       6.1364 6.1616    7.0149       7.2139          7.6866       7.6968       7.7861];
p1=  1.0e+003 *[             0.5480             0.6948    0.7086 0.7090    0.7101             0.7130    0.7131 0.7136    0.7143    0.7143       0.7149 0.7155    0.7159       0.7168 0.7171 0.7175 0.7178 0.7181    0.7183       0.7193    0.7195                0.7217    0.7304 0.7313    0.7375       0.7381       0.7383    0.7422    0.7444    0.7450       0.7467             0.7590       0.7604       0.76153 0.76798       0.7683          0.7698             0.7703       0.7710    0.7729    0.7747    0.7763       0.7777    0.7090          0.7101    0.7111       0.7120       0.7129       0.7120    0.7129       0.7136       0.7143          0.7145       0.7149          0.7168 0.7171 0.7175 0.7178 0.7181 0.7183       0.7193    0.7195       0.7198       0.7202 0.7206 0.7217       0.7304    0.7313       0.7375    0.7400       0.7428 0.7460    0.7492    0.7525    0.7584    0.7609    0.7669       0.7676       0.7682 0.7687       0.7690       0.7692       0.7694    0.7695    0.7698       0.7699 0.7699 0.7700 0.7701 0.7703    0.7705       0.7708 0.7712    0.7719       0.7729          0.8041    0.8111          0.8643 ]
p2=          [          634.4721             750.0815 752.1115 755.1815 755.2815             757.0006    757.0010 757.0015 757.0025  757.0045    757.0105  757.0115    757.0215    757.0315  757.0415 757.0425  757.0435  757.0445 757.0455    757.0515    757.0600          757.0715 757.0725  757.0730 757.0735    757.0740    757.0757    757.0762  757.0771    757.0785    757.0795             757.0836    757.0844    758.0615    759.0001    759.0005       759.0015          759.0115    759.0321 759.0351    759.0415    759.0617    759.0815    759.1231    759.1285    759.2115    760.0135    760.0144       760.0151 761.2015    762.215       763.1159       764.8416    765.8436          768.8405  769.8406  770.8416 768.8411  769.8433  769.8506 771.7391    772.2567       772.2567    772.2567  772.2567  772.2567    771.7391 768.1159       771.7391 771.7391    771.7391 772.2567 772.2567    772.2567 772.2567 773.7391 778.7391    780.7391    784.4721 785.4721    787.4721    788.1965    789.1965    790.1965 792.1965    793.1965 798.1965  800.1965  800.19715  801.014  803.0426    804.6228 811.0252 826.1824    826.1538       833.5285 833.7310    861.1340 ]
p3= 1.0e+002 *  [             5.8032             6.4934    6.5313    6.6716    6.8229             7.0804    7.0885 7.0954    7.1011 7.1056       7.1089 7.1113       7.1129       7.1139 7.1138 7.1139 7.1233 7.1311    7.1341       7.1411    7.1513                7.2107    7.2114 7.2122    7.2129       7.2535       7.2754    7.3001 7.3149    7.3252       7.3369             7.4645       7.5157    7.5160       7.5161    7.5164          7.5166             7.5167    7.5170    7.5171       7.5260    7.5376       7.5479    7.5481       7.5483    7.5485       7.5561       7.5669       7.5725    7.5733    7.5738       7.5750          7.5793       7.5789             7.5843    7.5848 7.5850 7.5860    7.5862 7.5873    7.5880       7.5881       7.5899       7.5925    7.5928    7.5935    7.5940    7.5969       7.5970    7.5985    7.5986 7.5992    7.6017    7.6019    7.6039    7.6063    7.6138       7.6162    7.6380    7.6552       7.6675    7.6753       7.6798       7.6822 7.6824       7.6832 7.6834    7.6838 7.6839 7.6840 7.6842       7.6845    7.6854 7.6863       7.6865          7.7527       7.7647       7.791]

p=s;
t=[p1;p2;p3];

[pn,meanp,stdp,tn,meant,stdt] = prestd(p,t);
[ptrans,transMat] = prepca(pn,0.02);
net = newff(minmax(ptrans),[15 3],{'tansig','purelin'},'trainlm');
net.trainParam.show=200;
net.trainParam.lr=0.05;
net.trainParam.epochs=30000;
net.trainParam.goal=1e-3;
net = init(net);

net = train(net,ptrans,tn);
an = sim(net,ptrans);
a = poststd(an,meant,stdt);
[m,b,r] = postreg(a,t);

运行出现

Warning: Rank deficient, rank = 4,  tol = 1.6593e-010.
> In postreg at 57
??? Error using ==> unknown
Matrix dimensions must agree.
Error in ==> postreg at 72
an = a - mean(a);

请高手指点？
谢谢！

[ Last edited by nono2009 on 2009-9-21 at 22:23 ]

回复此楼