二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf

二维有限体积法计算热传导及源码

//#include "stdafx.h"

#include

using namespace std;

#define q 500

#define k 1000

void main ()

{ //input the value you want:

double L,dx,dy,T,Ax,Ay,d;

int m,n,i,j,kk,mm ;

//char str1[20];

char *str1="conf.dat";

//cout<>str1;

//string str1("data.");

ifstream fin(str1);

if(!fin)

{

cout<

exit(1);

}

cout<

fin>>L>>m>>n>>T;

cout<

//creat the maxtrix

mm=m*n;

dx=0.3/m;

dy=0.4/n;

vector ae(mm),ap(mm),aw(mm),an(mm),as(mm),su (mm),sp(mm) ;

Ax=L*dy;

Ay=L*dx;

//结¨¢点Ì?0,左Á ¨®下?角?顶£¤点Ì?

aw[0]=0.0;

ae[0]=k*Ax/dx;an[0]=k*Ay/dy;as[0]=0.0;

sp[0]=0.0;su[0]=q*Ax*k;

ap[0]=aw[0]+ae[0]+an[0]+as[0]-sp[0];

for(i=1;i

{

if(i

{

aw[i]=0.0;

ae[i]=k*Ax/dx;

an[i]=k*Ay/dy;

as[i]=k*Ay/dy;

su[i]=q*Ax*k;

sp[i]=0.0;

ap[i]=aw[i]+ae[i]+as[i]+an[i]-sp[i];

}

else if(i==n-1)

//左Á ¨®上¦?角?顶£¤点Ì?

{

as[i]=k*Ay/dy;

ae[i]=k*Ax/dx;

aw[i]=0.0;

an[i]=0.0;

sp[i]=-2*as[i];

su[i]=q*Ax*k+2*as[i]*T;

ap[i]=aw[i]+ae[i]+as[i]+an[i]-sp[i];

}

else if(i==mm-n) //右® ¨°下?角?顶£¤点Ì?

{ ae[i]=0.0;

as[i]=0.0;

aw [i]=k*Ax/dx;

an[i]=k*Ay/dy;

su[i]=0;

sp[i]=0;

ap[i]=aw [i]+ae[i]+as[i]+an[i]-sp [i];

}

else if(i==mm-1) //右® ¨°上¦?角?顶£¤点Ì?

{

ae[i]=0.0;

aw[i]=k*Ax/dx;

as[i]=k*Ay/dy;

an[i]=0.0;

su[i]=2*as[i]*T ;

sp[i]=-2*as[i];

ap[i]=aw

二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf相关推荐

微信小程序：二开版优化新紫色UI云开发新款壁纸小程序源码
这是一款由旧版https://airymz.com/3387.html进行优化二开的一个新版本本壁纸表情包头像小程序采用(dcloud云开发)所以无需服务器与域名无需服务器.无需域名.云开发直接上 ...
二开版优化新紫色UI云开发新款壁纸小程序源码支持用户投稿在线审核
本壁纸表情包头像小程序采用(dcloud云开发)所以无需服务器与域名无需服务器.无需域名.云开发直接上线特点:支持用户投稿,后台审核后会发订阅消息给用户提示作品审核状态,增加用户粘性,支持后端修改 ...
【Matlab条形码识别】二维条形码识别【含GUI源码 607期】
一.代码运行视频(哔哩哔哩) [Matlab条形码识别]二维条形码识别[含GUI源码 607期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 蔡利梅.MAT ...
【Matlab图像拼接】块匹配全景图像拼接【含源码 742期】
一.代码运行视频(哔哩哔哩) [Matlab图像拼接]块匹配全景图像拼接[含源码 742期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 蔡利梅.MATLA ...
【Matlab图像检索】综合特征图像检索【含GUI源码 395期】
一.代码运行视频(哔哩哔哩) [Matlab图像检索]综合特征图像检索[含GUI源码 395期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 蔡利梅.MAT ...
【Matlab肌电信号】肌电信号处理【含GUI源码 966期】
一.代码运行视频(哔哩哔哩) [Matlab肌电信号]肌电信号处理[含GUI源码 966期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] 包子阳,余继周,杨 ...
【Matlab语音处理】声音信号频谱分析仪【含GUI源码 325期】
一.代码运行视频(哔哩哔哩) [Matlab语音处理]声音信号频谱分析仪[含GUI源码 325期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1]韩纪庆,张磊, ...
金字塔图像融matlab,金字塔变换的图像融合Matlab源码.pdf
金字塔变换的图像融合Matlab源码.pdf 附录1 金字塔变换图像融合方法程序 % 拉普拉斯金字塔融合函数 function Y = fuse_lap(M1, M2, zt, ap, mp) % M ...
针对大量数据，MATLAB生成EXCEL文件并进行排版处理的源码
#MATLAB生成EXCEL文件并进行排版处理的源码项目要求:给定某海洋站潮位数据在一年之年的每分钟潮位数据(一共有24×60×365=525600个数),txt格式,要求制作出海洋站该年份的潮汐表 ...

二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf

二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf相关推荐

最新文章

热门文章