| 查看: 1753 | 回复: 7 | ||
[求助]
现在用MATLAB求解一个循环后画出曲线遇到问题，求大神指点！！ 已有3人参与
|
|
现在用MATLAB求解一个循环后画出曲线遇到一个小问题； 源程序如下： Dj=6.93; k=0.28; t=1/(2*k); angle=45.6*pi/180; d=16; a=angle; b=pi/2; n=10; h=(b-a)/n; for i=1:1:n+1; psaiI=angle+(i-1)*h; z=fsolve(@(z)tan(psaiI)*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))-1/((2*z-1)^t)-sqrt((2*z-1)/z^2),1); Dm=Dj*sqrt(z^2/(2*z-1)); xx=fsolve(@(x)cos(x)/sin(psaiI-x)-sqrt(z^2/(2*z-1)),0); plot(psaiI,xx) %绘制曲线，psaiI为横坐标，xx为纵坐标 end 为什么运行总是错误！感觉前面计算都可以，就是在plot(psaiI,xx)这步出问题，求大神指点啊。对matlab不熟悉啊。 |
回复此楼» 猜你喜欢
化工京区271求调剂
已经有4人回复
-
295求调剂。一志愿报考郑州大学化学工艺学硕，总分295分
已经有5人回复
-
求调剂
已经有10人回复
-
材料284求调剂，一志愿郑州大学英一数二专硕
已经有11人回复
-
0856化工专硕求调剂
已经有15人回复
-
303求调剂
已经有5人回复
-
289求调剂
已经有4人回复
-
江苏省农科院招调剂1名
已经有4人回复
-
290分材料工程085601求调剂 数二英一
已经有3人回复
-
0856调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解非线性方程组报错，请各位大神指点！
已经有8人回复
- matlab画图问题！！！急！希望会的大神来看看！！
已经有11人回复
- 求高手解方程组（应该不难）
已经有22人回复
- 如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！
已经有11人回复
- Tecplot新手求各位大神指导
已经有6人回复
- matlab为啥出现以下这种错误
已经有6人回复
- matlab画图问题: 怎样用小图标或自定义图形代替+ x . o等做图
已经有10人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- matlab的plot画图（如何画出与取值数目不一样的标记符）
已经有4人回复
- matlab解决重心法选址问题，运行不了啊，求助各位大神
已经有10人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- matlab可以进行曲线插值吗？
已经有3人回复
- 用matlab绘曲线，但总是出现问题
已经有5人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 如何用matlab画出两组实验测定数据概率密度曲线，求两条曲线重叠部分面积
已经有8人回复
- 求教一个origin中的问题，求大神指点
已经有4人回复
- 【求助】如何用Matlab求联合概率曲线
已经有3人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
aazzyyhello
铜虫 (初入文坛)
- 应助: 5 (幼儿园)
- 金币: 74.6
- 红花: 1
- 帖子: 43
- 在线: 93.8小时
- 虫号: 3332448
- 注册: 2014-07-21
- 专业: 电磁场与电路
815292578
2楼2014-10-23 07:57:49
3楼2014-10-23 08:08:10
aazzyyhello
铜虫 (初入文坛)
- 应助: 5 (幼儿园)
- 金币: 74.6
- 红花: 1
- 帖子: 43
- 在线: 93.8小时
- 虫号: 3332448
- 注册: 2014-07-21
- 专业: 电磁场与电路
4楼2014-10-23 08:11:33
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6335.8
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 206.6小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
【答案】应助回帖
感谢参与，应助指数 +1
|
Dj=6.93; k=0.28; t=1/(2*k); angle=45.6*pi/180; d=16; a=angle; b=pi/2; n=10; h=(b-a)/n; for i=1:1:n+1; psaiI(i)=angle+(i-1)*h; z=fsolve(@(z)tan(psaiI)*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))-1/((2*z-1)^t)-sqrt((2*z-1)/z^2),1); Dm=Dj*sqrt(z^2/(2*z-1)); xx(i)=fsolve(@(x)cos(x)/sin(psaiI-x)-sqrt(z^2/(2*z-1)),0); end plot(psaiI,xx) %绘制曲线，psaiI为横坐标，xx为纵坐标 |
» 本帖已获得的红花（最新10朵）
5楼2014-10-23 08:42:10
6楼2014-10-23 16:55:08
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
815292578: 金币+10, ★★★★★最佳答案, 十分感谢啊..... 2014-10-27 10:39:13
感谢参与，应助指数 +1
815292578: 金币+10, ★★★★★最佳答案, 十分感谢啊..... 2014-10-27 10:39:13
|
[code] z=fsolve(@(z)tan(psaiI)*sqrt((1-(2*z-1)/z^2)*(1-1/((2*z-1)^t)^2))-1/((2*z-1)^t)-sqrt((2*z-1)/z^2),1); xx(i)=fsolve(@(x)cos(x)/sin(psaiI-x)-sqrt(z^2/(2*z-1)),0); [\code] 这两句中的psaiI改为psaiI(i) |
7楼2014-10-24 10:20:09
8楼2014-10-27 10:38:55
