双二次Lagrange 有限元计算特征值程序(基于iFEM)
function lambda = c0P2(h)
%% Mesh
[node,elem] = squarequadmesh([,,,],h);
elem = elem(:,[,,,]);
showmesh(node,elem);
findnode(node);
findquadelem(node,elem);
%% Construct Data Structure
[elem2dof,edge,inDof] = c0dofP2(elem);
elem2dof=double(elem2dof);
N = size(node,); NT = size(elem,);
Ndof = N+NT+size(edge,);
A=sparse(Ndof,Ndof);
B=sparse(Ndof,Ndof);
%% Assemble stiffness matrix
% Since Dphi_i*Dphi_j is quadratic,
% numerical quadrature rule is used here
option.quadorder = ; % default order
[pts,weight] = quadquadpts(option.quadorder);
pts=pts*-;
x=pts(:,);y=pts(:,);
h1=h/;h2=h1;
area=h2*h1;
%% Assemble Matrix
for i=:
for j=i:
DuDv=area*(Dxphi(x,y,h1,i).*Dxphi(x,y,h1,j)...
+Dyphi(x,y,h2,i).*Dyphi(x,y,h2,j))'*weight;
uv=area*(phi(x,y,i).*phi(x,y,j))'*weight;
if i==j
A = A + sparse(elem2dof(:,i),elem2dof(:,j),DuDv,Ndof,Ndof);
B = B + sparse(elem2dof(:,i),elem2dof(:,j),uv,Ndof,Ndof);
else
A = A + sparse([elem2dof(:,i);elem2dof(:,j)],...
[elem2dof(:,j);elem2dof(:,i)],DuDv,Ndof,Ndof);
B = B + sparse([elem2dof(:,i);elem2dof(:,j)],...
[elem2dof(:,j);elem2dof(:,i)],uv,Ndof,Ndof);
end
end
end
%% Solve Ax = lambda Bx, and its first solution is *pi^
lambda=eigs(A(inDof,inDof),B(inDof,inDof),,'sm');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% subfunction Dxphi
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function s = Dxphi(xi,eta,L1,i) % gradient of basis phi for x
switch i
case
s = (eta.*(.*xi - ).*(eta - ))/(*L1);
case
s = (eta.*(.*xi + ).*(eta - ))/(*L1);
case
s = (eta.*(.*xi + ).*(eta + ))/(*L1);
case
s = (eta.*(.*xi - ).*(eta + ))/(*L1);
case
s = -(eta.*xi.*(eta - ))/L1;
case
s = -((eta.^ - ).*(*xi + ))/(*L1);
case
s = -(eta.*xi.*(eta + ))/L1;
case
s = -((eta.^ - ).*(.*xi - ))/(*L1);
case
s = (*xi.*(eta.^ - ))/L1;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% subfunction Dxphi
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function s = Dyphi(xi,eta,L2,i) % gradient of basis phi for y
switch i
case
s = (xi.*(.*eta - ).*(xi - ))/(*L2);
case
s = (xi.*(.*eta - ).*(xi + ))/(*L2);
case
s = (xi.*(.*eta + ).*(xi + ))/(*L2);
case
s = (xi.*(.*eta + ).*(xi - ))/(*L2);
case
s = -((.*eta - ).*(xi.^ - ))/(*L2);
case
s = -(eta.*xi.*(xi + ))/L2;
case
s =-((*eta + ).*(xi.^ - ))/(*L2);
case
s = -(eta.*xi.*(xi - ))/L2;
case
s = (*eta.*(xi.^ - ))/L2;
end
end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% subfunction phi
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function s = phi(xi,eta,i) % gradient of basis phi
switch i
case
s = (eta.*xi.*(eta - ).*(xi - ))/;
case
s = (eta.*xi.*(eta - ).*(xi + ))/;
case
s = (eta.*xi.*(eta + ).*(xi + ))/;
case
s = (eta.*xi.*(eta + ).*(xi - ))/;
case
s = -(eta.*(xi.^ - ).*(eta - ))/;
case
s = -(xi.*(eta.^ - ).*(xi + ))/;
case
s = -(eta.*(xi.^ - ).*(eta + ))/;
case
s = -(xi.*(eta.^ - ).*(xi - ))/;
case
s = (eta.^ - ).*(xi.^ - );
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end
lambda = c0P2(h)
function [elem2dof,edge,inDof] = c0dofP2(elem)
totalEdge=sort([elem(:,[,]);elem(:,[,]);elem(:,[,]);elem(:,[,])],);
[edge, i2, j] = myunique(totalEdge);
N = max(elem(:));
NT = size(elem,);
NE = size(edge,);
elem2edge = reshape(j,NT,);
elem2dof = uint32([elem N+elem2edge (N+NE+:N+NE+NT)']);
i1(j(*NT:-:)) = *NT:-:;
i1 = i1';
bdEdgeIdx = (i1 == i2);
isBdDof = false(N+NE+NT,);
isBdDof(edge(bdEdgeIdx,:)) = true; % nodal
idx = find(bdEdgeIdx);
isBdDof(N+idx) = true;
inDof = find(~isBdDof);
end
[elem2dof,edge,inDof] = c0dofP2(elem)
http://files.cnblogs.com/files/wangshixi12/%E5%8F%8C%E4%BA%8C%E6%AC%A1Lagrange%E6%9C%89%E9%99%90%E5%85%83.rar
双二次Lagrange 有限元计算特征值程序(基于iFEM)的更多相关文章
- (原)使用mkl计算特征值和特征向量
转载请注明出处: http://www.cnblogs.com/darkknightzh/p/5585271.html 参考文档:mkl官方文档 lapack_int LAPACKE_sgeev(in ...
- 第十六届“二十一世纪的计算”学术研讨会 密西根州立大学教授Anil K. Jain主题演讲
Biometrics---How Do I Know Who You Are? 密西根州立大学教授Anil K. Jain主题演讲" title="第十六届"二十一世纪的 ...
- 第十六届“二十一世纪的计算”学术研讨会 图灵奖获得者Butler W. Lampson主题演讲
Personal Control of Digital Data 图灵奖获得者Butler W. Lampson主题演讲" title="第十六届"二十一世纪的计算&qu ...
- 三种计算c#程序运行时间的方法
三种计算c#程序运行时间的方法 第一种: 利用 System.DateTime.Now // example1: System.DateTime.Now method DateTime dt1 = S ...
- 二 ASP.NET MVC 第一个程序 hello world
二 ASP.NET MVC 第一个程序 hello world https://blog.csdn.net/xmroom/article/details/51335917 我使用的Visual s ...
- Senparc.Weixin.MP SDK 微信公众平台开发教程(二十一):在小程序中使用 WebSocket (.NET Core)
本文将介绍如何在 .NET Core 环境下,借助 SignalR 在小程序内使用 WebSocket.关于 WebSocket 和 SignalR 的基础理论知识不在这里展开,已经有足够的参考资料, ...
- 第十六届“二十一世纪的计算”学术研讨会 牛津大学肿瘤成像学教授Michael Brady主题演讲
Computing and Healthcare 牛津大学肿瘤成像学教授Michael Brady主题演讲" title="第十六届"二十一世纪的计算"学术研讨 ...
- Spark练习之通过Spark Streaming实时计算wordcount程序
Spark练习之通过Spark Streaming实时计算wordcount程序 Java版本 Scala版本 pom.xml Java版本 import org.apache.spark.Spark ...
- 微信小程序-基于canvas画画涂鸦
代码地址如下:http://www.demodashi.com/demo/14461.html 一.前期准备工作 软件环境:微信开发者工具 官方下载地址:https://mp.weixin.qq.co ...
随机推荐
- 改写jquery.validate.unobtrusive.js实现气泡提示mvc错误
个人对于这个js.css不是很擅长,所以这个气泡提醒的样式网上找了下,用了这个http://www.cnblogs.com/wifi/articles/2918950.html当中的第一种写法. 对于 ...
- emacs配置eslint 语法检查.找不到node解决
使用emacs配置eslint 当调用语法检查时报错 Suspicious state from syntax checker javascript-eslint: Checker javascrip ...
- 116、Android获取Manifest中<meta-data>元素的值
在AndroidManifest.xml中,<meta-data>元素可以作为子元素, 被包含在<activity>.<application> .<serv ...
- linux内核调试技巧之addr2line
addr2line工具是一个可以将指令的地址和可执行影像转换为文件名,函数名和源代码行数的工具.这在内核执行过程中出现崩溃时,可用于快速定位出出错的位置,进而找出代码的bug. 用法 addr2lin ...
- JavaScript的toString()和valueof()方法
toString()方法: 函数:函数 (function(){}).toString(); //返回"function(){}" typeof((function(){}).to ...
- 编译器zynq设置环境变量
设置临时环境变量export CROSS_COMPILE=arm-xilinx-linux-gnueabi- export PATH=/work/tool/cross_compiler/bin:$PA ...
- Erlang初学
这篇文章主要介绍了Erlang初学:Erlang的一些特点和个人理解总结,本文总结了函数式编程.一切都是常量.轻量进程.进程端口映射及典型缺点等内容,需要的朋友可以参考下 我对 Erlang 编程理念 ...
- mac 10.9 安装 gevent
安装步骤: Gevent依赖libevent和greenlet,需要分别安装. 1,安装 macport (如已安装,可以跳过) 2,通过终端 键入: sudo port install libeve ...
- GHOST(幽灵)重大漏洞
cd /usr/local/srcwget https://webshare.uchicago.edu/orgs/ITServices/itsec/Downloads/GHOST.cgcc GHOST ...
- 查看最点CPU的语句
select c.last_execution_time,c.execution_count,c.total_logical_reads, c.total_logical_writes,c.total ...