mm nn\] = size(im);
w=im;
w = double(w);
%setting finite difference constants
alpha =0.5;
k = 1;
h = 1;
lambda = (alpha\^2)\*(k/(h\^2));
\[m n\] = size(w);
% A matrix form Ax=B linear system
A = zeros(m,m);
% this gen_vec would be rotated and used to populate the matrix A
gen_vec = zeros(1,m);
gen_vec(1,1) = lambda;
gen_vec(1,2) = (1-2\*lambda);
gen_vec(1,3) = lambda;
%filling in values of A matrix
for i=2:m
A(i,:) = gen_vec;
gen_vec = circshift(gen_vec,\[1 1\]);
end
A(1,1) = (1-2\*lambda);
A(1,2) = (lambda);
%making the top-right and bottom-left corners null
A(1:2,n-1:n) = 0;
A(m-1:m,1:2) = 0;
fprintf('size of w: %d\\n',\[size(w)\]);
fprintf('size of w: %d\\n',\[size(A)\]);
w_j_1 = w;
j=1;
figure
for i=1:200 %for each iteration
%multiplication by A on both sides results in diffences in both x and y
## ****🎉3**** ****参考文献****
> 文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。
\[1\]李俊盛,刘宗田.基于异性扩散-中值滤波的超声医学图像去噪方法\[J\].计算机应用与软件, 2009, 26(1):3.DOI:10.3969/j.issn.1000-386X.2009.01.028.
\[2\]张瞳,朱虹,张然,等.复小波域维纳滤波与偏微分扩散相结合的图像去噪方法\[J\].中国图象图形学报A, 2009.
\[3\]王译禾.基于非线性扩散滤波结构信息的图像去噪方法研究\[D\].南京信息工程大学,2016.DOI:10.7666/d.Y3169747.
\[4\]莫绍强.基于各向异性扩散滤波的图像去噪研究\[J\].内蒙古师范大学学报:自然科学汉文版, 2017, 46(1):4.DOI:10.3969/j.issn.1001-8735.2017.01.006.
## [🌈](https://mp.weixin.qq.com/mp/appmsgalbum?__biz=Mzk0MDMzNzYwOA==&action=getalbum&album_id=2591810113208958977#wechat_redirect "🌈")****4 Matlab代码实现****