【灵敏性】基于巴顿模型计算输入空间频率的对比敏感度值研究（matlab代码实现）

 👨‍🎓个人主页：研学社的博客  

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

使用Barten模型计算输入空间频率的对比敏感度值。输入参数包括空间频率（周期/度）、显示亮度 （cd/m2）、周围亮度 （cd/m2） 和视场（度）。

📚2 运行结果

部分代码：

function [ CSF_out ] = CSF( fq,L,Ls,X )
% Calculate contrast sensitivity (Barten model) for a series of spatial frequencies
%   fq: spatial frequency in cycle/degree, should be a scalar or a vector
%   L: display luminance in cd/m2
%   Ls: surround luminance in cd/m2
%   X: Field of view in degree
%   CSF: retured contrast sensitivity
%   CSF equation from: P. Barten, 揊ormula for the contrast sensitivity of the human eye,� Proc. SPIE 5294, 231�238 (2003).
%   By Zong Qin (qinzong.wnlo@gmail.com) from National Chiao Tung University, Taiwan

Itm1=5400*exp(-0.0016*fq.^2*(1+100/L)^0.08);
Itm2=1+144/(X^2)+0.64*fq.^2;
Itm3=63./(L^0.86)+1./(1-exp(-0.02*fq.^2));
CSF_uncorrected=Itm1./sqrt(Itm2.*Itm3);

Itm4=(Ls/L)*(1+144/(X^2))^0.25;
Itm5=(1+144/(X^2))^0.25;
Itm6=((log(Itm4))^2-(log(Itm5))^2)/(2*log(32)*log(32));
Surround_factor=exp(-Itm6);
CSF_out=CSF_uncorrected.*Surround_factor;

end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]夏鹏,杨特,王乐,杨智春.平稳随机载荷的灵敏度分析识别方法[J].振动与冲击,2022,41(23):300-306.DOI:10.13465/j.cnki.jvs.2022.23.035.

[2]P. Barten, “Formula for the contrast sensitivity of the human eye,” Proc. SPIE 5294, 231–238 (2003).

🌈4 Matlab代码实现