lzth.net
当前位置:首页 >> mAtlAB 积分函数作图问题 >>

mAtlAB 积分函数作图问题

在这个积分式子中,如果除y之外都是常量,那就存在一个问题:y是被积分的变量,积分之后得到的结果与y无关,只能是一个常数而已,怎么画图? 参考代码(各常数的值随便取的): M = 1;G = 2;H = 3;a = 4;b = 5;F = 1/M * ( G*quadl(@(y)1./(H/(M...

f = @(x,t)exp(-t).*(1-exp(-(x/2).^12.3.*t.^2.2));F = @(x)quadgk(@(t)f(x,t), 0, inf);fplot(F, [-6 6])

说明 这个问题和另外两个问题(编号2051722037141864067、1638082848257894860)基本上是重复的,我已经在那两个问题做了回答,主要原因是匿名函数f0在x比较大、s比较小的情况下会出现非数NaN,导致计算失败。更具体的分析与解决方法这里不再重...

这个问题的困难不在于积分限里有变量,而在于积分限里有无穷大。 如果用数值积分,quad系列函数都不能处理积分区域为无界的情况,但2012a引入了integral系列新的数值积分函数可以处理。 不知道你的函数F(x,y)表达式,这里举个例子来说明。 假设F...

函数Q不是数组,需要赋值 不采用sym定义,Conversion to double from “sym” is not possible. T=300:1:600; U=293./(273+T)*2100/15; Q=-(x-U).^2./(U.^2)/0.0032; p=int(1/(sqrt(2*pi)*8.4)*exp(Q),x,-inf,sqrt(2/3)*110); plot(T,Q);

下面是求sin(t)/t的一个函数:变上限积分,另外提醒你:这个函数没有对应的初等函数形式的原函数,也就是说不定积分没有解析解!! clear,clc t=linspace(0,10,100); f=@(t) sin(t)./(t+eps); F=t; n=length(t); for i=1:n F(i)=quad(f,0,t(i)); ...

用matlab绘制这类复杂的积分函数图形,应采用循环语句求出,当T为某值时的C值。 实现方法: x=20:1:70; for i=1:length(x) T=x(i); C(i)=fun_ex(T); end T=x; plot(x,C,'*-'),grid on; title('T--C(T)曲线图'); xlabel('T');ylabel('C(T)'); 运...

10sin(0.03k) Int exp(-0.08(r-1)^2) dr。 程序如下: %用 MATLAB 算 Int exp(-0.08(r-1)^2 dr 这个积分 F = @(r)exp(-0.08*(r-1).^2); R = 25; Q = quad(F,0,R); %然后再在结果上乘以 10sin(0.03k) 的部分 k = 1:628; F1 = 10*sin(0.03*k) * ...

clear all;clc; fx=@(x)(x.^(-1)-1).^(1/3); x=0:0.01:0.8; for ii=1:length(x) x0=x(ii); Fx(ii)=quadl(fx,x0,1); end plot(x,Fx);

syms x Ef Et y; Ei=0.5426; Ef=0.842; Et=linspace(1.1,1.4,100); %%%%好像是9.11e-31 n=((4*pi*(2*1.18*9.11*10^31)^(3/2)/(4.14*10^(-15))).*(sqrt(x-1.1)/(1+2.71.^((x-Ef)./0.026)))); for i=1:length(Et) i y(i)=double(int(n,x,1.1,Et(i)...

网站首页 | 网站地图
All rights reserved Powered by www.lzth.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com