第1个回答 推荐于2018-03-26
clear all
clc
syms va rho tau
f=@(va,rho,tau) (2.5-0.1*(va+rho*sin(tau)))./(exp(2.5-0.1*(va+rho*sin(tau)))-1);
n = 5000;%当想要结果更精确时,可以把n设置更大
tau = linspace(0,2*pi,n);
%自编的简单方法
Tn =@(va,rho) pi/n*(sum(f(va,rho,tau(1:end-1)))+sum(f(va,rho,tau(2:end))));
%Tn积分结果变成参数的函数,结果由参数确定。
va = 2;rho = 5;
Tn(va,rho) %这就是数值积分带入参数后的结果,这时积分结果成为参数的函数
%利用matlab自带函数的方法,先确定参数值然后积分。两种方法都可以,后者精度高。
h = @(tau) f(va,rho,tau) %
integral(h,0,2*pi)本回答被网友采纳