clear; close all; clc; L = 1; x = linspace(0/L,L/L,100); t = [0.01 0.02 0.05 0.10 0.20]; m = 0; sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t); plot(sol(1,:)); hold on; plot(sol(2,:)); plot(sol(3,:)); plot(sol(4,:)); plot(sol(5,:)); legend(st…

clear; close all; clc; L = 1; x = linspace(0/L,L/L,100); t = [0.01 0.02 0.05 0.10 0.20]; m = 0; sol = pdepe(m,@heatpde,@heatic,@heatbc,x,t); plot(sol(1,:)); hold on; plot(sol(2,:)); plot(sol(3,:)); plot(sol(4,:)); plot(sol(5,:)); legend(st…

以下を参考に。 www.mathworks.com 方程式の定義 以下の熱伝導方程式を解く 上記を以下と比較し、パラメータm、f,、c、sを決定する。 コード（mはメインの関数で用いるのでここでは使わない） function [c,f,s] = heatpde(x,t,u,dudx) c = 1; f = dudx; s = …

燃焼ガスの熱伝達係数はNakkaの値（1000）を用いた。 www.nakka-rocketry.net clear; close all; clc; L = 2e-3; %[m]. 2mm Tin = 1000; %K Tinf = 300; %K alpha_in = 1000;% 燃焼ガスの熱伝達係数 by Nakka % https://www.nakka-rocketry.net/therm.html#:…