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