| 查看: 1825 | 回复: 0 | ||
[求助]
求画出matlab中用surf程序画出的两曲面的交线的程序
|
|
clc clear all; d=0.046;C_1=0.002;t=0.4;p_1=7870;p_2=865;u=0.00003;b=0.0001;l=5;A_2=1; w=0:0.5:25; z_n=0:2.5/75:2.5; Q=zeros(51,76); for i=1:51 for j=1:76 Q1(i,j)=(2*b.*exp((3.*z_n(j))./(2.*l) - (t.*w(i))./2))/(w(i)+ (6.*d.^3.*u)./(C_1.^3.*((3.*p_2)./5000 + p_1.*(d - 3./100).^2))) - (2.*b.*w(i).*exp((3.*z_n(j))./(2.*l) - (t.*w(i))./2))./(w(i)+ (6.*d.^3.*u)./(C_1.^3.*((3.*p_2)./5000 + p_1.*(d - 3./100).^2))).^2 - (b.*t.*w(i).*exp((3.*z_n(j))./(2.*l) - (t.*w(i))./2))/(w(i) + (6.*d.^3.*u)./(C_1.^3.*((3.*p_2)./5000 + p_1.*(d - 3./100).^2))); Q2(i,j)=(3.*A_2.*b.*w(i).*exp((3.*z_n(j))./(2.*l) - (t.*w(i))./2))./(l.*(w(i)+ (6.*d.^3.*u)./(C_1.^3.*((3.*p_2)./5000 + p_1.*(d - 3./100).^2)))); end end [z_nz_nw,www]=meshgrid(z_n,w); subplot(1,2,1) surf(Q1); view(50,40) grid on xlabel('z_n'); ylabel('w'); zlabel('a/A_2'); hold on; surf(Q2); view(50,40) grid on xlabel('z_n'); ylabel('w'); zlabel('a/A_2'); samedot=(abs(Q1(i,j)-Q2(i,j))<=0.5); z_nz_n=samedot.*z_n; ww=samedot.*w; zz=samedot.*Q1(i.j); plot3(z_nz_n(samedot~=0),ww(samedot~=0),zz(samedot~=0),'r-'); subplot(1,2,2) h1=plot3(z_nz_n(samedot~=0),ww(samedot~=0),zz(samedot~=0),'r-'); |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有35人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有15人回复
-
基金中了
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
应用matlab求两个曲面的交线
已经有7人回复
- 如何在surfer里画出两个面的交线？
已经有4人回复
- 【求助完成】matlab画图surf
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50