最近在处理多b值MRI曲线拟合的数据,每组的数据结构如下
b:[0 20 50 80 100 150 200 400 600 800 1000];
S:[297 283.8 265.2 257.2 256.1 225.8 215.2 169.9 138.5 109.8 101.5];

DWI-MRI成像中b值和信号强度的比值关系如下:

  1. 单指数模型:S(b)/S(0) = exp(-b*ADC);
  2. 双指数模型:
  3. 拉伸指数:
  4. DKI

有了b值和S值,怎么拟合得到这些参数呢?利用MATLAB的曲线拟合工具箱可以很好的解决这些问题;

举个例子DKI拟合,对应的Matlab函数

ft = fittype( 'exp(-a*x+x^2*a^2*b/6)', 'independent', 'x', 'dependent', 'y' );
opts = fitoptions( ft );
opts.Algorithm = 'Levenberg-Marquardt';
opts.Display = 'Off';
opts.Lower = [-Inf -Inf];
opts.StartPoint = DKI_start;
opts.Upper = [Inf Inf];
[fitresult, gof] = fit( xData, yData, ft, opts );

[fitresult, gof]存储拟合的结果参数!但是曲线拟合工具箱对参数的初始值非常敏感,有没有一种拟合方法可以全局搜索出全局最优值?

http://cn.mathworks.com/help/gads/globaloptimsolution-class.html

其实题主也在试探性的研究这个问题,根据网上搜索的资料,初步可以认定Matlab最优化工具箱具有解决这个问题的潜在可能:
Matlab中有两个全局最长沙铁路段曹辉优化函数:MultistartGlobalSearch

fun = @(a,x) exp(-a*x);
problem = createOptimProblem('lsqcurvefit','objective', fun,'xdata',xdata,'ydata',ydata, 'x0',0.001);
ms = MultiStart;
% ms = GlobalSearch;
[a,fval,exitflag,output,solutions] = run(ms, problem,50);

MATLAB利用全局优化曲线拟合-段曹辉相关推荐

  1. matlab曲线拟合的应用,MATLAB软件在曲线拟合中的应用

    2007 年第 3 期 福 建 电 脑 MATLAB软件在曲线拟合中的应用 冯元珍 1, 屠小明 2, 罗建平 2 (1. 南京人口管理干部学院基础部 江苏 南京 210042 2. 南京医科大学数学 ...

  2. Matlab利用M_map和mapshow绘制网格地图

    Matlab利用M_map和mapshow绘制网格地图 最近在做利用网格数据绘制'Lambert'投影的地图,之前没有系统的利用matlab做过相关的工作,花了几天时间,问了专业的同学,最终整理了两种 ...

  3. Matlab利用膨胀腐蚀计算硬币数目

    Matlab利用膨胀腐蚀计算硬币数目 这段时间一直在学习Matlab图像处理,浅浅的记录一下一个简单的膨胀腐蚀算法. sample04.jpg 代码如下: %%计算图中硬币的数量I=imread('s ...

  4. matlab pup,matlab利用bar函数画不同颜色直方图

    matlab利用bar函数画直方图,参考文献[1]是matlab官方提供的help文档.里面提供了bar函数的基本用法,但是没有说明如何在同一张图中,为每个bar设置不同的颜色. 例子代码: myda ...

  5. Matlab 利用M文件产生模糊控制器

    Matlab 利用M文件产生模糊控制器 有位朋友咨询了一个问题,想要通过粒子群优化模糊控制的论域范围,在这种情况下,就无法直接通过matlab的fuzzy工具箱来自动生成fis模糊文件,需要通过M文件 ...

  6. ML之RF:基于Matlab利用RF算法实现根据乳腺肿瘤特征向量高精度(better)预测肿瘤的是恶性还是良性

    ML之RF:基于Matlab利用RF算法实现根据乳腺肿瘤特征向量高精度(better)预测肿瘤的是恶性还是良性 目录 输出结果 实现代码 输出结果 更新-- 实现代码 %RF:RF实现根据乳腺肿瘤特征 ...

  7. Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分、篮板、助攻

    Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分.篮板.助攻 目录 输出结果 预测勒布朗詹姆斯的20 ...

  8. matlab的exec程序,C++调用Matlab画图的一段程序

    劳动节闲来无事,写了一天程序,just for fun. 看,这是C++调用Matlab画图的一段程序.暂时不想多解释了,有兴趣的话,看看下面的代码吧. 以下几段代码由上到下,越来越旧.最上面的是最新 ...

  9. matlab求递归问题,matlab利用递归求解差分方程

    导航:网站首页 > matlab利用递归求解差分方程 时间:2018-5-23 matlab利用递归求解差分方程 function y = recur(a,b,n,x,x0,y0); % % y ...

最新文章

  1. 一个冷僻的知识点try直接返回finally里的设置null其实无效
  2. HTTP与服务器的四种交互方式
  3. WebApi项目创建CURD
  4. Android启动initlogo.rle制作
  5. gwt 同步和异步_GWT Spring和Hibernate进入数据网格世界
  6. python写算法求最短路径,Python实现迪杰斯特拉算法并生成最短路径的示例代码
  7. 【技巧】搜狗输入法特殊技巧
  8. NHibernate笔记Ⅲ--关于EF和NH的比较
  9. 【LeetCode】【数组】题号:*4238. 除自身以外数组的乘积
  10. 如何利用卡诺云系统管理早教机构?昆明收银系统还有此妙用!
  11. C++之struct
  12. 历年美赛赛题和特等奖论文集合,文中直接获取!!!
  13. 腾讯云/阿里云域名注册解析备案
  14. 如何选择Python版本2还是3
  15. Node js 开发之Postman
  16. 【C语言练习4】根据公式计算π的值,π=4-4/3+4/5-4/7+4/9+... 打印出一个表格来显示,用公式中的1项、2项、3项...计算出来π的近似值
  17. UCK区块链私享沙龙第二期圆满结束,感恩晚宴再聚共识
  18. C#如何读取身份证阅读器
  19. 计算机专业读书推荐卡图片大全,如何完成实践作业《好书推荐卡》,大家谈
  20. python图像运算、图像融合

热门文章

  1. java 输入提示_Java实现输入自动提示与补全功能
  2. 【读书笔记】-《程序员的思维锻炼-开发认知潜能的九堂课》
  3. day03:静态网页的练习
  4. [附源码]计算机毕业设计JAVA网易云音乐推荐系统
  5. box2dflash flash物理引擎
  6. [转载自 李国杰] 对计算机科学的反思
  7. 2022年我的22个感悟
  8. Android OpenGL ES FrameBuffer离屏渲染
  9. android 以太网工具,Android使用ethtool工具查询/设置以太网网卡状态
  10. java中使用tika,Apache Tika使用报告