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

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 乌班图ip地址配置及远程SSH
  • ¥15 怎么让点阵屏显示静态爱心,用keiluVision5写出让点阵屏显示静态爱心的代码,越快越好
  • ¥15 PSPICE制作一个加法器
  • ¥15 javaweb项目无法正常跳转
  • ¥15 VMBox虚拟机无法访问
  • ¥15 skd显示找不到头文件
  • ¥15 机器视觉中图片中长度与真实长度的关系
  • ¥15 fastreport table 怎么只让每页的最下面和最顶部有横线
  • ¥15 java 的protected权限 ,问题在注释里
  • ¥15 这个是哪里有问题啊?