Сеточные методы в задачах аэродинамики

% UIWAIT makes Polet wait for user response (see UIRESUME)

% uiwait(handles.figure1);

 

 

% --- Outputs from this function are returned to the command line.

function varargout = Polet_OutputFcn(hObject, eventdata, handles)

% varargout  cell array for returning output args (see VARARGOUT);

% hObject    handle to figure

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

 

% Get default command line output from handles structure

varargout{1} = handles.output;

 

 

% --- Executes on button press in calculate.

function calculate_Callback(hObject, eventdata, handles)

% hObject    handle to calculate (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

clc

% Начальные условия

x0 =0;

     Продолжение ПРИЛОЖЕНИЯ А

     Реализация  метода конечно-разностных схем 

y0 =1000;

dt = 0.001;

a = 10;

m= 20;

%v0 = 6;

g=9.8;

Vpr = 7.5;

fi=pi/3;

% считывание данных  с полей интерфейса

Q1=get(handles.edit1,'String');

Q2=get(handles.edit2,'String');

Q3=get(handles.edit3,'String');

Wv=str2double(Q1);

Vv=str2double(Q2);

v0=str2double(Q3);

x(1) = x0;

y(1) = y0;

v(1) = v0;

 

u(1) = v(1)*cos(a);

w(1) =v(1)*sin(a);

i = 2;

t=false;

while t==false

x(i) = x(i-1)-dt*(u(i-1)+Vv);

y(i) = y(i-1)-dt*(w(i-1)+Wv);

u(i) = u(i-1)-dt*(   (  (-g./Vpr.^2)  * ( v(i-1) * (u(i-1)+Vv))   ) - ( (g./Vpr.^2) * ( ( u(i-1)*cos(m)+w(i-1)*sin(m) )./( u(i-1)*sin(m)-w(i-1)*cos(m) )  ) *( v(i-1)*(w(i-1)+Wv ))) )  ;

w(i) = w(i-1)-dt*( (-g-((g./Vpr.^2) * (v(i-1)* (w(i-1)+Wv))  ))+ ( (g./Vpr^2) * ( ( u(i-1)*cos(m)+w(i-1)*sin(m) )./( u(i-1)*sin(m)-w(i-1)*cos(m) )  ) * (v(i-1)*(u(i-1)+Vv))))   ;

v(i) = ((u(i-1)+Vv).^2 + (w(i-1)+Wv).^2).^0.5;

 

if y(i)<=0

    t=true;

end

i=i+1;

end

 

 

DL = x(i-1); %Длина пути

r(1)=0;

l(1)=1000

DL/1000

for k=2:1000

    r(k)=(k-1)*DL/1000;

    l(k)=l(1)-k+1;

end; 

     Продолжение ПРИЛОЖЕНИЯ А

     Реализация  метода конечно-разностных схем 

plot(r/1000,l,'-');

rass=DL/1000;

set(handles.rass,'String',rass);

 

function edit1_Callback(hObject, eventdata, handles)

% hObject    handle to edit1 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

 

% Hints: get(hObject,'String') returns contents of edit1 as text

%        str2double(get(hObject,'String')) returns contents of edit1 as a double

 

 

% --- Executes during object creation, after setting all properties.

function edit1_CreateFcn(hObject, eventdata, handles)

% hObject    handle to edit1 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

 

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

 

 

 

function edit2_Callback(hObject, eventdata, handles)

% hObject    handle to edit2 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

 

% Hints: get(hObject,'String') returns contents of edit2 as text

%        str2double(get(hObject,'String')) returns contents of edit2 as a double

 

 

% --- Executes during object creation, after setting all properties.

function edit2_CreateFcn(hObject, eventdata, handles)

     Продолжение ПРИЛОЖЕНИЯ А

     Реализация  метода конечно-разностных схем 

% hObject    handle to edit2 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

 

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

 

 

 

function edit3_Callback(hObject, eventdata, handles)

% hObject    handle to edit3 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

 

% Hints: get(hObject,'String') returns contents of edit3 as text

%        str2double(get(hObject,'String')) returns contents of edit3 as a double

 

 

% --- Executes during object creation, after setting all properties.

function edit3_CreateFcn(hObject, eventdata, handles)

% hObject    handle to edit3 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    empty - handles not created until after all CreateFcns called

 

% Hint: edit controls usually have a white background on Windows.

%       See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))

    set(hObject,'BackgroundColor','white');

end

 

 

% --- Executes during object creation, after setting all properties.

function rass_CreateFcn(hObject, eventdata, handles)

% hObject    handle to rass (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

     Продолжение ПРИЛОЖЕНИЯ А

     Реализация  метода конечно-разностных схем 

% handles    empty - handles not created until after all CreateFcns called

% rass=DL;

% set(handles.rass,'String',rass);

 

 

% --- Executes on button press in pushbutton2.

function pushbutton2_Callback(hObject, eventdata, handles)

% hObject    handle to pushbutton2 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

close; 
 

Рисунок 8 – Интерфейс пользователя