作者johnaqu (Dade文)
看板MATLAB
標題[繪圖] 急!老師開大殺人了,nlinfit的初始值
時間Fri Dec 7 12:20:57 2012
各位版中高手大家好,
小弟前幾天才剛接觸matlab這好玩的東西
現在困惑於沒辦法做出令老師滿意的平面圖,故來此求助各位高手
function mx=Mx(tau,mO,t1); %定義function
Mo=input('hi,Mo?');
T1=input('T1=');
tau=(1:1000); %1000組數據
y=Mo*(1-2*exp(-tau/T1));
nois=rand(1,1000);cnois=nois*1; %提高或增強noise強度
Y=y+cnois; %noise曲線圖
hold on;
ft=inline('beta(1)+2*beta(1)*exp(-tau/beta(2))','beta','tau'); %自變函數
beta=nlinfit(tau,Y,ft,[Mo T1]);
Tau=1:1000;
mO=beta(1);t1=beta(2);
mx=mO*(1-2*exp(-Tau/t1));
plot(Tau,mx,['r'],tau,Y,['b']);
title('T1 relaxation');
xlabel('Time')
ylabel('Magnetic vector')
end
小弟的作業要做兩張圖，tau vs. Y 那張圖沒問題，現在的問題是如何用tau vs. Y
的這張圖來nlinfit出像似的曲線
在本版有找到利用迴圈的方式找出適合的初始值，問題是小弟實在不懂該怎麼將
迴圈的語法與nlinfit做結合?
不知道可否請任何高手大大給予幫助，小弟感激不盡0..0
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.120.19.1
1F：推 crc2121:就讓老師殺吧XDDDDD 今年沒畢業,了不起明年當學第XDDDD 12/07 12:29
2F：推 uu4avzd:樓上回文好屌 12/07 22:32
3F：→ johnaqu:屌是屌 可是小弟很無奈ˊˋ..唉 12/08 11:45
4F：推 alamabarry:建議你把問題簡化~~再來問 12/08 12:35
5F：→ alamabarry:而不是整個貼上來 12/08 12:35
6F：→ johnaqu:嗯謝謝 可是我的問題就是怎麼利用迴圈來找初始值0..0 12/08 13:42
7F：→ panpan19:看起來像是用某種演算法，去除雜訊，以得到主要訊號? 12/08 21:00
8F：→ johnaqu:hi大大 是故意製造有雜訊的資料,然後非線性迴歸出原本的 12/08 21:30
9F：→ johnaqu:理論曲線,可是小弟我不會設定最佳初始值ˊˋ 12/08 21:32
10F：推 crc2121:衝吧~~ 12/08 22:56
11F：推 crc2121:我是你就會"退"選了XDDDDD老師擺明要砍人還菇在裡面@@"a 12/09 00:16
12F：→ tomichy:12月還可以退選 ?? 12/10 02:47
13F：推 kurt28:有些可以 12/10 08:39
14F：推 crc2121: 停修 12/10 13:43