| 查看: 707 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
s084730木虫 (正式写手)
|
[求助]
【求助】matlab中的几个公式含义
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% To fit a 2-D gaussian %% m = Image %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [cx,cy,sx,sy,PeakOD] = Gaussian2D(m,tol); [sizey sizex] = size(m); [x,y] = MeshGrid(1:sizex,1:sizey); fit = abs(PeakOD)*(exp(-0.5*(x-cx).^2./(sx^2)-0.5*(y-cy).^2./(sy^2))); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% a function to fit a thermal cloud 2-D function [cx,cy,sx,sy,PeakOD] = Gaussian2D(m,tol); %% m = image %% tol = fitting tolerance options = optimset('Display','off','TolFun',tol,'LargeScale','off'); [sizey sizex] = size(m); [cx,cy,sx,sy] = centerofmass(m); pOD = max(max(m)); mx = m(round(cy),  ;x1D = 1:sizex; ip1D = [cx,sx,pOD]; fp1D = fminunc(@fitGaussian1D,ip1D,options,mx,x1D); cx = fp1D(1); sx = fp1D(2); PeakOD = fp1D(3); my = m(:,round(cx))'; y1D = 1:sizey; ip1D = [cy,sy,pOD]; fp1D = fminunc(@fitGaussian1D,ip1D,options,my,y1D); cy = fp1D(1); sy = fp1D(2); PeakOD = fp1D(3); [X,Y] = meshgrid(1:sizex,1:sizey); initpar = [cx,cy,sx,sy,PeakOD]; fp = fminunc(@fitGaussian2D,initpar,options,m,X,Y); cx = fp(1); cy = fp(2); sx = fp(3); sy = fp(4); PeakOD = fp(5); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % PURPOSE: find c of m of distribution function [cx,cy,sx,sy] = centerofmass(m); [sizey sizex] = size(m); vx = sum(m); vy = sum(m'); vx = vx.*(vx>0); vy = vy.*(vy>0); x = [1:sizex]; y = [1:sizey]; cx = sum(vx.*x)/sum(vx); cy = sum(vy.*y)/sum(vy); sx = sqrt(sum(vx.*(abs(x-cx).^2))/sum(vx)); sy = sqrt(sum(vy.*(abs(y-cy).^2))/sum(vy)); 其中cx = sum(vx.*x)/sum(vx); cy = sum(vy.*y)/sum(vy); sx = sqrt(sum(vx.*(abs(x-cx).^2))/sum(vx)); sy = sqrt(sum(vy.*(abs(y-cy).^2))/sum(vy));是什么含义 |
回复此楼» 猜你喜欢
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 技术路线里面的公式太复杂，符号太多，需要每个都解释什么含义么？
已经有8人回复
- MATLAB中有意思的why函数
已经有5人回复
- 请教公式中大括号的含义
已经有6人回复
- 理解数字信号处理的三把钥匙
已经有37人回复
- 【求助】MATLAB中Normalization的数据预处理方法具体是指什么意思啊
已经有12人回复
s084730
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4431.5
- 散金: 644
- 红花: 1
- 帖子: 534
- 在线: 130.1小时
- 虫号: 774019
- 注册: 2009-05-19
- 性别: GG
- 专业: 光学信息获取与处理
3楼2013-04-30 11:19:28
2楼2013-04-27 23:32:24
