作业帮关于MATLAB方程组绘图
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/20 19:00:27
关于MATLAB方程组绘图
三个方程如下
tan(mq*h)/tan(mp*h)+4*k^2*mp*mq/(mq^2-k^2)^2=0
mp=sqrt(w^2/cl^2-k^2)
mq=sqrt(w^2/ct^2-k^2)
要画出其中k和w的图像
请求知道的人告诉我吧,此题追加100分.
要用数值法就解,其中要用到迭代求根算法.
cl=4000;ct=1970;h=0.00384;
多谢您的回答,但是代入我自己的参数以后.这个图像有些不对,貌似要用数值法来求解.麻烦您再帮我看看好吗?
三个方程如下
tan(mq*h)/tan(mp*h)+4*k^2*mp*mq/(mq^2-k^2)^2=0
mp=sqrt(w^2/cl^2-k^2)
mq=sqrt(w^2/ct^2-k^2)
要画出其中k和w的图像
请求知道的人告诉我吧,此题追加100分.
要用数值法就解,其中要用到迭代求根算法.
cl=4000;ct=1970;h=0.00384;
多谢您的回答,但是代入我自己的参数以后.这个图像有些不对,貌似要用数值法来求解.麻烦您再帮我看看好吗?
再检查一下公式
clc;clear;
cl=4000;ct=1970;h=0.00384;
mp=sym('sqrt(w^2/cl^2-k^2)')
mq=sym('sqrt(w^2/ct^2-k^2)')
mp=subs(mp)
mq=subs(mq)
f=sym('tan(mq*h)/tan(mp*h)+4*k^2*mp*mq/(mq^2-k^2)^2')
f1=subs(f)
ezplot(f1,[0 0.001 0 0.0000001])
运行结果：
f1 =
tan(12/3125*(1/3880900*w^2-k^2)^(1/2))/tan(12/3125*(1/16000000*w^2-k^2)^(1/2))+4*k^2*(1/16000000*w^2-k^2)^(1/2)*(1/3880900*w^2-k^2)^(1/2)/(1/3880900*w^2-2*k^2)^2
clc;clear;
cl=4000;ct=1970;h=0.00384;
mp=sym('sqrt(w^2/cl^2-k^2)')
mq=sym('sqrt(w^2/ct^2-k^2)')
mp=subs(mp)
mq=subs(mq)
f=sym('tan(mq*h)/tan(mp*h)+4*k^2*mp*mq/(mq^2-k^2)^2')
f1=subs(f)
ezplot(f1,[0 0.001 0 0.0000001])
运行结果：
f1 =
tan(12/3125*(1/3880900*w^2-k^2)^(1/2))/tan(12/3125*(1/16000000*w^2-k^2)^(1/2))+4*k^2*(1/16000000*w^2-k^2)^(1/2)*(1/3880900*w^2-k^2)^(1/2)/(1/3880900*w^2-2*k^2)^2