matlab 求已知概率密度函数的随机数生成

N=10000; %需要随机数的个数
a=zeros(N,1); %存放随机数的数列
n=0;
f1=@(t) 1./(1.2*pi*(1+5*(t-7.3).^2));
f2=@(t) 1./(1.05*pi*(1+6*(t-8.2).^2));
tt=linspace(0,24,1000);
ff=f1(tt).*(tt<=7.5)+f2(tt).*(tt>7.5);%根据公式计算概率密度
s=trapz(tt,ff);  %计算整个区间概率密度的积分
ff=ff/s;         %归一化概率密度
 
 
while n