关于磁共振图像重建margosian的实现所用到的matlab函数（小白向）

Posted 2021-02-19 w-yan

1、phantom函数
简介：生成头部幻影图像
举例：p=(phantom(N)); imshow(p);%会出现一个头部幻影图像，其中N表示图像矩阵的行列数。
其他：p=phantom(def,n)，这个是标准格式，上面的例子是助教代码中的例子。

2、zeros函数
简介：生成全是零的矩阵
举例：phase=zeros(N,N);%矩阵的行列数均为N，zeros(m,n)生成的是m*n的矩阵

3、exp(1i*2*pi*(iter_x/N - 1/2).^2)的理解
介绍：1i表示复数，虚部为1；.^表示乘幂计算

4、矩阵的点乘（.*）
介绍：矩阵对应位置的数值相乘，比如二阶矩阵[1,2;2,3]，点乘自身得到的是[1,4;4,9]

5、fftshift函数
简介：将零频分量移到频谱中心
举例：fftshift(X,dim);%dim=1时，是列的上下变换；dim=2时，是行的左右变换。
补充：fftshift(fftshift(X,1),2)表示将矩阵中心变换

6、fft函数
简介：快速傅里叶变换
举例：fft(X,n,dim);%表示将X进行快速傅里叶变换，n表示采样点数量，dim表示维度，1、2表示行与列；
补充：助教的代码是fft(X,[],dim)表示自动选取采样点数量

7、randn函数
简介：正态分布随机数
举例：randn(N);%返回N*N矩阵的正态分布随机数

8、complex函数
简介：创建复杂数组
举例：c=complex(a,b);%返回元素是复数的矩阵，c=a+bi

9、round函数
简介：四舍五入到最接近的小数或者整数
举例：round(x);%四舍五入到最接近的整数
补充：round(x,N)表示四舍五入到小数位数为N的小数

10、ifft函数
简介：傅里叶逆变换
用法：和fft极其相似

11、flipud函数
简介：矩阵上下翻转
补充：fliplr函数表示矩阵左右翻转

12、hamming函数
简介：生成hamming fliter
用法：hamming（n）

14、real函数
简介：求复数的实部
用法：real(x)

15、repmat函数
简介：复制矩阵
用法：repmat（X,m,n)，将矩阵X复制m行n列

16、angle函数

简介：求相位

用法：你猜！