matlab雷诺曲线,关于用有限差分法解雷诺方程

这是我编写的程序，希望对你有用处！

function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot)

Fx=0;

Fy=0;

epsilon=0.6;

L=0.03;

d=0.06;

lambda=L/d;

m=40;

n=30;

phi1=0;

phi2=2*pi;

delta_phi=(phi2-phi1)/m;

delta_lambda=2/n;

k=1;

P0=zeros(m+1,n+1);

A=zeros(m+1,n+1);

B=zeros(m+1,n+1);

C=zeros(m+1,n+1);

D=zeros(m+1,n+1);

E=zeros(m+1,n+1);

F=zeros(m+1,n+1);

H=zeros(m+1,n+1);

for i=1:1:m+1

% theta(i)=(i-1)*delta_phi;

for j=1:1:n+1

H(i,j)=1+epsilon*cos((i-1)*delta_phi);

end

S=0;

T=0;

ERR=1e-3;

GAP=1;

while GAP>ERR

k=k+1

for i=1:1:m+1

for j=1:1:n+1

if (i==1)|(j==1)|(i==m+1)|(j==n+1)

P(i,j)=0;

else

A(i,j)=(1+epsilon*cos((i+1/2-1)*delta_phi))^3;

B(i,j)=(1+epsilon*cos((i-1/2-1)*delta_phi))^3;

C(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;

D(i,j)=(d/L)^2*(delta_phi/delta_lambda)^2*(1+epsilon*cos((i-1)*delta_phi))^3;

E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j);

F(i,j)=6*delta_phi*((1+epsilon*cos((i+1/2-1)*delta_phi))-(1+epsilon*cos((i-1/2-1)*delta_phi)))+12*(delta_phi)^2*(X_dot*cos(phi1+(i-1)*delta_phi)+Y_dot*sin(phi1+(i-1)*delta_phi));

P(i,j)=(A(i,j)*P0(i+1,j)+B(i,j)*P0(i-1,j)+C(i,j)*P0(i,j+1)+D(i,j)*P0(i,j-1)-F(i,j))/E(i,j);

if P(i,j)<0

P(i,j)=0;

else

end

% for i=2:1:m

% for j=2:1:n

S=sum(sum(abs(P-P0)));

T=sum(sum(abs(P)));

% end

GAP=S/T;

P0=P;

end

for i=1:1:m+1

for j=1:1:n+1

Fx=Fx+P(i,j)*cos(phi1+(i-1)*delta_phi)*delta_phi*delta_lambda;

Fy=Fy+P(i,j)*sin(phi1+(i-1)*delta_phi)*delta_phi*delta_lambda;

end

matlab雷诺曲线,关于用有限差分法解雷诺方程相关推荐

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...
这是MATLAB有限差分法求解雷诺方程,可以参考下.希望对你有用 function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot) Fx=0; Fy=0; epsilon= ...
matlab解薛定谔方程,有限差分法解薛定谔方程与MATLAB实现
第30卷第3期高师理科学刊 Vol. 30 No.3 2010年 5月 Journal of Science of Teachers′ College and University Ma ...
matlab做拉普拉斯曲线,Matlab 拉普拉斯滤波以及 imfilter 详解（多图）
Matlab 拉普拉斯滤波以及 imfilter 详解(多图) 08-08栏目:技术 TAG:imfilter imfilter 代码: jhua.orgclc clear all A=imread ...
数值方法3：偏微分方程1：使用有限差分法解一维热传导（扩散）方程
初学练习,看b站课程,教学为matlab代码,自己写的Python代码,后面会放b站课程链接,感兴趣的同学可以学习观看. 说明:Python初学者,代码可能不够漂亮,欢迎大家批评指正.本系列代码用no ...
matlab将求解sin隐式解,Matlab隐式符号方程求解和赋值
近日处理了一个隐式方程的求解,由于方程含有较多的未知数,而且这些参数均是跟实验相关的一些参数,所以,必须得到需要求解的解与这些参数之间的一个表达式.之前是考虑用的Maple推导求解了该隐私方程,求解 ...
车牌识别与计算机编程,基于MATLAB的车牌识别程序详解.ppt
基于MATLAB的车牌识别程序详解自定义一个字符函数,用来从车牌区域中提取出7个字符,其中利用切割函数来进行切割. 程序:function [word,result]=getword(d) word ...
matlab 剪毛刺,MATLAB消除曲线毛刺Outlier Detection and Removal [hampel]
MATLAB消除曲线毛刺Outlier Detection and Removal [hampel] function [YY, I, Y0, LB, UB, ADX, NO]=hampel(X, Y ...
MATLAB解隐函数方程时符号表达式转化为数值的方法-用vpa函数
今天在解决一个小问题时,遇到解隐函数方程,中间涉及一个解的传递问题,才好好研究了一下 syms这个语句的一些语法规则,最终用vpa这个函数解决了符号表达式到double数值的转化. syms是符号函数 ...
雷诺手表怎么换电池？雷诺手表更换电池的方法
雷诺手表怎么换电池?雷诺手表更换电池的方法现在很多人都戴机械手表,石英表或者电子表,无论是把控时间还是作为饰品非常流行,而每块手边的电池寿命也就一两年,有的表本身不是很昂贵去修表店有些值不当的,但是 ...

matlab雷诺曲线,关于用有限差分法解雷诺方程

matlab雷诺曲线,关于用有限差分法解雷诺方程相关推荐

最新文章

热门文章