matlab用solve函数解方程时不能得到数值解

5641456

用matlab中solve函数解一个方程时给出了很多限制条件，并不能得到数值解，请问怎么解决？谢谢啦代码如下，就是公式里除了alpha其余全是已知量，为了输入一个delta的值解出一个alpha的值。

delta=10;
n_o=1.528;
n_e=1.7297;
lambda=632.8e-9;
d=30e-6;

alpha=solve('tan(alpha)=((n_e^2-n_o^2)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2))^0.5)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))^0.5))/(n_o^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))+n_e^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2)))','alpha')；

TouAkira

1. alpha = vpasolve([tan(alpha)==((n_e^2-n_o^2)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2))^0.5)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))^0.5))/(n_o^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))+n_e^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2)))],[alpha])

复制代码

alpha = 0.083898634335274686247749002355

5641456

TouAkira 发表于 2019-9-15 21:41
alpha = 0.083898634335274686247749002355

您好，首先非常感谢您的回答，我复制了您写的代码发现了如下错误，请问是什么原因导致的？
错误提示：
错误使用 alpha
输出参数太多。

chuyanwu2012

5641456 发表于 2019-9-15 21:50
您好，首先非常感谢您的回答，我复制了您写的代码发现了如下错误，请问是什么原因导致的？
错误提示：
错 ...

试一下这个：

clear;clc
syms alpha
delta=10;
n_o=1.528;
n_e=1.7297;
lambda=632.8e-9;
d=30e-6;

d = vpasolve(tan(alpha)==((n_e^2-n_o^2)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2))^0.5)*(((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))^0.5))/(n_o^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_e^2)/(n_o^2-n_e^2))+n_e^2*((((2*pi*d*n_o*n_e)/((2*pi*d*n_o+delta*lambda)*cos(alpha)))^2-n_o^2)/(n_e^2-n_o^2))),alpha)

5641456

chuyanwu2012 发表于 2019-9-16 08:26
试一下这个：

clear;clc

非常感谢