版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

OlivarD

新蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 712
帖子: 5
在線: 30.1小時
蟲號: 1854391
注冊: 2012-06-10
專業(yè): 生物信息學(xué)

[求助] GM（1,1）的matlab程序

1. 這個程序只能運行一次，第二次就會出現(xiàn)維度不對的問題，改了很多遍也都這樣，不知道怎么改才能讓他正常運行
2. GM(1,1)的a，b這里沒有輸出，本人菜鳥一枚，不知道應(yīng)該怎么編寫這個程序，還望大神指點一二
3. 以上每一步計算出來的數(shù)據(jù)如果輸出保存
很急，謝謝各位了~

h=xlsread('Microsoft.xlsx');%數(shù)據(jù)保存在excel里面調(diào)用
global k;
for k=1:1:1254 %一行一行的計算
y=[h(k,: )];
n=length(y);
yy=ones(n,1);
yy(1)=y(1);
clear
for i=2:n
yy(i)=yy(i-1)+y(i);
end
B=ones(n-1,2);
for i=1: (n-1)
B(i,1)=-(yy(i)+yy(i+1))/2;
B(i,2)=1;
end
BT=B';
for j=1:n-1
YN(j)=y(j+1);
end
YN=YN';
A=inv(BT*B)*BT*YN; %出錯的地方
a=A(1);
u=A(2);
t=u/a;
t_test=1;
i=1:t_test+n;
yys(i+1)=(y(1)-t).*exp(-a.*i)+t;
yys(1)=y(1);
for j=n+t_test:-1:2
ys(j)=yys(j)-yys(j-1);
end
x=1:n;
xs=2:n+t_test;
yn=ys(2:n+t_test);
det=0;
for i=2:n
det=det+abs(yn(i)-y(i));
end
det=det/(n-1);
disp(['百分絕對誤差為：',num2str(det),'%']);
disp(['預(yù)測值為： ',num2str(ys(n+1:n+t_test))]);
end

[ Last edited by OlivarD on 2013-5-12 at 01:41 ]

回復(fù)此樓

» 猜你喜歡

329求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何求這個matlab程序畫出的legend曲線的法矢量？已經(jīng)有6人回復(fù)
急求，求幫忙寫一個matlab程序！已經(jīng)有3人回復(fù)
最優(yōu)化方法理論的MATLAB程序典型例子已經(jīng)有158人回復(fù)
想寫matlab程序之前不妨看看matlab的fileexchange 已經(jīng)有4人回復(fù)
算法程序用matlab，C語言，軟件平臺，硬件平臺運行耗時之間的關(guān)系已經(jīng)有9人回復(fù)
求解相算法的相關(guān)matlab程序塊已經(jīng)有3人回復(fù)
求窗口傅里葉變換的matlab程序已經(jīng)有6人回復(fù)
關(guān)于振動學(xué)方面的MATLAB程序怎么發(fā)核心？已經(jīng)有10人回復(fù)
中心差分求解二維熱傳導(dǎo)的matlab程序已經(jīng)有12人回復(fù)
如何知道一個MATLAB程序在計算中占的內(nèi)存是多大？謝謝！已經(jīng)有3人回復(fù)
matlab程序中調(diào)用m文件的問題已經(jīng)有13人回復(fù)
求助多輸出支持向量機的matlab程序已經(jīng)有6人回復(fù)
模式搜索MATLAB程序已經(jīng)有9人回復(fù)
matlab程序問題，總是出錯，麻煩大家看看問題在那里？已經(jīng)有6人回復(fù)
William T. Freeman的基于例子超分辨率matlab程序跑不起來��！已經(jīng)有6人回復(fù)
關(guān)于貝葉斯線性判別分析的matlab程序已經(jīng)有4人回復(fù)
同樣的MATLAB程序，為何在別人的電腦上運行不了已經(jīng)有21人回復(fù)
matlab 串口通信的問題已經(jīng)有3人回復(fù)
【求助】一組光譜數(shù)據(jù)的多元散射校正（MSC），求具體的MATLAB程序及用法已經(jīng)有6人回復(fù)
【求助】多目標pso算法的MATLAB程序已經(jīng)有3人回復(fù)

1樓 2013-05-11 23:22:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時
蟲號: 1530793
注冊: 2011-12-09
性別: GG
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

★ ★
1592203609: 金幣+2, 應(yīng)助指數(shù)+1, 謝謝參與 2013-05-16 12:05:58

1.把數(shù)據(jù)發(fā)上來，給你試試，如果不方便，就把數(shù)據(jù)減少。發(fā)上來。
至少把出錯信息給出。

2.如果沒錯的話，a,u就是你想要的a,b.
a=A(1);
u=A(2);

3.用save 命令，存成mat文件，保存你的變量。你沒說明白每一步你要保存什么。

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

2樓2013-05-15 14:59:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 OlivarD 的主題更新