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

weiweipcb

如图所示，要在matlab中求这种非线性较强的积分，做出（x,t）的图像，程序应该如何写？我调用的quadv数值积分的函数，也没有成功，麻烦大神指点一下，非常感谢！
4 \- b% g4 _% l( y/ l9 E5 t

leleeda

可以用integral函数
9 C, |$ R6 Q" M

• 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-')
  & p- m; X. S: n. t% s

weiweipcb

leleeda 发表于 2020-6-19 13:27
) [; Q1 s, [, z' U( `可以用integral函数

syms t real

大神太牛了，非常感谢！
8 B/ O5 Q- d. c' Y* p  l4 t7 b