Zaloguj się
Blog na Matlab.pl
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Przepływ ciepła
Forum MATLAB Strona Główna-> MATLAB
Post Wysłany: 25 Grudnia 2016, Nie 12:43 pm Temat postu: Przepływ ciepła Odpowiedz z cytatem
 
AUTOR:
Fasol00005



Dołączył: 12 Paź 2015
Posty: 3


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam, przerabiam jeden z podręczników do nauku matlaba. Jestem na rodziale przepływ ciepła. Zapisany jest w podręczniku gotowy skrypt. Niestety nie działa. Prosze o pomoc.
k=1e-7;
n=10;
L=0.36;
T0 = 25;
TN = -15
z=ones(1,n-1)*(T0+TN)*0.25;
[t,T] = ode45(@rownanie,linspace(0,200000,100),z,[],n);
surf(T)
view([31 22]);
xlabel('warstwa - x')
ylabel('czas')
zlabel('temperatura')

function dy = rownanie(t,T,n)
dy=zeros(n-1,1);
global k L T0 TN
d=n*n*k/(L^2);
dy(1)=d*(T(2)-2*T(1)+T0);
for i=2:n-2
dy(i)=d*(T(i+1)-2*T(i)+T(i-1));
end
dy(n-1)=d*(TN-2*T(n-1)+T(n-2));


 

Post Wysłany: 26 Grudnia 2016, Pon 8:46 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
alhambra
Może pisać książki


Dołączył: 05 Sie 2008
Posty: 690


Ogląda profil użytkownika Wyślij prywatną wiadomość
tak na szybko po obejrzeniu kodu ...

w linijce
[t,T] = ode45(@rownanie,linspace(0,200000,100),z,[],n);

usuń []
bo w sumie dajesz 4 parametry z czego jeden pusty w a funkcji
function dy = rownanie(t,T,n)
masz możliwość implementacji tylko 3 parametrów



_________________
Korepetycje z Matlaba, pomoc ....
 

Post Wysłany: 23 Styczenia 2017, Pon 8:21 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
Fasol00005



Dołączył: 12 Paź 2015
Posty: 3


Ogląda profil użytkownika Wyślij prywatną wiadomość
Niestety to nie działa. Jakiś inny pomysł?


 

Forum MATLAB Strona Główna-> MATLAB
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 35024260
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Pomoc w napisaniu funkcji - operacje na macierzy  (25 Marca 2017, Sob 12:40 pm)
Porównianie wykresów w środowisku matlab i simulink  (24 Marca 2017, Pią 2:14 am)
Komunikat Struct contents reference from a non-struct array  (23 Marca 2017, Czw 9:06 pm)
Problem z zadaniem  (23 Marca 2017, Czw 6:53 pm)
Styczna do odpowiedzi skokowej  (22 Marca 2017, Sro 11:43 am)
Zadanie z Scilab. (Wykonanie skryptu.)  (21 Marca 2017, Wto 5:06 pm)
metody rozwiazywania rownan liniowych  (20 Marca 2017, Pon 10:19 pm)
Metoda Latin Hypercube Sampling  (19 Marca 2017, Nie 7:31 pm)
Kombinacje na zbiorach  (19 Marca 2017, Nie 11:48 am)
Wypełnianie macierzy wektorami  (18 Marca 2017, Sob 8:11 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum