用MATLAB计算函数的积分

很基础但是很实用，自己写的时候知道一些指令，但是搭配在一起就用的不太好

1 写一个函数句柄结构

A=@(x) 36*x+11;

N_quadrat=@(x) (c.^2)./(omega0.^4).*x.^4.*(1+((2.*D./omega0).*x).^2)./((1-((1./omega0).*x).^2).^2+((2.*D./omega0).*x).^2); 

%N_quadrat 是x的函数

 这种都是常见的语法格式，只不过函数的格式简繁不一。

2 如果有函数之间的运算记得构建新的函数

 s= @(x) 0.0000122*x.^(-2); % s是x的函数

 fun = @(x) N_quadrat(x).*s(x); %fun也是x的函数（或表达式），但是是两个函数的成绩，为了避免错误尽量都使用点运算

3 积分，使用指令 integral

N_z_eff = integral(fun,0.1,1000) %语法格式integral （函数名称，积分下限，积分上限）

 

4 搭配fplot指令来把结果形象化

例子

我们来计算y=x^2 在0到10上的积分

附代码

clear all;
close all;
clc;

y=@(x) x.^2;
Res=integral(y,0,20)
fplot(y,[0,20])