plot Bethe-Bloch equation

Question

Jeroen 2012-2-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/30411-plot-bethe-bloch-equation

评论： Kunal Kumar 2021-5-25

Dear people,

For a project I have to plot the Bethe-Bloch equation. That is this equation: dedx = (-K*z^2*Z)/(A*beta2)*(0.5*log((2*me*1*10^9*beta2.*gamma2.*Tmax)/I^2)-beta2);

Every time i want to run the program it says the matrix dimension must agree. i do not know what i am doing wrong. I am trying this since the morning and I would like to move on. The only thing that chances is the eng. This has to run from 0 to 0.5 GeV.

clc; clear; z = 1; re = 2.817*10^-15; Z = 29; A = 63.546; Na = 6.022*10^23; mp = 0.93827231; %GeV eng = 0.001:0.001:0.5; %GeV gamma = (eng + mp)/mp; gamma2 = gamma.*gamma; beta = (1.0 - (1.0/(gamma2)')).^0.5; beta2 = beta.*beta; me = 0.51*10^-3; c = 3*10^8; I = 10*Z; K = 4*pi*Na*re^2*me*c^2;

Tmax = (2.0 * me*1e9*beta2'*gamma2)/(1.0 + (2.0*gamma*(me/mp) + (me/mp)^2)); Tmax = Tmax';

dedx = (-K*z^2*Z)/(A*beta2)*(0.5*log((2*me*1*10^9*beta2.*gamma2.*Tmax)/I^2)-beta2);

figure; plot(eng,dedx);

3 个评论
显示 1更早的评论隐藏 1更早的评论

Jeroen 2012-2-27

clc;

clear;

z = 1;

re = 2.817*10^-15;

Z = 29; A = 63.546;

Na = 6.022*10^23;

mp = 0.93827231; %GeV

eng = 0.001:0.001:0.5; %GeV

gamma = (eng + mp)/mp;

gamma2 = gamma.*gamma;

beta = (1.0 - (1.0/(gamma2)')).^0.5;

beta2 = beta.*beta;

me = 0.51*10^-3;

c = 3*10^8;

I = 10*Z;

K = 4*pi*Na*re^2*me*c^2;

Tmax = (2.0 * me*1e9*beta2'*gamma2)/(1.0 + (2.0*gamma*(me/mp) + (me/mp)^2));

Tmax = Tmax';

dedx = (-K*z^2*Z)/(A*beta2)*(0.5*log((2*me*1*10^9*beta2.*gamma2.*Tmax)/I^2)-beta2);

figure;

plot(eng,dedx);

Kunal Kumar 2021-5-25

Jeroen, I am having this same problem. Can you help me with the code. I tried running your code but it shows error. Can you share the code that ran on MATLAB.

请先登录，再进行评论。

请先登录，再回答此问题。