I am using this program and I get a message saying that my licence is not valid
1 view (last 30 days)
Show older comments
Alain Gibaud
on 29 Aug 2017
Answered: Walter Roberson
on 29 Aug 2017
clf
Npt=128;
d=0.2;
%calcul de la differencielle de l'integrale de delta
%v=repmat(0,[Npt,Npt,Npt]);
v=zeros(Npt,Npt,Npt);
v(5:Npt-5,5:Npt-5,5:Npt-5)=0.9;
v(30:70,30:70,30:70)=0.05;
v(14:17,14:17,14:17)=0.5;
%v(15:16,15:16,15:16)=1;
v=v*1e-6;
pix=1.4e-6;
T=sum(v,3)*pix;
A=diff(T,2,1)/pix^2;
[N1A, N2A]=size(A);
A(N1A:Npt,:)=0;
B=diff(T,2,2)/pix^2;
[N1B, N2B]=size(B);
B(:,N2B:Npt)=0;
%calcul de mu
mu=v*1e-9;
u=sum(mu,3);
%calcul du second terme
Z=u-d*(A+B);
figure(100)
I0=1e10;
I=I0./(exp(Z));
IA=mean(mean(I));
surf(I-IA)
view(2)
axis tight
colormap(gray)
colorbar
figure(200)
surf(u)
view(2)
axis tight
colormap(gray)
colorbar
0 Comments
Accepted Answer
More Answers (0)
See Also
Categories
Find more on Subplots in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!