| 查看: 371 | 回复: 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) |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神用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人回复
2楼2014-11-04 22:03:13