| 查看: 1663 | 回复: 6 | |||
mathislhc新虫 (正式写手)
|
[求助]
求大神告知Matlab最小二乘法拟合分布函数 已有2人参与
|
|
看一篇文献 其中，腐蚀损伤数据的常见分布函数有：对数正态、极值I型Gumbel分布和三参数威布尔分布； 现在只知道以上三种分布函数F(x)和腐蚀深度x 然后通过最小二乘法可以估计参数（位置参数、尺度参数和形状参数） 但小白不知道F(x)是多少啊？？ Matlab 怎么拟合呢？？？ 直接用webfit 拟合腐蚀深度，估计的参数也和文献不一样啊... 求指导！！ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有17人回复
-
为什么nbs上溴 没有产物点出现呢
已经有6人回复
-
求推荐博导
已经有4人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
FMStation
至尊木虫 (知名作家)
- 仿真EPI: 1
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18189.7
- 红花: 97
- 帖子: 8891
- 在线: 1487.4小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
【答案】应助回帖
感谢参与，应助指数 +1
|
https://www.mathworks.com/exampl ... ng-the-fit-function Fit Exponential Models Using the fit Function https://www.mathworks.com/examples/curvefitting Curve Fitting Toolbox |
2楼2016-11-16 10:39:08
mathislhc
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3838.5
- 红花: 3
- 帖子: 311
- 在线: 337.8小时
- 虫号: 4034713
- 注册: 2015-08-23
- 性别: GG
- 专业: 机械结构强度学
3楼2016-11-16 12:14:23
mathislhc
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3838.5
- 红花: 3
- 帖子: 311
- 在线: 337.8小时
- 虫号: 4034713
- 注册: 2015-08-23
- 性别: GG
- 专业: 机械结构强度学
4楼2016-11-16 19:48:17
mathislhc
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3838.5
- 红花: 3
- 帖子: 311
- 在线: 337.8小时
- 虫号: 4034713
- 注册: 2015-08-23
- 性别: GG
- 专业: 机械结构强度学
5楼2016-11-18 19:35:35
huab1984666
新虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 241.5
- 散金: 2094
- 红花: 12
- 帖子: 1320
- 在线: 210.3小时
- 虫号: 516828
- 注册: 2008-03-03
- 性别: GG
- 专业: 胶体与界面化学
【答案】应助回帖
|
以对正态分布为例 以一种情况：损伤系数X与一个自标量（例如温度T） 目标函数：Mifit.m function f=Myfit(beta) data=load('XX.txt'); T=data(:,1); x=data(:,2); %% mu=beta(1); sigma=beta(2) FX=lognpdf(T,mu,sigma) PX=(FX-x).^2 MT= sum(PX); %% F-test检验参数 f=MT./(length(T)-2); %2为未知参数的个数； 优化拟合工具：以模式搜索法patternsearch为例 lb=[0 0]; %未知参数mu和sigma最小估计 ub=[inf inf];%未知参数mu和sigma最大估计 options = psoptimset('Display','Iter','MaxIter',500,'TolFun',1e-6,'TolX',1e-6,'CompleteSearch','on'); beta0=[0.76]; [beta,fval,exitflag,output] = patternsearch(@Mifit,beta0,[],[],[],[],lb,ub,options ); beta0=beta; [beta,fval,exitflag,output] = patternsearch(@Mifit,beta0,[],[],[],[],lb,ub,options ); %% mu=beta(1); sigma=beta(2); %% data=load('XX.txt'); T=data(:,1); x=data(:,2); FX=lognpdf(T,mu,sigma) plot(T,data,'r.',T,FX,'k-'); 第二种情况：之检验损伤系数X的最佳分布；X的数据做图为Y轴 以log-normal分布为例 你用“”help lognfit“” 就知道怎么用了 parmhat = lognfit(data) [parmhat,parmci] = lognfit(data) [parmhat,parmci] = lognfit(data,alpha) [...] = lognfit(data,alpha,censoring) [...] = lognfit(data,alpha,censoring,freq) [...] = lognfit(data,alpha,censoring,freq,options)； 其中data就是你的X值，alpha就是代表置信的α值 |
6楼2016-11-21 09:28:21
huab1984666
新虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 241.5
- 散金: 2094
- 红花: 12
- 帖子: 1320
- 在线: 210.3小时
- 虫号: 516828
- 注册: 2008-03-03
- 性别: GG
- 专业: 胶体与界面化学
7楼2016-11-21 09:30:07