[繪圖] 急!老師開大殺人了,nlinfit的初始值
各位版中高手大家好,
小弟前幾天才剛接觸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
推
12/07 12:29, , 1F
12/07 12:29, 1F
推
12/07 22:32, , 2F
12/07 22:32, 2F
→
12/08 11:45, , 3F
12/08 11:45, 3F
推
12/08 12:35, , 4F
12/08 12:35, 4F
→
12/08 12:35, , 5F
12/08 12:35, 5F
→
12/08 13:42, , 6F
12/08 13:42, 6F
→
12/08 21:00, , 7F
12/08 21:00, 7F
→
12/08 21:30, , 8F
12/08 21:30, 8F
→
12/08 21:32, , 9F
12/08 21:32, 9F
推
12/08 22:56, , 10F
12/08 22:56, 10F
推
12/09 00:16, , 11F
12/09 00:16, 11F
→
12/10 02:47, , 12F
12/10 02:47, 12F
推
12/10 08:39, , 13F
12/10 08:39, 13F
推
12/10 13:43, , 14F
12/10 13:43, 14F