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

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

24小時(shí)熱門版塊排行榜

返回列表

OlivarD

新蟲 (初入文壇)

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

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

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

h=xlsread('Microsoft.xlsx');%數(shù)據(jù)保存在excel里面調(diào)用
global k;
for k=1:1:1254 %一行一行的計(jì)算
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; %出錯(cuò)的地方
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(['百分絕對(duì)誤差為：',num2str(det),'%']);
disp(['預(yù)測(cè)值為： ',num2str(ys(n+1:n+t_test))]);
end

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

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

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

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

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

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

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

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

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

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

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