| 查看: 732 | 回复: 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));是什么含义 |
回复此楼» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 技术路线里面的公式太复杂，符号太多，需要每个都解释什么含义么？
已经有8人回复
- MATLAB中有意思的why函数
已经有5人回复
- 请教公式中大括号的含义
已经有6人回复
- 理解数字信号处理的三把钥匙
已经有37人回复
- 【求助】MATLAB中Normalization的数据预处理方法具体是指什么意思啊
已经有12人回复
2楼2013-04-27 23:32:24
s084730
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4431.5
- 散金: 644
- 红花: 1
- 帖子: 534
- 在线: 130.1小时
- 虫号: 774019
- 注册: 2009-05-19
- 性别: GG
- 专业: 光学信息获取与处理
3楼2013-04-30 11:19:28
100