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

查看数: 581 | 评论数: 2 | 收藏 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2020-6-19 11:30

如图所示，要在matlab中求这种非线性较强的积分，做出（x,t）的图像，程序应该如何写？我调用的quadv数值积分的函数，也没有成功，麻烦大神指点一下，非常感谢！% R5 K6 R# T% G. q & U: l, h- |7 X ...

回复

weiweipcb 发表于 2020-6-19 18:10

leleeda 发表于 2020-6-19 13:27
; k# X5 T7 t. p' f# l可以用integral函数# o$ m2 ?: U b0 S7 L# C4 h

syms t real

大神太牛了，非常感谢！

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-')- M% c" ]2 J% s) }0 c