Graph Plotting the differential equation

1 次查看（过去 30 天）
Dhananjay Singh
Dhananjay Singh 2020-11-23
评论： Stephan 2020-11-23
How I plot qSol(r) wrt r and TSol(r) wrt r ?
%heat conduction with an electrical heat source in copper wire
clc; clear;
close all;
%parameters
R = 0.2; % radius in m
L = 10; % length of wire in m
ke = 60000000; % electrical conductivity of copper in S/m
k = 401; % thermal conductivity of copper in W/m.K
I = 10; % electrical current in Amp
Se = (I*I)/ke; % rate of heat production per unit volume
T0 = 20 % Temp in Celcius
%Energy Balance in cylindrical shell of thickness del(r) and length L
% Doing the energy balance and taking del(r) common and limiting it to 0, a
% linear differential equation is formed.
% SOLVING DIFFERENTIAL EQUATION
syms q(r) %q is the heat flux
ode = diff(r*q,r) == Se*r
cond = q(0)==0 %B.C at r = 0,q~=inf gives a expression which states that at r = 0 , q = 0
qSol(r)=dsolve(ode,cond)
%Temperature Profile
syms T(r)
ode = diff(T,r) ==(-Se*r/(2*k))
cond = T(R)==T0;
TSol(r)=dsolve(ode,cond)
%Plot Flux profile

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

采纳的回答

Stephan
Stephan 2020-11-23
编辑：Stephan 2020-11-23
yyaxis left
fplot(qSol,[0, 100])
yyaxis right
fplot(TSol,[0, 100])

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by