Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

jilly120 发表于 2005-12-8 10:54

哪位看到的大虾帮帮忙啊？谢谢！

哪位大虾可以帮小女子看看下面的程序呀？
p=
t=; 
pp=; 
tt = ; 
LearnRate = 0.01; LearnGoal = 0.01; LearnEpoch = 4000; HiddenLayer = 3; OutputLayer = 1;
 = fuc_training(p,t,pp,tt,LearnRate,LearnGoal,LearnEpoch,HiddenLayer,OutputLayer);
a=0; =SaveDlg(a); switch Save_or_not case 'Yes', save qq p t net TR R2 pred pp tt;%DemoResult is the filename to be saved otherwise end
function =fuc(p,t,pp,tt,speed,goal,epoch,hiden_layer,output_layer)
%preprocessing the input vectors = prestd(p); = prepca(p,0.02); =premnmx(p,t); pp = trastd(pp,meanp,stdp); pp = trapca(pp,transMat); pp =tramnmx(pp,minp,maxp); Pt = minmax(p); net=newff(Pt,,{'tansig','purelin'},'trainlm'); fprintf('\n\n'); echo off ; net.trainParam.lr=speed; net.trainParam.epochs=epoch; net.trainParam.goal=goal; net.trainParam.min_grad=1e-15; net.trainParam.show=500; %net.trainParam.max_fail=100; %net.trainParam.mem_reduc=10; = train(net,p,t); pred = sim(net,pp); pred = postmnmx(pred,mint,maxt); cv = (pred-tt)./tt; R = corrcoef(pred,tt); %R2 = R(1,2); % = size(tt);
fprintf('\nthe predicted value is '); fprintf('%12.4f',pred); fprintf('\nthe measured value is'); fprintf('%12.4f',tt); fprintf('\nthe error is '); fprintf('%12.4f',cv); %fprintf('\nthe correlate coefficient is %6.4f', R2);
 pic(tt,pred,R);
 function =Save_Dlg(a) s=questdlg('是否保存训练结果?','保存训练结果','Yes','No','No');

我想把里面的P矩阵换成下面的数据：
p=;
t=; 
pp=; 
tt = ; 

明明已经换了，可为什么出来的结果还是预测的原来的数据的值呢？？？是不是要建立一个叫做“MATLAB data file"的文件呢？这个文件在哪里创建呢？
还有一个问题就是我只有四个样本，是不是对于神经网络预测来说太少了呢？
哪位高手明白，请帮帮忙吧！谢谢呀！！
或者直接回复到我的邮箱<a href="mailtjjpurple120@tom.com" target="_blank" >jjpurple120@tom.com</A>也可以，谢谢！！

suffer 发表于 2005-12-8 20:57

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！

明明已经换了，可为什么出来的结果还是预测的原来的数据的值呢？？？
这句话不是很明白，能否说清出点
4组样本确实很少
另外你给的程序中少了pic这段代码，运行不了

[ 本帖最后由 ChaChing 于 2009-7-11 23:27 编辑 ]

jilly120 发表于 2005-12-8 21:39

多谢suffer的提醒！ pic的代码如下： function a=fun(a,b,R) hold off; =size(a); for i=1:col plot(a(i),b(i),':ro'); hold on; end temp1=min(a); temp2=min(b); x=zeros(1,2); x(1)=min(temp1,temp2); temp1=max(a); temp2=max(b); x(2)=max(temp1,temp2); y=x; plot(x,y); hold off; title('The Prediction'); xlabel('Meas Data'); ylabel('Pred Data'); hold off; 可否帮我看看整个程序是否有问题呢？因为我运行时只训练1epoch，怎么就停止了呢？

suffer 发表于 2005-12-9 15:46

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！

我运行上述程序的结果 TRAINLM, Epoch 0/4000, MSE 0.304479/0.01, Gradient 11.2279/1e-015 TRAINLM, Epoch 500/4000, MSE 0.018452/0.01, Gradient 0.00154877/1e-015 TRAINLM, Epoch 1000/4000, MSE 0.0184515/0.01, Gradient 6.7753e-005/1e-015 TRAINLM, Epoch 1500/4000, MSE 0.0184514/0.01, Gradient 0.000540231/1e-015 TRAINLM, Epoch 2000/4000, MSE 0.0184513/0.01, Gradient 4.81969e-005/1e-015 TRAINLM, Epoch 2500/4000, MSE 0.0184513/0.01, Gradient 0.000500366/1e-015 TRAINLM, Epoch 3000/4000, MSE 0.0184512/0.01, Gradient 6.42612e-005/1e-015 TRAINLM, Epoch 3500/4000, MSE 0.0184512/0.01, Gradient 4.57205e-005/1e-015 TRAINLM, Epoch 4000/4000, MSE 0.0184512/0.01, Gradient 4.13985e-005/1e-015 TRAINLM, Maximum epoch reached, performance goal was not met.
 the predicted value is 143.0908 80.4107 121.4202 156.2426 75.6621 the measured value is 247.4000 283.4000 199.2000 99.3000 286.4000 the error is -0.4216 -0.7163 -0.3905 0.5734 -0.7358

jilly120 发表于 2005-12-10 15:41

呵呵，这正是我想问的问题！ 我把程序开始的P矩阵换成了： 
p=;
t=; 
pp=; 
tt = ; 在运行时只训练1epoch，就停止了 结果如下： the predicted value is 4.7 4.5 3.5 4.2 5.2 3.7 the measuredvalue is 6.9 6.2 6.2 4.3 4.9 5.2 误差很大，而且只运行了一代！

happy 发表于 2005-12-10 18:40

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！

t明显和上面的不对应阿

页: [1]

声振论坛's Archiver

哪位看到的大虾帮帮忙啊？谢谢！

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！

回复：（jilly120）哪位看到的大虾帮帮忙啊？谢谢！