matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫

这是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

end，

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...相关推荐

matlab imfill holes,到底怎么实现，bwfill和imfill - 程序语言 - 小木虫 - 学术科研互动社区...
认真研读一下MATLAB的help文档吧,解释最权威: BWFILL Fill background regions in binary image. BWFILL is a grandfathere ...
matlab读取超大txt文件,fortran如何读取超大文本文件 - 程序语言 - 小木虫 - 学术科研互动社区...
有这样一个外部力的文件,需要每个时间步调用 Fbrk.txt 文件内容如下,第一列是时间 1.07800 13.73559 -16.39831 ...
matlab风力机叶片仿真教程,我用matlab编写了一个风力机叶片的优化设计程序，有点问题 - 程序语言 - 小木虫 - 学术科研互动社区...
主程序 clear; global r; global R; global phi; global lamb; V1=7;P=400;lamb0=6;B=3; D=2.5; R=D/2; N=60*l ...
matlab中并行用不了,请教matlab在linux下的并行问题 - 程序语言 - 小木虫 - 学术科研互动社区...
刚在网上搜了一下,matlabpool open local 8 设置节点数,但是必须在matlab下面才行,比如我打开linux, >>matlab 进入matlab系统,然后matla ...
matlab实现kuan滤波,用Matlab设计一个0.5-20HZ的带通滤波器！ - 程序语言 - 小木虫 - 学术科研互动社区...
这是我在网上看到的,希望对你有所帮助. ----------------------------------------------------------------------------- fu ...
c语言调用dll实例结构体传参,Matlab调用Dll，对于自定义结构体数组如何传参？ - 程序语言 - 小木虫 - 学术科研互动社区...
%[例子1] % Matlab调用DLL,整形数组传参 v = [1,22,33;44,-56,67]; pv = libpointer('int16Ptr',v); m = get(pv,'Valu ...
matlab如何取等距点,matlab编程等距取点 - 程序语言 - 小木虫 - 学术科研互动社区...
各位大神,我想用matlab编一个程序,目的是对于一组数据等距取点,但不会写,请各位赐教! 首先,我有一个21×4的矩阵A[21,4],希望对每一列数据每间隔3个点取一个数,然后将这些取出来的点保存在 ...
matlab简单的程序,一段简单的matlab程序 - 程序语言 - 小木虫 - 学术科研互动社区...
原程序是可以运行的,为全面理解程序内容,我将分以下几个部分进行分析: 1."for x=varx"怎么理解? 请参看Matlab关于for函数的帮助文件: Syntax:for i ...
matlab 摄氏度符号怎么打,【求助】弱弱的问下，Matlab如何在公式中输入希腊字母？ - 程序语言 - 小木虫 - 学术科研互动社区...
很多时候都要在matlab画图的时候添加一些公式符号之类的,有一些特殊的字符并不能直接从键盘上输入,比如希腊字母等等.但是有想用,因为这样使图看起来漂亮而且容易理解. 例如:我想输入摄氏度的符号,怎么 ...

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...相关推荐

最新文章

热门文章

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术 科研 互动社区...

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术 科研 互动社区...相关推荐

最新文章

热门文章

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...

matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术科研互动社区...相关推荐