MATLAB应用

MATLAB中图形矩阵如何计算坐标，使其能够带进公式进行计算？

这个是我做ifft2之后，提取的虚部，现在要对他进行进一步处理，用公式。但是我不知道怎样把公式写进代码里？在此求助各路大神！

撤销：Ctrl/Command + Z
重做：Ctrl/Command + Y
加粗：Ctrl/Command + B
斜体：Ctrl/Command + I
标题：Ctrl/Command + Shift + H
无序列表：Ctrl/Command + Shift + U
有序列表：Ctrl/Command + Shift + O
检查列表：Ctrl/Command + Shift + C
插入代码：Ctrl/Command + Shift + K
插入链接：Ctrl/Command + Shift + L
插入图片：Ctrl/Command + Shift + G
查找：Ctrl/Command + F
替换：Ctrl/Command + G

代码

I1 = imread(‘E:\work\MATLAB\HRTEM.tif’);
I2 = im2double(I1);
I3= I2.255;
subplot(4,4,1),imshow(I3,[]);
title(‘I3’);
H = fft2(I3);%做fft
K = fftshift(H);%零频分量移到频谱中心，便于后期处理
S = log(abs(K)+1);%返回数组S中每个元素的自然对数ln（x），显示中心低频部分，加对数是为了更好地显示。
subplot(4,4,2),imshow(S,[]);
title(‘S’);
Z=real(K);
x=[900,900,1100,1100];
y=[650,850,850,650];
mask1=roipoly(Z,x,y);%从Z中选取一个感兴趣的多边形区域，根据x，y
mask2=double(mask1);%转换数据类型为精度类型
subplot(4,4,3),imshow(mask2,[]);
title(‘mask2’);
% imagesc(mask2);%将数组mask2的数据显示为一个图像，该图像使用颜色图中的全部颜色
% title(‘mask2’);
S1 =Z. mask2;%选出根据所选区域mask的衍射点
% display(S1);%显示变量S1的值
S2 = ifft2(K);
S3 = fftshift(S2);%将S1的零频分量移到频谱中心并做二维快速傅里叶反变换
subplot(4,4,4),imshow(S3,[]);
title(‘S3’);%添加标题
% disp(S2);%显示变量的值
x = real(S2);%提取实部
subplot(4,4,5),imshow(x,[]);% subplot在平铺位置建立坐标系，imshow显示图像
title(‘x’);
% disp(x);%显示变量x的值
y = imag(S2);%提取虚部
y1=fftshift(y);
subplot(1,1,1),imshow(y1,[]);% subplot在平铺位置建立坐标系，imshow显示图像
title(‘y’);
% disp(y);%显示变量y的值
% P®= -2gpi*u®;%未扭曲晶格的倒易晶格向量
% plot(y);
以上是我的代码，求各路大神指教。在此谢过！

ttps://mermaidjs.github.io/
[3]: https://mermaidjs.github.io/
[4]: http://adrai.github.io/flowchart.js/