请问如何用matlab解这个积分方程组??
发布网友
发布时间：2022-04-26 04:38
共1个回答
热心网友
时间：2022-06-20 20:22
我在想一个问题，你对N、δ、η的数量级有没有一个认识我编程之后一直无法求出，按理应该可以求出的
clc;clear[x,~,jud]=fsolve(@myfun,[1000,2,1e-3]);
while(jud~=1)
x0=rand(3,1);
[x,~,jud]=fsolve(@myfun,x0);
end
disp(x)
disp(x0)
function ff=myfun(x)
ff=[32300-pi/6*2.66*x(1)*integral(@(d)funt(d,x(1),x(2)),0,0.43e-3);
26100-pi/6*2.66*x(1)*integral(@(d)funt(d,x(1),x(2)),0.43e-3,0.65e-3);
33100-pi/6*2.66*x(1)*integral(@(d)funt(d,x(1),x(2)),0.65e-3,1.1e-3);
];
end
function f=funt(x,a,b)
f=x.^2.*exp(-((log(x)-a).^2/2/b^2))/sqrt(2);
end
追问老师提示选择一个算法，什么遗传算法好像是，具体值算不出来，只能求出最接近的解，也就是用特殊值代进去看看等式最接近哪个，用matlab代入特殊值，特殊值里面好像有一个7左右的还有一个8左右的值，老师就提示这些！
老师提示选择一个算法，什么遗传算法好像是，具体值算不出来，只能求出最接近的解，也就是用特殊值代进去看看等式最接近哪个，用matlab代入特殊值，特殊值里面好像有一个7左右的还有一个8左右的值，老师就提示这些！
