Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2011-05-09, 13:28   #1
likeshu1103
初级会员
 
注册日期: 2011-04-11
年龄: 36
帖子: 8
声望力: 0
likeshu1103 正向着好的方向发展
默认 请问这段代码错误在哪

t=0:10^-12:40000*10^-12; dt=10^(-12);N=40000;
l=1;xs=0.2;
a=2.5*10^(-3);b=9.345*10^(-3);
H=5.76*10^7; %对应文献符号:H=σ
erel=2.5;e0=8.85*10^(-12);u0=4*pi*10^(-7);
L=u0/2*pi*log(b/a);C=(2*pi*e0*erel)/(log(b/a));R=(1/pi*H)*(1/a+1/2*b*0.6*10^(-3));
RL1=100;RL2=10;
x=xs;
G=(1/2*pi)*sqrt(u0/H)*(1/a+1/b); %对应文献符号:G=ξ
v=1/sqrt(L*C);
a1=(G/2)*sqrt(C/L);a2=(R/2)*sqrt(C/L);
for k=0:N
dk=dt*k;
if dk<x/v
g1(k)=0;
else
g1(k)=exp(-a2.*x).*(a1.*x/(2.*sqrt(pi.*(k.*dt-x/v)^3)).*exp(-(a1.*x)^2/(4.*(k.*dt-x/v)));
end
end
有哪位大侠知道的指出来,不胜感激。
likeshu1103 当前离线   回复时引用此帖
旧 2011-05-09, 16:04   #2
yhcode
高级会员
 
注册日期: 2008-05-13
年龄: 44
帖子: 232
声望力: 21
yhcode 正向着好的方向发展
默认 回复: 请问这段代码错误在哪

MATLAB中对于数据或矩阵的第一个元素的下标是1而不是0,你的for循环应该改成从1开始,对于改后其他数据是否也要改动,你自己看吧
__________________
工科'985'博士,(图像处理、分析及理解;模式识别;运动估计;数据分析等)定做程序、算法实现--qq:752105755
yhcode 当前离线   回复时引用此帖
旧 2011-05-10, 19:39   #3
likeshu1103
初级会员
 
注册日期: 2011-04-11
年龄: 36
帖子: 8
声望力: 0
likeshu1103 正向着好的方向发展
默认 回复: 请问这段代码错误在哪

谢谢啦,我再试试吧。
likeshu1103 当前离线   回复时引用此帖
旧 2011-05-10, 19:44   #4
likeshu1103
初级会员
 
注册日期: 2011-04-11
年龄: 36
帖子: 8
声望力: 0
likeshu1103 正向着好的方向发展
默认 回复: 请问这段代码错误在哪

??? g1(k)=exp(-a2.*x).*(a1.*x/(2.*sqrt(pi.*(k.*dt-x/v)^3)).*exp(-(a1.*x)^2/(4.*(k.*dt-x/v)));
|
Error: Incomplete or misformed expression or statement.
请问这是啥意思?
likeshu1103 当前离线   回复时引用此帖
旧 2011-05-10, 21:39   #5
yhcode
高级会员
 
注册日期: 2008-05-13
年龄: 44
帖子: 232
声望力: 21
yhcode 正向着好的方向发展
默认 回复: 请问这段代码错误在哪

问题是:表达式不完整或错误表达
检查括号是否匹配,以及各个变量的size是否匹配。如果不嫌麻烦可以将该表达式拆成多个小的表达式一步一步验证,最后在整合
__________________
工科'985'博士,(图像处理、分析及理解;模式识别;运动估计;数据分析等)定做程序、算法实现--qq:752105755
yhcode 当前离线   回复时引用此帖
旧 2011-05-12, 12:04   #6
风湿百合
初级会员
 
注册日期: 2011-05-12
帖子: 7
声望力: 0
风湿百合 正向着好的方向发展
默认 回复: 请问这段代码错误在哪

您好!请问您是否可以在百忙之中帮我写该一段matlab的程序!十分感激您的帮助!
[r,q]=size(x);
[s2,q]=size(t);
s1=6;
[w1,b1]=rands(s1,r);
[w2,b2]=rands(s2,s1);
df=100; %学习过程显示频率
me=3000; %最大训练步数
eg=0.001;%误差平方和指标
lr=0.1;
lr_inc=1.05;
lr_idec=0.95;
mom_const=0.1;
err_ratit=1.04;
tp=[df me eg lr lr_inc lr_idec mom_const err_ratio];
[w1,b1,w2,b2,epochs,errors]=trainbpx(w1,b1,'logsig',w2,b2,'purelin',x,t,tp);
x=[0.00 0.252 0.439 0.626 0.813 1.000 0.00 0.190 0.392 0.595 0.798 1.000 0.00 0.252 0.439 0.626 0.831 1.000,
1.000 0.795 0.592 0.378 0.204 0.000 1.000 0.815 0.626 0.415 0.189 0.000 1.000 0.835 0.637 0.451 0.235 0.000];
a=simuff(x,w1,b1,'logsig',w2,b2,'purelin');
风湿百合 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 22:31


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.