参考代码：
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
温馨提示：答案为网友推荐，仅供参考
第1个回答 2014-04-01
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=[];
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的数组呢？
本回答被提问者和网友采纳
相似回答