士郎小天使 2022-04-20 11:40 采纳率: 81.1%
浏览 180
已结题

matlab求积分,如何用for循环或者向量求一组结果

如下代码,想把4个y代入这个fun函数中求积分,但一个个求太麻烦了,想一次性求出来,要如何改

y=[2.916;88.631;0.68;11.137];
fun= @(x)exp(y(1)*cos(x))/(2*pi);
y1 = integral(fun,0,360) ;
fun= @(x)exp(y(2)*cos(x))/(2*pi);
y2 = integral(fun,0,360) ;
fun= @(x)exp(y(3)*cos(x))/(2*pi);
y3 = integral(fun,0,360) ;
fun= @(x)exp(y(4)*cos(x))/(2*pi);
y4 = integral(fun,0,360) ;

我试着用for循环,但是这种写法会把4个值加在一起,是哪里有问题吗

y=[2.916;88.631;0.68;11.137];
n=numel(y)
for i=1:1:n
fun= @(x)exp(y(i)*cos(x))/(2*pi);
y0 = integral(fun,0,360) ;
end
  • 写回答

1条回答 默认 最新

  • czy1219 2022-04-20 12:08
    关注

    这样就可以了:

    %方法1
    y=[2.916;88.631;0.68;11.137];
    ret = [];
    for i=1:length(y)
        fun= @(x)exp(y(i)*cos(x))/(2*pi);
        ret=[ret integral(fun,0,360)];
    end
    
    %方法2
    %直接利用矩阵运算
    y=[2.916;88.631;0.68;11.137];
    fun= @(x)exp(y*cos(x))/(2*pi);
    ret=integral(fun,0,360,'ArrayValued',true);
    
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
    1人已打赏

报告相同问题?

问题事件

  • 系统已结题 4月28日
  • 已采纳回答 4月20日
  • 创建了问题 4月20日

悬赏问题

  • ¥15 如何在scanpy上做差异基因和通路富集?
  • ¥20 关于#硬件工程#的问题,请各位专家解答!
  • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
  • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
  • ¥30 截图中的mathematics程序转换成matlab
  • ¥15 动力学代码报错,维度不匹配
  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥15 c++头文件不能识别CDialog