matlab 求变上限积分方程?

作者&投稿:迪疮 (若有异议请与网页底部的电邮联系)
matlab求变上限积分方程~

syms x t

solve(int(t,0,x)-x)
ans =
0
2

参考代码:
>> syms x>> g=sin(x)/x;>> h=int(g,0.1,x);>> h=int(g,0.1,x)h =sinint(x) - sinint(1/10) >> ezplot(g)>> figure>> ezplot(h)>> diff(subs(h,x,x^2),x,2)ans =4*cos(x^2) - (2*sin(x^2))/x^2

syms x y y1 t f1 f2 f3 f4 f5 f22 f32f1= exp(-(t-50).^2./200) ;f2=int(f1,-Inf,y);f3=int(f1,-Inf,0.6*y);f22=int(f2,0,y1);f32=int(f3,0,0.6*y1);i=1;for x1 = 0.0: 0.2: 1.0,f4 = (1-x1)*f22 + x1*f32 -75*sqrt(2*pi);f5 = f4;y11(i) = solve(f5,'y1');x11(i) = x1;i=i+1;endyout = [x11; single(y11)];fprintf('X is %4.2f and Y is %10.4f
', yout)plot(x11,single(y11))计算结果:X is 0.00 and Y is 55.7412 X is 0.20 and Y is 58.2157 X is 0.40 and Y is 61.9204 X is 0.60 and Y is 68.5814 X is 0.80 and Y is 86.7197 X is 1.00 and Y is 141.6185 用2010a以下版本计算,可能有些误差。可以修改for x1 = 0.0: 0.2: 1.0来取得更多的计算结果,但电脑计算时间较长。这是for x1 = 0.0: 0.2: 1.0时的图形:



这是for x1 = 0.0: 0.01: 1.0时的图形:






英德市17758172368: matlab求变上限积分方程 -
惠兔丝白: f=@(x)quad(@(t)(1./sin(t).^1.5),x,pi/2)f(x)即为所要的变下限积分; 如f(0.1)即是下限为0.1 上限为pi/2的积分值fplot(f,[0.01,pi/2])给出了f的部分图像.

英德市17758172368: MATLAB求解变上限积分方程 -
惠兔丝白: 例:积分下限为0,上限是T,T的取值范围是[5,125],步长为5; syms y T x y=int('f(x)',0,T);T=5:5:125;subs(y)

英德市17758172368: 如何用matlab求解变上限积分 -
惠兔丝白: 直接使用CFTOOL工具箱 命令行输入cftool即可,然后选择拟合类型 x=[6.69:0.02:7.53]; y=[4.2,3.7,3.3,2.95,2.63,2.33,2.11,1.87,1.65,1.47,1.32,1.17,1.04,0.925,0.82,0.735,0.653,0.582,0.52,0.462,0.412,0.366,0.325,0.289,0.258,0.23,0.205,0.182,0....

英德市17758172368: matlab中如何求解如下变上限积分方程? -
惠兔丝白: x=0.6时,y=68.581359956511627574834151414947 close all clear,clc x=0:0.05:1; y=zeros(1,length(x)); syms y Y ex1=int(int('exp(-(t-50)^2/200)',-inf,y),0,Y); ex2=int(int('exp(-(t-50)^2/200)',-inf,0.6*y),0,0.6*Y); ex1=ex1/(75*sqrt(2)*pi^(1/2)); ex2=ex2/(...

英德市17758172368: matlab求变上限积分方程等式的解该怎么求 -
惠兔丝白: 这个可以这样做,用符号积分 >> syms x y; >> y=int('x^2',1,x)%求出了该积分的表达式 solve('x^3/3 - 1/3','x')求解这个表达式就行了~~ ans = 1- 1/2 - (3^(1/2)*i)/2- 1/2 + (3^(1/2)*i)/2

英德市17758172368: MATLAB求解变上限积分方程
惠兔丝白: clc; clear all; close all;syms r neq = 1/(sqrt(2*pi)*50)*exp(-(r-500)^2/5000);eqin = int(eq, r, 0, n);eqin = eqin - 5/8;eqin = vpa(eqin, 5)solve(eqin, n) 结果eqin = 0.5*erf(0.014142*n - 7.0711) - 0.125 Warning: Explicit solution could not be found. > In solve at 81 In gm at 7 ans = [ empty sym ] >>

英德市17758172368: 怎样求用MATLAB求变上线积分
惠兔丝白: 举个例子:比如说sin(t)/t在0到x上的积分可以通过以下代码实现: syms y; y=int('sin(t)/t','t',0,'x'); 求得的y就是积分上限为x的积分!

英德市17758172368: matlab中求变限积分的方法,并画图!!! -
惠兔丝白: 假如求函数名为mat的函数的在以x=0为下限的不定积分for i=1:20x(i)=0.1*i;y(i)=quad('mat',0,x(i));end,plot(x,y) 祝你学习愉快啊.

英德市17758172368: 如何用matlab求变上限的定积分 -
惠兔丝白: 这个可以这样做,用符号积分>> syms x y;>> y=int('x^2',1,x)%求出了该积分的表达式 solve('x^3/3 - 1/3','x')求解这个表达式就行了~~ ans =

英德市17758172368: matlab程序编写变上限积分的函数我遇到一个公式,需要用matlab粗计算,但是我不会使用matlab语言编写变上 -
惠兔丝白: ^举个例子,f(x)为t^2,t从0积分到x f=@(x)quadl(@(t)t.^2,0,x); 然后你就可以代入任何一个x求f(x)了,比如 f(2) 你甚至可以画出f(x),比如 fplot(f,[-2,2])

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网