24小时热门版块排行榜    

查看: 277  |  回复: 2
当前主题已经存档。

yujunhui

银虫 (著名写手)

[交流] 【求助】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 ]
回复此楼

» 猜你喜欢

安静
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)

请高手指点啊啊!
谢谢了!
安静
2楼2009-09-21 20:18:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)

3楼2009-09-27 15:51:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yujunhui 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见