| 查看: 1770 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
mathislhc新虫 (正式写手)
|
[求助]
求大神告知Matlab最小二乘法拟合分布函数 已有2人参与
|
||
|
看一篇文献 其中，腐蚀损伤数据的常见分布函数有：对数正态、极值I型Gumbel分布和三参数威布尔分布； 现在只知道以上三种分布函数F(x)和腐蚀深度x 然后通过最小二乘法可以估计参数（位置参数、尺度参数和形状参数） 但小白不知道F(x)是多少啊？？ Matlab 怎么拟合呢？？？ 直接用webfit 拟合腐蚀深度，估计的参数也和文献不一样啊... 求指导！！ |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有43人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
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
FMStation
至尊木虫 (知名作家)
- 仿真EPI: 1
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.2小时
- 虫号: 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
