| 查看: 366 | 回复: 1 | |||
sunjie2668木虫 (小有名气)
|
[求助]
关于调用函数解方程组(带系数)的问题,求大神帮忙已有1人参与
|
|
%程序为 clc clear l=0.51; f=2; m=0.0055; p=m/(2*l); T=2.695; z=sqrt(p*T); w=0:0.1:10; k=w/(z/p); y=amplitude(f,k,l); y=abs(y); w=w/(2*pi); figure plot(w,y,*k*,*LineWidth*,1.1) %函数为 function y=amplitude(f,k,l) syms f k l A1 A2 B1 B2 C1 C2 D1 D2 M11 M12 M21 M22 M31 M32 M41 M42; [A1,A2,B1,B2,C1,C2,D1,D2,M11,M12,M21,M22,M31,M32,M41,M42] = solve(A1+A2==0,B1+B2==0,A1==M11*exp(-1i*k*l),A2==M12*exp(1i*k*l), ... B1==M21*exp(-1i*k*l),B2==M22*exp(1i*k*l),C1+C2==0,C1==M31*exp(1i*k*l),C2==M32*exp(-1i*k*l),D1+D2==f,D1==M41*exp(1i*k*l),D2==M42*exp(-1i*k*l), ... M11+M12==M21+M22,M21+M22==M31+M32,M31+M32==M41+M42,M12+M22+M41+M31==M11+M21+M42+M32,A1,A2,B1,B2,C1,C2,D1,D2,M11,M12,M21,M22,M31,M32,M41,M42); y= A1*exp(1i*k*l)+A2*exp(-1i*k*l); y = simplify(y); 错误为 Error using plot Conversion to double from sym is not possible. Error in uM_freq_k_l__plot (line 21) plot(w,y,*k*,*LineWidth*,1.1) |
» 猜你喜欢
谈谈两天一夜的“延安行”
已经有15人回复
EST投稿状态问题
已经有6人回复
职称评审没过,求安慰
已经有15人回复
垃圾破二本职称评审标准
已经有11人回复
投稿Elsevier的Neoplasia杂志,到最后选publishing options时页面空白,不能完成投稿
已经有16人回复
毕业后当辅导员了,天天各种学生超烦
已经有4人回复
聘U V热熔胶研究人员
已经有10人回复
求助文献
已经有3人回复
投稿返修后收到这样的回复,还有希望吗
已经有8人回复
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐,对您同样有帮助:
求大神用matlab求下图阴影部分的面积和周长!!
已经有4人回复
现在用MATLAB求解一个循环后画出曲线遇到问题,求大神指点!!
已经有7人回复
关于根据一个系统模型建立复杂微分方程解的问题,一般都是设其解为某种形式。
已经有3人回复
非线性方程的线性化问题,大神们进来指教下!!!!感谢了!!!!!!!!!
已经有9人回复
用MATLAB回归方程的参数,是个隐函数方程
已经有10人回复
求高手解方程组(应该不难)
已经有22人回复
求大神帮忙解题啊!!!
已经有4人回复
有没有大神看一下能不能用FORTRAN求解下面的式子
已经有14人回复
求大神帮忙拟合一个非线性方程,求出模型参数
已经有15人回复
1stopt解方程组
已经有12人回复
关于病态方程组求解的疑问
已经有5人回复
微分方程求解(通解中含hypergeom函数)
已经有17人回复
Mathematica关于绘制隐函数方程组的问题,跪求,快疯了
已经有11人回复
求两个值中最大的值但编写的代码不知道哪里错了,求大神帮忙。
已经有10人回复
求助-帮我用mathematics计算下下面的方程组
已经有13人回复
matlab ode15s方程中两个参数怎么调用其他函数带入求解
已经有6人回复
求助!matlab用fsolve函数求解非线性方程组的问题!
已经有19人回复

【答案】应助回帖
★ ★ ★ ★ ★
感谢参与,应助指数 +1
sunjie2668: 金币+5, ★有帮助, 谢谢,已解决,plot前面加上eval函数就可以了 2014-11-04 23:00:58
感谢参与,应助指数 +1
sunjie2668: 金币+5, ★有帮助, 谢谢,已解决,plot前面加上eval函数就可以了 2014-11-04 23:00:58













回复此楼