作者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