如何绘制复杂的隐函数?

20次浏览(最近30天)
Rozana阿拉姆
Rozana阿拉姆 2020年5月10日
评论道: 马修•诺伊 2021年8月30日6点50分
PLmax = A /1+aexp(−b[arctan(h /R)−A)+ 10log(h2 + R2)+ b
在那里,PLmax = 10分贝
一个= 12.0810;
b = 0.1139;
B = 194.2648
一个= -21.4000
h R是隐式的。
图应该是这样的,但是我的代码没有得到任何输出。有人能修改它吗?
clc;
清晰;
f1 = 2000 e + 06;
c = 3 * 10 ^ 8;
itaNLOS = 23分贝;
《= 1.6 db;
y =(4 *π)/ c;
z = log10 (f1);
B1 = (20 * z) + (2 * log10 (y)) + 23
一个= itaLOS-itaNLOS
一个= 12.0810;
b = 0.1139;
P = 10分贝;
f = @ (h, R) (A / 1 + * exp ((- b) * (((h / R)每股)——)))+ 10 * log10 (h。^ 2 + R ^ 2) + B1
fimplicit (f, [500:500:4500 500:500:4500])
5个评论
明星黾
明星黾 2021年8月28日14:30
对于任何想知道情节的人来说……
f1 = 2000 e + 06;
c = 3 * 10 ^ 8;
itaNLOS = 23;
《= 1.6;
B2 = 20 * log10 (f1) + 20 * log10(4 *π/ c) + itaNLOS;
一个= itaLOS-itaNLOS;
一个= 12.0810;
b = 0.1139;
f = @ (h, R) (a / (1 + . * exp ((- b)。* (atand (h / R)))) + 10。* log10 (h。^ 2 + R ^ 2) + b2 - 110);
G =fimplicit(f,[0 3000 0 3000]);
包含(“h”
ylabel (“R”

登录评论。

答案(1)

马修•诺伊
马修•诺伊 2021年5月19日
你好
我做了一些尝试来了解问题出在哪里。我不得不说,我不是这个领域的专家,我不知道是否有更好的方程,但至少我发现了一些小错误,至少 隐式的 函数可以工作,但与预期的P目标值和预期的形状都不一样。所以我想知道是否仍然存在一个问题,无论是由于常量还是我们如何执行方程(我围绕“官方”做了一些变化,但没有人给出令人满意的行为)
所以为了知道对于给定的h,R对方程会返回什么典型的P值,我首先通过生成h和R的网格来绘制等高线图。
你可以看到P在100到130分贝之间,这就解释了为什么 隐式的 如果P的目标是10db,就不会有解了,这就超出了我们的预期,所以,如果P应该是10db,那么我们就有问题了要么是常数,要么是方程的写法;函数求值时也缺少P 隐式的
clc;
clearvars;
f1 = 2000 e + 06;
c = 3 * 10 ^ 8;
itaNLOS = 23;%数据库;
《= 1.6;%数据库;
B1 = 20 * log10 (f1) + 20 * log10(4 *π/ c) + itaNLOS;
一个= itaLOS-itaNLOS;
一个= 12.0810;
b = 0.1139;
P = 10;%数据库;
%%第一次绘图使用网格
x = 50:50:4500;
y = 50:50:4500;
[h R] = meshgrid (x, y);
f = a / (1 + * exp (- b * ((h / R)——每股)))+ 10 * log10 (h。^ 2 + R ^ 2) + B1;
图(1)中,轮廓(x, y, f, 10); colorbar (“绿色”);
colormap(飞机)
集(gca),“YDir”“正常”
%返回
% % % % % % % % %
% P = 10;%数据库;
P = 125;%数据库;
f = @ (h, R) a / (1 + * exp (- b * ((h / R)——每股)))+ 10 * log10 (h。^ 2 + R ^ 2) + B1 - P
f,[500 4500 500 4500]
2的评论
马修•诺伊
马修•诺伊 2021年8月30日6点50分
我的荣幸!

登录评论。

标签

社区寻宝

在MATLAB中心找到宝藏,并发现社区如何可以帮助你!

开始狩猎!