Классическая теория рассеяния

@progger228 · Регистрация: 15.10.2015

Author24 — интернет-сервис помощи студентам

Всем приветики. Имеется программа для моделирования рассеяния классической частицы в кулоновском потенциале. Вроде всё верно, но matlab выдаёт ошибку связанную с ode45

Error using feval
Undefined function 'CenterPole' for input arguments of type 'double'.

Error in odearguments (line 87)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ode45 (line 115)
odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options,
varargin);

Error in rasseivanie (line 19)
[T Y]=ode45('CenterPole',[0:dt:Tmax],[x0 vx0 y(i) vy0]); %#ok<NBRAK>

Сам текст программы:

Matlab M

clear all
global Q q M
Q = 79; %заряд рассеивающего центра
q = 2; %заряд рассеиваемой частицы
M = 4; %масса частицы
 
Tmax = 2*10^7; %длительность временного интервала, на котором
%ищется численное решение системы ОДУ
dt = Tmax/2000; %шаг интегрирования
x0=-10^6;%
vx0=0.1;
vy0=0;
Ymin=0;
Ymax=10^5;
N=100;
dy=(Ymax-Ymin)/N;
for i=1:N+1
    y(i)=(i-1)*dy;
    [T Y]=ode45('CenterPole',[0:dt:Tmax],[x0 vx0 y(i) vy0]); %#ok<NBRAK>
    M1=length(T);
    xi(i)=acos(Y(M1,2)./(Y(M1,2).^2+Y(M1,4).^2).^0.5);
end;
figure(1); plot(y,xi);
x1=min(xi);
x2=max(xi);
Nxi=1000;
dxi=(x2-x1)/Nxi;
i=1:Nxi+1;
Yi=interp1(xi,y,Xi,'splain');
dro=abs(diff(Yi)./dxi);
j=1:Nxi;
dsigma(j)=2*piYi(j).*dro(j);
xi1(j)=Xi(j);
plot(xi1,dsigma);

В чем может быть ошибка?

@AlMih · 15.12.2019, 22:14

А где функция CenterPole?

@progger228 · 15.12.2019, 23:21 **[ТС]**

Matlab M

function dy=CenterPole(t, z)
global Q q M
A = 1.53*q*Q/M;
dy=zeros(4,1);
dy(1)=z(2);
dy(2)=A*z(1)/(z(1)^2+z(3)^2)^(3/2);
dy(3)=z(4);
dy(4)=A*z(3)/(z(1)^2+z(3)^2)^(3/2);

Добавлено через 13 минут
Нашёл ошибку. Этой ошибкой оказался я.

@AlMih 416 / 200 / 69 Регистрация: 20.01.2019 Сообщений: 713
	15.12.2019, 22:14	2
	А где функция CenterPole? 0

	15.12.2019, 23:21