EDA365电子论坛网

标题: 请问程序应该如何写？感谢！ [打印本页]

作者: weiweipcb 时间: 2020-6-19 11:30
标题: 请问程序应该如何写？感谢！
如图所示，要在matlab中求这种非线性较强的积分，做出（x,t）的图像，程序应该如何写？我调用的quadv数值积分的函数，也没有成功，麻烦大神指点一下，非常感谢！

作者: leleeda 时间: 2020-6-19 13:27

可以用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-')
& p) J. T3 p6 s5 Y1 Z$ x+ J

作者: weiweipcb 时间: 2020-6-19 18:10

leleeda 发表于 2020-6-19 13:27
/ u% e& e( p8 A( P) y) P$ N可以用integral函数! ]5 _% M2 G( o0 a) _2 F9 }

syms t real

大神太牛了，非常感谢！

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)