资源简介
求解流体力学中有限体积(FVM)法的matlab代码
代码片段和文件信息
function datagen
clc;
clear all;
format long;
h=0.01; %%每隔h秒传感器测一次温度,滤波器工作一次
nc=1000; %%温度传感器测量次数
time=linspace(0.0110.0nc);
%%%%%%%%%%%%%%%%%初边界条件设置%%%%%%%%%%%%%%%%%%%%%
cq=1.0;%%热流
rcp=1.0; %%密度和比热的乘积
SL=1.0; %%传热长度
%%%%%%%%%%%%%%%%%初边界条件设置完毕%%%%%%%%%%%%%%%%%
nx=50; %%控制体积的个数--nx
nt=10; %%每个滤波步长内积分步数--nt
T0=0.0*ones(1nx); %%每个控制体积的初始温度
%%%%%%%%%%%%%%%%%%热传导系数确定%%%%%%%%%%%%%%%%%%%
theta=zeros(13);
theta(1)=1.0;
theta(2)=0.0e-003;
theta(3)=0.0e-005;
%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%%%
Tout=zeros(1nx); %%t+dt时刻的各控制体积的温度分布
TC=zeros(1nc); %%各时刻的温度测量值
CK=zeros(1nc); %%各刻传系数的记录值
realTC=zeros(1nc);
M_ran=random(‘norm‘0.01.01nc);
%%%%%%%%%%%%%%%%%%%%开始温度测量%%%%%%%%%%%%%%%%%%%%
Tin=T0;
for i=1:nc
[ToutTCK] = fvm(TinthetanxSLnthrcpcq); %%[ToutTCK] = feval(‘fvm‘TinthetanxSLnthrcpcq);
display(i);
Tin=Tout;
TC(i)=TCK+0.001*M_ran(i);
realTC(i)=TCK;
CK(i)=theta(1)+theta(2)*TCK+theta(3)*TCK*TCK;
end
%%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%解析解%%%%%%%%%%%%%%%%%%%%%%%
Tjx=zeros(1nc);
for i=1:nc
vn=0.0;
for j=1:1000
vn=vn+2.0*exp(-j*j*pi*pi*time(i))*cos(j*pi)/(j*j*pi*pi);
end
Tjx(i)=T0(nx)+time(i)-1.0/6.0-vn;
end
%%%%%%%%%%%%%%%%%%%%%%%%end%%%%%%%%%%%%%%%%%%%%%%%%%
figure;
plot(timeTC‘b‘timeTjx‘r‘);
figure;
plot(timeTC-Tjx‘r‘);
%%%%%%%%%%%%%%%%%%%%测量值存储%%%%%%%%%%%%%%%%%%%%%%
save TC.dat TC -ascii;
save CK.dat CK -ascii;
save realTC.dat realTC -ascii;
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 16002 2014-04-03 09:03 有限体积法data_gen\data_gen\CK.dat
文件 1418 2009-09-07 09:51 有限体积法data_gen\data_gen\datagen.asv
文件 1745 2009-09-21 22:13 有限体积法data_gen\data_gen\datagen.m
文件 3336 2009-09-08 11:29 有限体积法data_gen\data_gen\fvm.m
文件 16002 2014-04-03 09:03 有限体积法data_gen\data_gen\realTC.dat
文件 16002 2014-04-03 09:03 有限体积法data_gen\data_gen\TC.dat
目录 0 2014-04-04 10:50 有限体积法data_gen\data_gen
目录 0 2014-04-04 10:50 有限体积法data_gen
----------- --------- ---------- ----- ----
54505 8
相关资源
- 蓄电池模型 (power battery)
- PSO SVM SVM用于分类时的参数优化
- 流体力学中simple算法的matlab代码
- 最小二乘算法(LMS)处理滤波并预测
- chirp信号的frft的峰值搜索仿真
- matlab编写的蛙跳聚类算法(SFLA)
- LBP特征提取的MATLAB实现
- 数字图像滤波程序
- GCC时延估计
- 基于kaiser窗的双谱线插值FFT谐波分析
- 加窗插值FFT算法的matlab 512点仿真
- GM_PHD_Filter_v110
- pls 可以用于计算排列熵
- 用α-β滤波器结合卡尔曼滤波器对相位
- matlab实现遗传算法 (Matlab Genetic Al
- 切比雪夫(Chebyshev)多项式展开式的
- SIR病毒传播模型的matlab程序
- STATCOM无功补偿不平衡负载,阶梯调制
-
matlab/simuli
nk仿真永磁无刷直流电机 - matlab 开发光伏阵列通用仿真模型
- 激光雷达点云滤波的topscan算法实现
- 比较分析协作通信三个节点时系统 (
- 基于HOUGH变换的车道线检测和识别的程
- matlab画庞加莱(poincare)映射
- 高斯过程回归算法工具箱
- 增广拉格朗日乘子法matlab源码
- 外国人开发的地磁数据处理全套程序
- 求解电力系统经济调度等优化问题的
- matlab频率估计
- 计算脉搏信号的方法 PPG (Blood Puls
评论
共有 条评论