matlab求解复杂隐函数,MATLAB复杂隐函数求解!求大神指导!
问题描述:
MATLAB复杂隐函数求解!求大神指导!
syms D1 D2 G1 G2 al E1 E2 Km w6
G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;
for D2=1.3:0.1:1.8;
ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));
ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)));
w=2*(pi^(1/2))*E1*(G1^(D1-1))*D1*(al^(D1/2))/3;
r=(2/(3-2D1))*(al^((3-2D1)/2)-ac1^((3-2D1)/2))-((al^((2-D2)/2))*2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac1^((5-2D1-D2)/2)));
t=Km*w6*D1*(al^(D1/2))/2;
u=(2/(2-D1)*(ac1^((2-D1)/2))-((al^((2-D2)/2)))*2/(4-D2-D1)*(ac1^((4-D2-D1)/2)));
P2=w*r+t*u;
w1=2*(pi^(1/2))*E2*(G2^(D2-1))*D1*(al^((D1+D2-2)/2))/3;
r1=2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac2^((5-2D1-D2)/2));
t1=Km*w6*D1*(al^((D1+D2-2)/2))/2;
u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));
P1=w1*r1+t1*u1;
P=P1+P2;P=3*10^4;
%想在已知P和D1,D2的情况下求解al,D1和D2是一组数据,最后画出三维图
1个回答
分类:
数学
2014-12-16
问题解答:
我来补答
题目写的比较乱.明确几点:1、是不是说对于D1=1.3:0.1:1.8,D2=1.3:0.1:1.8的每一种组合,按照P=P1+P2=3*10^4来求解a1?2、对于这样的复杂方程,求解析解是不太现实的,应该求数值解.3、方程可能存在多个解,怎么处理?例如,取D1=D2=1.5,画出P-3E4的曲线如下:显然,在0-30范围内有两个解.
再问: 你好,谢谢你的回答。以上函数中的al在我的设定中是有一定物理意义的,应该是个小于1大于0的正数。是否能够给出:在已知P的值时,得出D1,D2,al的三维关系图。麻烦大神把程序传上来!我快被这事弄崩溃了!
再答: 参考代码:syms D1 D2 G1 G2 al E1 E2 Km w6
G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);
[DD1,DD2] = meshgrid(1.3:0.1:1.8);
A1 = DD1*NaN;
for ii=1:size(DD1,1)
for jj=1:size(DD2,2)
D1 = DD1(ii,jj);
D2 = DD2(ii,jj);
ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));
ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1))); (因字数限制,此部分代码不变,删除)
t1=Km*w6*D1*(al^((D1+D2-2)/2))/2;
u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));
P1=w1*r1+t1*u1;
P=P1+P2;
try
A1(ii,jj)=fzero(@(a1)subs(P-3E4,a1),[eps 1]);
catch
fprintf('Error: D1=%g, D2=%g, %s\n',D1,D2,lasterr)
end
end
end
mesh(DD1,DD2,A1)
view(3)
其中,当D1=1.3, D2=1.4时,无解(用ezplot(P-3E4,[0 1])容易看到,曲线始终位于实轴之下)。
再答: 参考代码:syms D1 D2 G1 G2 al E1 E2 Km w6
G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);
[DD1,DD2] = meshgrid(1.3:0.1:1.8);
A1 = DD1*NaN;
for ii=1:size(DD1,1)
for jj=1:size(DD2,2)
D1 = DD1(ii,jj);
D2 = DD2(ii,jj);
ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1))); (此处有部分代码无改动,因字数限制省略)
u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));
P1=w1*r1+t1*u1;
P=P1+P2;
try
A1(ii,jj)=fzero(@(a1)subs(P-3E4,a1),[eps 1]);
catch
fprintf('Error: D1=%g, D2=%g, %s\n',D1,D2,lasterr)
end
end
end
mesh(DD1,DD2,A1)
view(3)
注意,有些情况下(例如D1=1.3, D2=1.4)方程无解。
展开全文阅读
matlab求解复杂隐函数,MATLAB复杂隐函数求解!求大神指导!相关推荐
- matlab object是什么意思,求大神指导一下代码大概是什么意思。。
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function varargout = apple(varargin) % APPLE M-file for apple.fig % APPLE, by ...
- matlab四条曲线围成面,matlab中怎么给四条曲线作出图例啊?求大神指导
clear all: clc disp('插值求取平滑曲线后的数值' ) t = [5 10 15 20 30 40 60 80 120] ; A = [29.84 40.01 46.71 53.87 ...
- matlab 赋空值,未对输出参数赋值 求大神帮忙解惑
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 x=imread('F:/flower.bmp'); wname='db5' [Lo_D,Hi_D] = wfilters(wname); lf=leng ...
- matlab运行一直正忙,matlab一直显示正忙,没反应,求大神帮助
代码如下,运行一直处于正忙状态,求大神帮忙看看代码有什么问题:Q %% 定义 (x,t) 平面上的网格点坐标 clear;clc; dx=0.05; % x 方向的步长 dy ...
- MATLAB中出现 索引超出矩阵维度,老是提示索引超出矩阵维度,求大神指导
function main clear L=32;%基底尺寸的设定: step1=500;%生长层数的设定, niu=1;deltt=0.05;deltx=1;dd=16;aa=(2*dd/deltx ...
- ubuntu14.04扩展屏幕后,打开matlab就死机。求大神帮忙解决!!!万分感谢
ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢
- 求大神帮忙,想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码)
求大神帮忙,想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码) 已知y是一个1*101的数组如下,函数如下,想求得x 135 130 134 128 129 13 ...
- matlab需要的软硬件,软件运行硬件环境 软件著作权-用MATLAB写了个小程序,想申请软件著作权,求大神指点...
您好,请问可以看看您的软件.硬件环境是怎么描述... 硬件环境:就是指你软件运行所需要的的硬件配置: 软件环境:就是指你软件运行所需要的系统的配置. 示例如下: 硬件环境:Intel® Pentium ...
- matlab 连续两个if,求大神解答一个matlab中的for循环嵌套if选择语句
公告: 为响应国家净网行动,部分内容已经删除,感谢读者理解. 话题:求大神解答一个matlab中的for循环嵌套if选择语句回答:1.最后缺少一个end;2.这句输出没什么意义:disp A(i)=A ...
最新文章
- ERROR LNK2019:无法解析的外部的符号 _sscanf或者_vsprintf
- 关于Unity中调试C#的方法
- windows中卸载Jenkins
- Distributed Systems-leader based分布式一致性协议
- numba 让python速度提升百倍
- 【备忘录】Product cost collectors – 成本收集器
- 做任务一定要看测试用例
- JDK+Tomcat+MySql环境配置—linux
- 基于matlab的捕食算法,【优化求解】基于matlab细菌觅食算法的函数优化分析【含Matlab源码 217期】...
- 加密设备攻防(二)- 智能设备篇
- 基于QT实现简单的音乐播放器
- 基于用户的协同过滤推荐算法研究
- 山东科技大学计算机研究生学院官网,山东科技大学-计算机科学与工程学院
- 中国互联网创业工具库Startup Tools
- SaaS启动阶段增长指南(上)
- 斐讯音箱控制扫地机器人_【斐讯 R1 智能音箱使用体验】操作|功能|音质|识别_摘要频道_什么值得买...
- 人生不怕晚,就看敢不敢|优锘科技 X《无尽攀登》专场见面会
- win10, cuda 9.0, python 3.5环境下复现 Flow-Guided Feature Aggregation for Video Object Detection 问题总结
- proposal中文翻译_proposal 是什么意思_proposal 的翻译_音标_读音_用法_例句_爱词霸在线词典...
- 华衫科技-实训课程-小滴服务(Html-Css-Javascript)
热门文章
- 使用Aspose.cells截图并去除截图中水印_含资源下载链接
- 学Python爬虫,不看看m3u8文件如何加密?i春秋 m3u8 文件加密解析
- 【阿里云】企业自建 Hadoop 迁移至 MaxCompute 最佳实践学习
- 考勤不是非得按指纹 刷脸操作更亲民!
- 2022年读过的书单
- python中的第三方库 cdk_用Jython调用CDK(Chemistry Development Kit)显示分子并高亮子结构...
- 在线!在线!在线 !疫情推动传统企业数字化转型!
- 怎么用dwl做个html语言模板,dwl是什么文件
- Excel转PDF最简单的方法
- 运维堡垒机深信服平台操作-安装控件调用客户端工具