| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 877 人关注过本帖
标题:[求助]急!关于积分问题!
只看楼主 加入收藏
help_me
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2006-6-5
收藏
 问题点数:0 回复次数:3 
[求助]急!关于积分问题!

最近在作论文的时候,碰到一个积分,由于积分式很长,好像还是不可积的,用了好几种方法也没有积出来,希望大家能帮我看看,赐教一下,小弟先谢谢大家了!程序在下面
-----------------------------


Q=10; % 冷却率
T=980; %初始温度
Ae3=870;
time=1; % 时间步
T=T-Q*time; %温度的改变


% 此处开始if语句
if T<=Ae3
Tc=Ae3-T; % 过冷度

syms T; %声明符号变量

% 计算Dr=Dr0*exp(-Qr/(R*T))的参数
Dr0=1.75/10^5;
Qr=143.32;
R=8.3144;

%K1,K2
K1=2.07/10^5;
K2=6.33/10^15;

%k是Boltzmann常数(玻尔兹曼)
k=1.38/10^23;

%超元素S的自由能中的相关参数值
Tmsi=0;
Tnsi=-3;
Tmmn=-39.5;
Tnmn=-37.5;
%Xsi,Xmn为置换型合金元素si,mn的摩尔分数
Xsi=0.2;
Xmn=0.3;

%超元素S的活度系数的相关参数值
%Wr为奥氏体中碳原子间的交换能
Wr=1250;

%Xc是奥氏体中碳元素的初始摩尔分数
Xc=0.5;


%计算:超元素S的活度系数as
Zr=14-12*exp(Wr/(R*T));

%奥氏体到铁素体相变中超元素S的自由能
Gfe=20853.06-466.35*T-0.046304*T^2+71.147*T*log(T);
Gs=141*(Xsi*(Tmsi-Tnsi)+Xmn*(Tmmn-Tnmn))+Gfe;


%铁素体形核的驱动力
Gn=Gs-R*T*(log((1-Zr*Xc)/(1-Xc)))/(Zr-1);

%奥氏体中碳溶质的扩散率
Dr=Dr0*exp(-Qr/(R*T));

% I是在某一温度T时奥氏体晶粒表面单位面积的铁素体形核率
I=K1*(1/sqrt(k*T))*Dr*exp(-K2/(k*T*(Gn)^2));

%f为被铁素体占据的奥氏体晶粒表面分数
f=0;

%积分函数
ft=(I*(1-f))/Q;

%此处是一个关于ft的积分函数,积分限是(0,Tc),积分变量是(Ae3-T),该如何对 ft 函数进行积分呢?????

%Sr 为该元胞的表面积
%Lca,两个近邻元胞的距离。
Lca=20;

%此处元胞类型取为四边形元胞(正方体)
L=Lca;
Sr=Lca*Lca;

nuc_pf=Np*Sr;

end

------------------------------------------

麻烦大家帮我看一下了。谢谢了

搜索更多相关主题的帖子: 积分 变量 论文 温度 
2006-06-05 15:35
help_me
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2006-6-5
收藏
得分:0 

顶!helpme。thanks
2006-06-06 09:50
abingchem
Rank: 6Rank: 6
等 级:贵宾
威 望:24
帖 子:716
专家分:0
注 册:2004-12-30
收藏
得分:0 

Matlab中有几个积分的函数,trapz,quad,quadl,非常好用,不过要编程,你的程序就要大动一下,把ft写成一个function,这样调用就方便了。


此人已消失
2006-06-06 11:30
help_me
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2006-6-5
收藏
得分:0 
谢谢大家了。问题已经解决了
2006-06-06 19:21
快速回复:[求助]急!关于积分问题!
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.030104 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved