tnjk.net
当前位置:首页 >> mAtlAB高斯拟 >>

mAtlAB高斯拟

新建一个函数,代码如下: function [fitresult, gof] = createFit(x, y) [xData, yData] = prepareCurveData( x, y); % Set up fittype and options. ft = fittype( 'gauss1' ); opts = fitoptions( 'Method', 'NonlinearLeastSquares' ); opts...

你这个题很简单啦。不用使用cholesky分解就可以直接做了~ 因为你给的两个协方差矩阵都是对角线矩阵,这意味着两次二维变量的各维之间是不相关的,对于高斯过程也就是相互独立的,那么你就可以直接利用randn函数来生成数据。对角线上的值就是各自...

使用matlab中的拟合工具箱吧,很简单的。 在command windows下输入cftool进入拟合工具箱,Data中输入你的数据,然后使用Fitting进行拟合,里面可以选择Gaussian,这样可以得出完整的高斯函数表达式。

在拟合前可以先对x数据进行标准化: %% x,y是题目中的数据plot(x, y, 'or')hold on[t, mu, sigma] = zscore(x);cfun = fit(t(:), y(:), fittype('a * exp(-((x - b)/c)^2)'));xx = 1:0.1:60;t = (xx-mu)/sigma;yy = cfun.a * exp(-((t - cfun.b)...

就是要多峰值的高斯函数吧?下面的这个行不行? ezplot('0.5*exp(-(x+4)^2)+exp(-(x/2)^2)+0.7*exp(-(x-5)^2)')

你是要进行数值拟合吧?要有数据的,请把数据、拟合样式写下。

假设是32点的高斯窗,用以下matlab语句就可以看到高斯窗时域和频域的情况: N = 32; wvtool(gausswin(N)); 或者: window=gausswin(n); [h,w]=freqz(window,1); subplot(1,2,1) stem(window); subplot(1,2,2); plot(w/pi,20*log(abs(h)/abs(h(1...

1. WGN:产生高斯白噪声 y = wgn(m,n,p) 产生一个m行n列的高斯白噪声的矩阵,p以dBW为单位指定输出噪声的强度。 y = wgn(m,n,p,imp) 以欧姆(Ohm)为单位指定负载阻抗。 y = wgn(m,n,p,imp,state) 重置RANDN的状态。 在数值变量后还可附加一些标志...

关于MATLAB曲线拟合,我写了一系列的经验,为了相互统一,采用下面的数据: x=[0 0.3000 0.6000 0.9000 1.2000 1.5000 1.8000 2.1000 2.4000 2.7000 3.0000] y=[2.0000 2.3780 3.9440 7.3460 13.2320 22.2500 35.0480 52.2740 74.5760 102.6020 ...

可用:lsqcurvefit or nlinfit 进行拟合。格式为:b=lsqcurvefit(fx,b,x,y);b=nlinfit(x,y,fx,b);其中圆括号内的fx为拟合函数,以你的问题,可表示为:fx=@(b,x)b(1)*exp(-((x-b(2))./b(3)).^2);圆括号中的b为初值,其中的b(2)最好能给出峰值对...

网站首页 | 网站地图
All rights reserved Powered by www.tnjk.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com