Question

matlab 矩阵维度必须一致

clear all;
T = -15:5:10;
R = 10^(-4);
syms r1 r2;
ymin = @(r1) r1
f = @(r1,r2) r1.*r2.*exp(-R.*pi.*r2.^2).*exp(-R.*pi.*sqrt(T./(r2.^-4+r1.^-4)).*atan(sqrt(T.*r1.^4./(r2.^4+r1.^4))));
inum = integral2(f,0,Inf,ymin,Inf);
q = inum*(2*R*pi)^2
y = arrayfun(inum,T);
plot(T,p,'r.-');
报错 矩阵维度必须一致。请教一下该怎么改呢？

Answer 1

1、“矩阵维度必须一致。”的处理方法

你的问题应该用for循环语句，将T值代入积分表达式内后计算。

2、T = -15～10是有问题的，当T<0时，其积分值是无意义的（NaN）,所以T值应取大于等于0的值。

3、y = arrayfun(inum,T);这句也有问题，应该将 inum 改为 @(r1,r2)inum。

4、按照这三点修改后，运行可以得到图形。