Добавил:
debilX2
t.me
Фулл всегда есть, ФОЭ ТОЭ ТЭЦ Электроника, КЭТ ИиКГ и тд https://t.me/whitedevil752
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:lab2 вроде / func_weibull_hist_visual
.m function[g] = func_weibull_hist_visual(Y,b,a,K,nbins)
% check Y - uint8
[counts,xn] = imhist(Y,nbins); % n - counts; xn - binLocations;
figure(1), stem(xn,counts), grid;
y1 = counts/(numel(Y)*(xn(2)-xn(1))); % считаем плотность по гистограмме
t = linspace(1,256,1000);
y2 = pdf('Weibull',t,b,a);
% y2= gampdf(t,K,b);
figure(2), plot(t,y2,xn,y1,'r'), grid
legend('Теорет. Вейбулл-плотность','Гистограмма полученных чисел')
xlabel('x'), ylabel('Histogram');
title(['b = ',num2str(b),' a = ',num2str(a),' K = ',num2str(K)]);
% check Y - uint8
[counts,xn] = imhist(Y,nbins); % n - counts; xn - binLocations;
figure(1), stem(xn,counts), grid;
y1 = counts/(numel(Y)*(xn(2)-xn(1))); % считаем плотность по гистограмме
t = linspace(1,256,1000);
y2 = pdf('Weibull',t,b,a);
% y2= gampdf(t,K,b);
figure(2), plot(t,y2,xn,y1,'r'), grid
legend('Теорет. Вейбулл-плотность','Гистограмма полученных чисел')
xlabel('x'), ylabel('Histogram');
title(['b = ',num2str(b),' a = ',num2str(a),' K = ',num2str(K)]);
Соседние файлы в папке lab2 вроде