F分布噪声图像

clc,clear,close all
warning off
feature jit off
im = imread('coloredChips.png');
Z1 = imnoise_F(size(im,1),size(im,2),3,3);
Z1 = uint8(Z1);   % 类型转换
figure('color',[1,1,1]),
im(:,:,1) = im(:,:,1) + Z1;  % R
im(:,:,2) = im(:,:,2) + Z1;  % G
im(:,:,3) = im(:,:,3) + Z1;  % B
subplot(121); imshow(im);title('加F分布噪声图像')
subplot(122); imhist(Z1); title('加F分布噪声图像直方图')

function R = imnoise_F(M, N, a, b)
% input:
%       F分布，噪声的类型；
%       M,N：输出噪声图像矩阵的大小
%       a,b：各种噪声的分布参数
% output:
%       R: 输出的噪声图像矩阵，数据类型为double型
% 设定默认值
if nargin < 4
   a = 1;b=1;
end
   % 产生F分布噪声
   R1=zeros(M,N);
   R2=zeros(M,N);
   for i=1:a
        R1 = R1 + 5*randn(M, N).^2;
        R2 = R2 + 5*randn(M, N).^2;
        R = R1./R2;
   end
end

F分布的性质
　　1、它是一种非对称分布；
　　2、它有两个自由度，即n1 -1和n2-1，相应的分布记为F（ n1 –1， n2-1）， n1 –1通常称为分子自由度， n2-1通常称为分母自由度；
　　3、F分布是一个以自由度

   

和

   

为参数的分布族，不同的自由度决定了F 分布的形状。
　　4、F分布的 倒数性质：

   

。