请问程序应该如何写？感谢！

weiweipcb

如图所示，要在matlab中求这种非线性较强的积分，做出（x,t）的图像，程序应该如何写？我调用的quadv数值积分的函数，也没有成功，麻烦大神指点一下，非常感谢！
5 y) m. Z% `9 C! V$ j+ ]; ?* u6 ~
5 q4 O; m! x" b$ u

leleeda

可以用integral函数

• syms t real
• t_2 = 0.0976;
• L = 0.3;
• Ee = 2.068e+11;
• x = 0:1e-3:0.2135;
• t_1 = x - 0.0226;
• result = zeros(size(x));
• result = arrayfun(@(ii) integral(@(t) (t_2-t).*cos(t).*( sin(t_1(ii)) ).^2./(2.*Ee.*L.*( sin(t) + (t_2 - t).*cos(t) )), t_1(ii), t_2 ), [1:1:numel(x)]);
• plot( x, result, 'ko-')
  

weiweipcb

leleeda 发表于 2020-6-19 13:27
6 J" ?9 s! Y7 N1 z' E9 [可以用integral函数
! S; N, B# B& `* U1 l4 S  r

syms t real

4 a2 l& o: A* E' J+ |- ~* i6 N大神太牛了，非常感谢！