有限元 弧长法 matlab,有没有有关弧长法的程序!!!
% % % This is a code of Pseudo-Arc Length Continuation Method
% the method can be used for solving the nonlinear equations
% the principle can re found in some text books
clc
clear
% ------------------------开始预测--------------
k=1;
for ds=0.02:0.02:9
a0=1;
x0=3;
y0=2.6458;
a=a0;
x=x0;
y=y0;
step=5;
for i=1:1
DF(1,1)=2*a*x;
DF(1,2)=2*y;
DF(1,3)=x^2;
DF(2,1)=2*x;
DF(2,2)=-2*y;
DF(2,3)=0;
DF;
J1=DF(:,2:end);
J2=DF(:,1);
K=DF(:,3);
J2=[J2,K];
J3=DF(:,1:2);
J1=((-1)^(1+1))*det(J1);
J2=((-1)^(2+1))*det(J2);
J3=((-1)^(3+1))*det(J3);
vv=[J1,J2,J3]';
v=vv/norm(vv);
end
% --------------预测结束-------------------
% -------------开始校正----------------------
x=3;
y=3;
a=2;
p=[x,y,a]';
F(1,1)=a*x^2+y^2-16;
F(1,2)= x^2-y^2-2;
F(1,3)=(x-x0)*v(1,1)+(y-y0)*v(2,1)+(a-a0)*v(3,1)-ds;
delta=1e-6;
epsilon=1e-6;
for i=1:100
J=zeros(3,3);
J(1,1)=2*a*x; J(1,2)=2*y; J(1,3)=x^2;
J(2,1)=2*x; J(2,2)=-2*y; J(2,3)=0;
J(3,1)=v(1,1); J(3,2)=v(2,1); J(3,3)=v(3,1);
q=p-(J\F');
x=q(1);
y=q(2);
a=q(3);
F(1,1)=a*x^2+y^2-16;
F(1,2)= x^2-y^2-2;
F(1,3)=(x-x0)*v(1,1)+(y-y0)*v(2,1)+(a-a0)*v(3,1)-ds;
F;
Z=F;
err=norm(q-p);
relerr=err/(norm(q)+eps);
p=q;
Y=Z;
if (err
q
break
end
end
q
ds;
x0=q(1);
y0=q(2);
a0=q(3);
c(k)=a0;
b(k)=x0;
k=k+1;
a=q(3,1);
end
b
c
plot(c,b)
hold on
% a=c
% x=sqrt(18./(a+1))
% a=2:-0.1:-1
% y=sqrt(18./(a+1)-2)
% plot(a,x,'o')
% hold on
%
a=-1:0.1:1
x=sqrt(18./(a+1))
plot(a,x,'o')
hold on
我找到了这个程序,不知道对不对
有限元 弧长法 matlab,有没有有关弧长法的程序!!!相关推荐
- 图像相减的matlab仿真及光栅滤波法,图像相减的MATLAB 仿真及光栅滤波法实验实现...
图像相减的MATLAB 仿真及光栅滤波法实验实现 毕业设计(论文)中期报告题目图像相减的MATLAB仿真及光栅滤波法实验实现院(系)光电学院专业光信息科学与技术班级090106姓名陈凤学号090106 ...
- 图像相减的matlab仿真及光栅滤波法,图像相减的MATLAB 仿真及光栅滤波法实验实现.doc...
您所在位置:网站首页 > 海量文档  > 计算机 > 图形图像 图像相减的MATLAB 仿真及光栅滤波法实验实现.do ...
- matlab rbf函数_基于径向基函数(RBF)的无网格伪谱法与程序实现(2)——微分矩阵...
参考资料 Gregory E. Fasshauer. Meshfree Approximation Methods with MATLAB. P.387 P401 数值实现 Matlab 2019a ...
- matlab中存档算法代码,MATLAB 智能算法超级学习手册中程序代码
[实例简介] MATLAB 智能算法超级学习手册中程序代码 [实例截图] [核心代码] dc90ef43-7920-434e-bdb8-0636c31c0b44 └── MATLAB 智能算法超级学习 ...
- matlab 坐标不用科学计数法,matlab不用科学计数法
『壹』 matlab中怎么才能不是科学计数法表示结果.比如1.0e+003 * 2.7581,怎么使它显示为2758.1谢谢了,很急啊 format long (小数位14) 或 format sho ...
- 考虑交通网络流量的电动汽车充电站规划matlab 采用matlab软件参照相关资料完成电动汽车程序
考虑交通网络流量的电动汽车充电站规划matlab 采用matlab软件参照相关资料完成电动汽车程序,采用粒子群方法,程序本人编制,运行可靠 ID:5868638495393683快乐程序人
- matlab用jc法计算可靠度,基于MATLAB的截尾分布下JC法计算可靠度
摘 要: 在水工结构可靠度分析中,随机变量的分布形式常因几何尺寸.物理环境等条件限制,传统JC法已经不适用,因此需要对部分变量进行截尾分布处理.在此借助MATLAB丰富的函数资源,编制出截尾分布处理后 ...
- matlab中caitu_tiqu,Matlab平台基于颜色的车牌识别程序
Matlab平台基于颜色的车牌识别程序 matlab 2020-11-18 下载地址 https://www.codedown123.com/51252.html Matlab平台基于颜色的车牌识别程 ...
- matlab画一维波动方程,MATLAB编辑一维波动方程的模拟的程序.doc
MATLAB编辑一维波动方程的模拟的程序 一维波动方程的模拟 程序: function wave_equation() %一维线性齐次波动方程 options={'空间杆长L','空间点数N' ,'时 ...
最新文章
- PLM的关键点—实施篇
- 如何用Python从数据库里面获取数据?4个步骤就能轻松实现
- 0x84bb0001 sqlserver_sqlserver 2000 远程连接 服务器的解决方案
- Vue.js环境搭建
- go结构体初始化_golang中结构体的初始化方法(new方法) | 学习笔记
- DateTime::diff
- 英语不好,也能让你看懂man手册
- 学习笔记(01):Linux操作系统入门培训_Linux/Unix基础培训教程视频课程-01.Linux入门培训之_软硬件平台介绍...
- php 判断浏览器和设备
- 这个双十一,我被淘宝“套娃”游戏整懵了
- spring boot内置容器性能比较(Jetty、Tomcat、Undertow)
- Drug Discov. Today | 药物发现中的先进机器学习技术
- pdf 天线理论与技术 钟顺时_钟顺时. 天线理论与技术[M]. 北京:电子工业出版社, 2011: 290-296....
- SA-LOAM:具有语义辅助的回环检测LOAM系统
- 信用卡诈骗罪16个有效辩点
- 搭建git 服务器出现 Error executing action `create` on resource 'user[GitLab user and group]'
- 解决vue在IE11读取缓存的问题
- jquery.slides.js 幻灯片脚本使用方法
- 高新技术企业是什么?
- 英国旅游--庄园酒店
热门文章
- php+实现群发微信模板消息_使用php实现微信小程序发送模板消息(附代码)
- Excel表格如何输出为清晰图片
- 深度学习入门(三)——神经元激活值的计算方法
- c语言notify头文件,SendNotifyMessage()函数
- bzoj 4819 [Sdoi2017]新生舞会
- Java和JavaScript有什么区别?
- 计算机科学与未解决的难题,NP(未解难题)_百度百科
- (javascript)vue项目打包后,写的覆盖element的样式无效了。这是什么原因?
- mysql怎么连接,mysql怎么连接
- Linux命令如何进入screen,linux screen命令基础