Question

MATLAB求数组积分

K是一个数组，怎样求小k的积分，并且在MATLAB中画出图来

Answer 1

参考代码：

k=arrayfun(@(K)sin(150*quad(@(t)K*cos(t).^2,-pi/2,pi/2)),10:5:80)
k =
   1.0e-05 *
  Columns 1 through 8
    0.1505   -0.0822   -0.0475    0.0593   -0.2131    0.1482   -0.0692    0.0778
  Columns 9 through 15
    0.0045   -0.0050    0.0054   -0.0059    0.0063   -0.0068    0.0072

Answer 2

syms sita
k=[];
for K=10:5:80;
h=quad('cos(sita).^2',-pi/2,pi/2);
k=[k sin(100*K*h)];
end
k
k =

1.0e-004 *

Columns 1 through 12

-0.1014 -0.1520 -0.2027 -0.2534 -0.3041 -0.3547 -0.4054 -0.4561 -0.5068 -0.5575 -0.6081 -0.6588

Columns 13 through 15

-0.7095 -0.7602 -0.8109

追问

如果大K是一个带有sita的数组呢？