３つの変数の関係が、グラフが表示されない（枠外にも外れるし）

前提

plot3Dでグラフが表示されない

実現したいこと

３つの変数の関係（ｆとｒとｐ）を３次元グラフで表示させたい

発生している問題・エラーメッセージ

エラーメッセージ![イメージ説明](https://ddjkaamml8q8x.cloudfront.net/questions/2022-12-30/11c15619-699b-471d-ba2e-92060ae1aeb5.jpeg)
エラーは無いですが、各点のドットもなく、横にはみ出たような、添付なグラフ表示で。

該当のソースコード

library(rgl)   

# Sound Pressure of Ultrasound from the cylindrical small diaphragm

r <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)
f <- c(1,2,3,4,5,6,7)
p <- matrix(c(1:105),nrow=15,ncol=7)

a = 0.001
k =  4187 
gamma = 3.1415/6
z<- k*a*sin(gamma)
J1 = 0

# for (m in 1:100)
for (m in 1:100)
{
  fact = 1
  for (n in 1:m)
    fact= fact*n
  G = 1
  for (n in 1:(m+1))
    G = G*n
 
  Q = (((-1)^m)/(fact*G))*((z/2)^(2*m+1))
  cat("m & fact & G & z & Q & J1 = ",m,fact,G,z,Q,J1,"\n")
  J1 = J1 + Q
}
  
Vm = 12

for (h in 1:15)
{
  for (i in 1:7)
  {
    l = r[h]*0.01
    omega = 2*3.1415*10^f[i]
    p[h,i] = (((a^2)*Vm*omega*l)/l)*abs(J1/z)
  }
}

plot3d(p,r,f)

試したこと

plot3dの中のpとｒとｆの順番を入れ替えたりもしましたが、変わらず

補足情報（FW/ツールのバージョンなど）

RStudioです

melian

2022/12/30 05:14

何だかよく分かりませんが、 r <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) としているので、以下のループで l が 0 になり、結果、p が全て NaN になっています。 for (h in 1:15) { 　for (i in 1:7) 　{ 　　l = r[h]*0.01 　　omega = 2*3.1415*10^f[i] 　　p[h,i] = (((a^2)*Vm*omega*l)/l)*abs(J1/z) 　} }

pontoon

2022/12/30 08:22

コメントありがとうございます。チェックしてみます。

行動規範の内容に同意します

回答1件

library(rgl)   

# Sound Pressure of Ultrasound from the cylindrical small diaphragm

r <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15)
f <- c(1,2,3,4,5,6,7)
p <- matrix(c(1:105),nrow=15,ncol=7)
dB <- matrix(c(1:105),nrow=15,ncol=7)

a = 0.001
k =  4187 
gamma = 3.1415/6
z<- k*a*sin(gamma)
J1 = 0

for (m in 1:100)
{
  fact = 1
  for (n in 1:m)
    fact= fact*n
  G = 1
  for (n in 1:(m+1))
    G = G*n
 
  Q = (((-1)^m)/(fact*G))*((z/2)^(2*m+1))
  cat("m & fact & G & z & Q & J1 = ",m,fact,G,z,Q,J1,"\n")
  J1 = J1 + Q
}
  
Vm = 12

for (h in 1:15)
{
  for (i in 1:7)
  {
    rou = 1.04
    l = r[h]*0.01
    omega = 2*3.1415*10^f[i]
    p[h,i] = (((a^2)*Vm*omega*rou)/l)*abs(J1/z)
    dB[h,i] = log10(p[h,i]/0.001)
    cat("**** p h i r[h] l f[i] = ",p[h,i],h,i,r[h],l,f[i],"\n")
    cat("\n")
  }
}

print("*** J1 =")
print(J1)
print("*** p  =")
print(p)
print("*** dB  =")
print(dB)
plot3d(p,r,f)

と、プログラムを変えたら、少なくとも、添付のような、まともらしいものが
表示され出しましたが、これは、３次元を俯瞰する角度は変えられないのでしょうか？
（上側からでなく、もう少し、横方向から見たい）

それと、
plot(dat$r,dat$dB,xlim=c(0,15),ylim=c(0,10),col="green")
と、２次元表示でもやってみたのですが、この時は何も表示されません。
何故だろう。