EssayproSolutions

(Solved) : Gaussian Fitting Matlab Trying Gaussian Fitting Using Matlab Want Obtain Parameters Plot F Q30285653 . . .

gaussian fitting for Matlab/

I am trying gaussian fitting using matlab.
You want to obtain the parameters and plot the fitted gaussianfunction calculated over the original signal.
However, in this process, one of the parameters, mu2, is fixed to999. How do I fix the code?

clear all; close all; clc;
%% load DATA
load(‘SPR0056_base_waves_allPEAKS.mat’);
%% normalization, interpolation
cnt=51;

PPG=DATA.WAVE(:,2);
PT=DATA.LPEAK(:,1);
PTI=round(PT*300);
L=length(PTI);

x=PTI(cnt):PTI(cnt+1);
y=PPG(x);
x = x-PTI(cnt);

norm_y=((y-min(y)) / (max(y)-min(y)));
re_PPG=interpol3(x,norm_y,x(end));
n=[1:1:1000];

plot(re_PPG,’color’,[0.1 0.1 0.1],’LineWidth’,1);
title([‘segment ‘,num2str(cnt)])
xlabel(‘Number of samples’)
ylabel(‘Norm a.u.’)
axis tight

%p= segGaussFit(n, re_PPG, 2, 1)

%% Gaussian fitting

miny = min(re_PPG); % y offset
dn=length(n);
flag =0;
mu=[];
for step = 1:(dn-1)
for initmu = 0:dn
flag = flag +1;
stg1 = double(initmu+1);
stg2 = double(stg1 + step);
  
if (stg1 + step >= dn) || (stg2 >= dn)
break;
end
  
temp1 = find(re_PPG(1:stg2-1)>0);
  
if temp1 > 0
p1 = max(re_PPG(temp1)); %postive value   
mu1 = find(re_PPG==p1);
else
p1 = re_PPG(abs(re_PPG)==max(abs(re_PPG(1:stg2-1)))); % initalguess height
mu1 = find(re_PPG==p1);
end
  
p2 = std(re_PPG(1:stg2-1)); % inital guess
miny1 = min(re_PPG(1:stg2-1)); % y offset
  
  
temp2 = find(re_PPG(stg2:dn)>0);
  
if temp2 > 0
temp2 = temp2 + stg2-1;
p3 = max(re_PPG(temp2));
mu2 = temp2(find(re_PPG(temp2)==p3));
else
p3 = re_PPG(abs(re_PPG)==max(abs(re_PPG(stg2:dn)))); % inital guessheight
mu2 = find(re_PPG==p3);
end
  
p4 = std(re_PPG(stg2:dn)); % inital guess
miny2 = min(re_PPG(stg2:dn)); % y offset
  
p5 = mu1;
p6 = mu2;

gaumd1 = (p1.* exp(-(((n-p5).^2)/(2*p2^2)))+miny1)+ …
(p3.* exp(-(((n-p6).^2)/(2*p4^2)))+miny2);

p = [p1 p2 p3 p4 p5 p6];
end
end

clear al close al l clc; load DATA load( SPR0056 base-waves_al IPEAKS.mat) normal i zation, interpolation cnt 51 PPG-DATA. WAVE(,2) PT-DATA. LPEAK( : , 1 ) ; PT I round (PT+300); L-length(PTI) -PTI (cnt) :PTI (cnt +1 ) ; y-PPG(x) x = x-PTI (cnt ) ; norm_y ((y-min(y)) / (max(y)-min(y))); re-PPG-interpol 3(x, norm-y, x( end)) ; n [11:1000] plot (re_PPG, color,[0.1 0.1 0.11, LineWidth,1); title(I segment num2str(cnt)]) xlabel Number of samples’) ylabel Norm a.u.’) axis tight Show transcribed image text

Expert Answer

---

Answer to Gaussian Fitting Matlab Trying Gaussian Fitting Using Matlab Want Obtain Parameters Plot F Q30285653 . . .